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A  description  of  rolling  contact  failure  modes  is  given 
and  the  variables  affecting  the  life  of  a  rolling  contact  are 
identified,  A  mathematical  mode)  of  subsurface  and  surface 
crack  propagation  is  presented.  The  life  to  failure  of  volume 
elements  in  the  vieinity  of  a  defect  is  formulated.  A  term 
"severity”  of  a  :.ii cr odef ec t  has  been  defined.  The  model  is 
characterized  by  the  inclusion  of  bulk  material  parameters, 
defect  characteristics  and  parameters  of  geometry,  stress, 
lubrication  and  surface  topography.  A  statistical  expression 
for  the  life  of  on  entire  rolling  body  is  based  on  the  defect 
life  formula.  The  new  model  includes  current  standard  bearing 
life  prediction  formulas  os  a  special  case.  To  assist  in  inter¬ 
pretation  of  the  model,  the  stressed  volume  in  a  Hertzian  ellip¬ 
tical  stress  field  has  been  determined  from  the  computed  contours 
of  equal  reversing  shear  stress.  A  stress  analysis  has  been 
conducted  on  the  stresses  near  interacting  asperities  and  around 
a  surface  defect  (furrow). 
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EVALUATION 


1.  Tha  present  technique  for  predicting  the  life  of  a  group  of  rolling  element 
hearings  does  not  consider  the  many  factors  known  to  affect  bearing  Ilf.'  ,  and 
for  large  bearings  the  technique  is  clearly  inadequate .  This  contract  ie  the 
first  part  of  a  two  part  effort  to  develop  a  practical  engineering  tool  for  the 
determination  of  the  expected  life  of  any  group  of  bearings.  This  fir6t  contract 
wee  to  consider  all  the  variables  that  affect  bearing  life  and  the  possible 
failure  mechanisms  Involved,  and  then  to  develop  equations  which  would  contain 
parameters  to  account  for  those  variables  known  to  affect  bearing  life.  These 
objectives  have  been  accomplished  and  the  results  of  this  contract  have  pro¬ 
vided  a  number  of  equations  containing  paranssters  characterising  material, 
geometry,  load,  defect  severity,  and  environmental  variables. 

2.  The  above  mentioned  equations  contain  constants  end  function  signs  which 
swat  be  evaluated  and  determinad  from  test  and  field  data  before  the  technique 
can  be  used  as  an  engineering  tool.  This  la  to  be  accomplished  in  the  second 
effort.  The  results  of  these  efforts  will  be  included  in  &  mechanical  reli¬ 
ability  handbook  and  should  provide  an  improved  prediction  technique . 
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SECTION  I 


ffffSOOCCTIOW  AND  SVMMASY 


This  is  the  Final  iepori  issued  in  fulfillment  of  go^e 
Air  Development  Center  Contract  No,  F30602-6T-C=0i47  on 
''Development  of  Set  hera'niicnl  Models  Predicting  Life  of  Large 
Roller  B  arings". 

The  objective  of  this  work  is  to  determine  the  variable? 
which  cause  the  life  of  (large)  rolling  bearings  to  vary  froa 
the  life  predicted  by  the  existing  desian  methods;  to  determine 
the  effects  of  these  variables  on  the  life  of  rolling  bearings 
and  to  formulate  an  improved  mathematical  model  for  the  predic¬ 
tion  of  rolling  bearing  life. 

The  present  study  covers  the  first  year  of  this  effort  and 
has  led  to  a  general  model  of  bearing  failure.  A  subsequent 
effort  is  currently  underlay  (&ADC  Contract  No.  F3O6O2-60-C-O147 ) 
to  cover  future  development  of  the  model  and  to  make  the  formulae 
sufficiently  specific  for  engineering  use. 

This  report  is  divided  in  several  sections  summarized  as 
follows; 

Section  _II  is  a  description  of  the  principal  concepts  devel¬ 
oped  in  this  Contract.  Using  the  currently  accepted  formula  as  a 
startinq  point,  this  section  brings  together  all  the  nes  concepts 
generated  in  this  study  and  explores  the  usefulness  of  the  new 
model.  An  outlook  on  future  research  is  given. 

Section  III  presents  a  synopsis  of  the  currently  accepted 
Lundberg-Falmgren  bearing  fatigue  life  theory  which  forms  the 
basis  of  the  ASA  standard  for  bearing  rating  and  is  the  starting 
point  for  the  present  study.  This  section  is  included  in  recog¬ 
nition  of  the  feet  that  the  fundamental  work  of  Lundberg  and 
Palragren  may  not  be  easily  accessible  to  all  readers. 

Section  IV  extracts  from  recent  literature,  the  principles 
of  fatigue  failure  theory  required  for  this  study,  as  follows; 

1.  A  Survey  of  rolling  contact  failure  identities  among 
which  Is  spalling  failure.  This  failure  mode  is  the 
subject  of  the  present  study, 

2.  A  listing  is  given  of  variables  affecting  contact 
fatigue  life.  These  variables  are  grouped  into  four 
main  categories,  viz.  material  variables,  surface 
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si crogeaaetry  varishleg,  design  variables  and  operating 
vsrl tbles, 

3.  kn  evaluation  |g  preiesited  of  the  interdependence  of 
the  variables  and  their  effect  on  subsurface  tad  surface 
initiated  apallisg  cecu rr eaee a . 

4.  A  model  is  offered  of  fatigue  failures.  Subsurface 

and  surface  initiated  fatigue  failures,  are  distinguished 
which  compete  to  prosote  spalling  failures  in  rolling 
contact.  In  both  subsurface  and  surface  failures,  the 
fatigue  process  is  described  as  a  sequence  of  phases  of 
fatigue  crack  generation,  propagation,  and  final  fracture, 
(i.e.,  spalling)  at  a  "most  critical*  c^aek  in  the 
rolling  element,  Crack  generation  in  rolling  contact 
is  postulated  to  result  fro®  localized  plastic  strain 
concentration  ground  stress  raisers,. 

Section  V  crvers  the  formulation  of  an  expression  for  the 
crack  growth  rat  as  a  function  of  strength  and  stress  parameters, 
ductility,  and  p  istic  micro-strain. 

This  concept  is  applied  to  a  situation  where  defects  of 
known  ^severity*  exist  is  a  unifera  saafrix,  to  yield  8  foretula 
for  the  fatigue  life  of  a  defect. 

Section  VI  gives  a  statistical  theory  of  failure  for  an 
entire  rolling  body,  based  upon  the  relationship  between  life 
and  defect  severity  developed  in  Section  V.  In  this  node!,  a 
rolling  body  is  conceived  a§  being  built  up  of  a  large  nusfaer 
of  ssall  cells  each  of  which  contains  exactly  one  defect 
(including  "defects”  of  no  influence  at  all).  The  severity  of 
&he  defect  present  in  any  cell  is  a  randos  variable.  The 
distribution  of  life  over  identically  loestedcells  in  a 
population  of  rolling  bodies  similarly  sada  and  operated,  is 
found  through  a  traasforaation  of  distributions  between  severity 
and  life. 

The  distribution  of  rolling  body  life  is  expressed  ss  a 
compound  of  the  individual  cell  life  distributions.  The  asyap- 
totic  distribution  of  shortest  cell  life  is  derived  for  the  case 
where  all  cell  lives  are  coeeoaly  distributed. 

Section  ¥11  covers  a  required  stress  analysis  in  a  Hertzian 
contact.  The  aicro-straia  r ange  in  the  highly  stressed  volume 
(or  surface)  i@  calculated. 


Seefejga  fill  presents  a  specialised  elastic  stress  analysis 
for  determinating  ef  the  aisisue  shear  strati  near  as  idealized 
surface  asperity. 

Section  IK  pretests  g  stress  analysis  for  the  detersi nation 
ef  nasciaus  shea?  stress  beneath  e  f ur?©e=sh aped  surface  defect.. 

Analytical  details  are  ©applied  in  seyersl  AppeedAees. 


. . I . . ■*)  i.i . . . . . sunt . . . 1 . . . r 


! 

i 


SECTION  II 
PHI  NCI  PAL  CONCEPTS 


A  mathematics'  model  of  rolling  contact  fatigue  is  a  complex 
subject.  Numerous  aspect#  of  raa t al 1 u rg i c §1 ,  mechanical,  and 
statistical  nature  have  to  be  considered  in  its  development.  At 
the  present  stage  of  the  development,  many  of  these  details  are 
still  open.  Others  have  been  covered  in  quite  some  depth.  The 
piesent  report  is  a  Summary  of  the  studies  to  date.  It  will  take 
up  the  several  aspects  of  the  problem  in  turn,  at  whatever  depth 
is  currently  accessible.  There  is  the  possibility  with  this  pre¬ 
sentation  that  the  reader  may  be  diverted  from  the  principal 
underlying  concepts  by  the  complexity  of  detail.  In  order  to 
prevent  this  and  to  facilitate  evaluation  of  work  accomplished 
from  the  point  of  view  of  the  engineer,  who  will  ultimately  use 
the  theory  for  practical  life  predictions,  a  review  of  the  principal 
concepts  is  offered  in  the  present  Section.  No  proofs  or  references 
will  be  cited:  these  are  either  given  in  the  subsequent  Sections  or 
referenced  there. 


1  . 


THE  PRINCIPLE  OF  ROLLING  BEARING  LIFE  PREDICTION 


Rolling  bearing  life  is  defined  here  as  fatigue  life.  Causes 
of  failure  other  than  fatigue  are  considered  avoidable  and  are, 
therefore,  eliminated  from  life  prediction.  Fatigue  li  s  is  pre¬ 
dicted  on  the  basis  of  a  cumulative  damage  concept,  i.e.  that  with 
repeated  application  of  cyclic  stresses,  irreversible  material 
changes  take  place  which  ultimately  result  in  failure.  This  con¬ 
cept,  with  its  statistical  implications,  was  first  applied  to 
rolling  contact  life  prediction  by  Lundberg  and  Palagren.  The 
Lundberg-Palagren  concept,  universally  used  today,  revolves  around 
a  phenomenological  equation  of  the  following  form,  between  numbers 
of  cycles  to  failure,  and  macro. jcopic  mechanical  variables: 


log 


S(N) 


N  <P(To.*o>  V 


(2.1) 


where  S(N)  =  the  probability  of  survival  to  N  cycles 
t0  =  maximum  shear  stress  amplitude 
Zo  =  depth  eo-ordinate  of  t0 
V  =  stressed  volume 
tp  =  function  sign 
e  ■-  constant 


b 


This  equities  explicitly  eostilsi  the  number  of  cycles,  a 
moxiBUB  sheas  stress  aad  iti  depth  ee-ordiaate,  isd  the  Ritreifed 

vei uasr”  which  latte?,  howeseir,  is  sever  expressed  is  absolute 

teres,  oily  as  a  facts?  of  praportioial ity. 

Tie  Leadberg-Pglagrea  life  predietion  theory  eossista  of  the 
application  of  Equities  (2.1)  to  the  required  side  variety  e£  geo¬ 
metrical  sad  klaeasiie  eoBditioas  which  characterise  a  ceaplex 
aaseably  §ueh  gg  a  rolling  besriag. 

Physically,  Equities  (2.1)  teaches  that  the  cumulative  prob¬ 
ability  of  eurvivgl  decreases  with  ineregslag  number  j f  cyeles  N, 
sssd  with  isere.gsing  Use  of  the  rolling  contact  system  (stressed 
volume).  Tie  specific  choice  of  the  fuaetios  <p  ( v0  ,  )  mas  ease 

by  Laadberg  asd  Palagreu  oftci  gad  for  §11,  and  is  gives  la 

Equation  (2.2) 

qi  (to  ,  2<j  )  s  Tee  2fl'=h  (2.2) 


e,fe;  eeastaats  >  0 

This  squattest  states  that  the  survival  probability  decreases 
with  iitcre@§ls@  shear  stress  reage,  but  lacrosses  fer  greater 
values  of  the  depth  co-ordinate  of  the  ssxi@ua  shear  stress  reage. 

By  applying  elastic  analysis  to  the  contact  situa-iea, 

Lsadberg  aad  Palegron  derived  detailed  steteaeatg  reggidiag  the 
effect  of  the  pertiseht  aa@r©~geoaetry  parameters  iaflueaeing 
Hortaal  surface  pressures  is  the  easiest,  the  subsurface  shear 
streissf  resulting  f res  these  pressures  aad  the  kiseaatie  para- 
m stars  defersioiBg  ausbers  of  eyeles  ia  terms  of  bearing  ring 
revolutions.  Iqueties  (2.1)  is  readily  modified  to  take  account 
©?  fise-variable  or  §pase~variii?l©  losdiag  l?y  using  the  "Palsgrea- 
iitier  hypothesis”  ©£  damage  ©eeuaul aiioa,  statiag  that  fatigue 
daaege  aeeuaulates  at  a  rat©  depesdiag  only  on  lead  conditions  at 
the  current  tiae,  so  that  Equation  (2.1)  esa  be  srittea  in  the  fora 

log  -X-  ~  /  [J  (p  (  t0  ,  so)  dH]S  dV  (2.3) 
S(N)  (VXN) 

The  pheBoseBologi eal  nature  of  Equatioas  (2.1)  through  (2.3) 

results  ia  aa  i(sp§§§s  if  ese  attempts  te  incorporate  iale  life 
predietiea,  aessly  acquired  knowledge  regardiag  the  effect  sf 
parameters  other  than  eoateet  geometry  gad  kinematics,  since 
these  equations  offer  ao  clue  as  to  the  proper  role  of  such  para¬ 
meters  ia  (lefisisg  life.  Per  this  reggon,  past  attempts  st 


refining  the  Lundbery-Palmgren  theory  have  relied  on  the  fact 
that  Equations  (2.1)  and  (2.3)  are  proportionalities,  l.t.  they 
contain  a  freely  available  constant  multiplier  relating  the 
absolute  magnitude  of  life  to  the  probability  of  survival.  This 
multiplier  is  intended  by  Luadberg  and  Palmgres  as  the  material 
constant,  but  has,  from  time  to  time,  been  used  to  in  orporate  a 
variety  of  correction  factors. 

The  limi tationg of  this  approach  are  obvious,  end  It  has 

therefore  begs  decided  in  the  present  study  of  improved  life  pre¬ 
diction  methods  to  abandon  attempts  at  modifying  the  basic 
Lundberg-Palmgren  equation.  Bather,  it  was  decided  to  generate 
novel  equations  from  which  the  Lundberg-Palmgren  equations  can  be 
obtained  as  a  special  esse. 

The  new  equations  were  derived  using  a  more  detailed  physical 
model  of  fatigue  failure  rather  than  as  purely  phenomenological 
equations.  It  is  recognized  that  the  use  of  sueh  a  model  has 
many  pitfalls,  the  most  obvious  being  that  its  details  may  not  be 
verifiable.  However,  the  drawbacks  are  more  than  compensated  by 
the  heuristic  value  of  a  detailed  model  and  can  be  rendered  harm¬ 
less  by  insisting  that  aoa-veri f i abl e  details  of  the  model  should 
not  enter  into  the  final  engineering  formula  for  life  prediction. 

Before  leaving  this  brief  review  of  Lundberg-P almgren  theory, 
it  is  noted  that  Equation  (2.1)  is  equivalent  to: 

H(N)  =  I  -  S(N)  =  1  -  exp  t-(|*)e}  (2.4) 

where  =  constant  "scale  parameter"  of  the  life 
dis  tribution 

H(PJ)  =  cumulative  probability  of  f a i lure 
within  N  cycles 

i.e.  failures  are  distributed  according  to  a  Weibull  distribution 
with  zero  lower  bound,  characteristic  life  Ns  and  dispersion 
exponent  e.  This  distribution  appears  in  the  Lundberg-Palmgren 
formulation  as  a  result  of  deliberate  choice,  as  a  useful  distri¬ 
bution  for  the  description  of  fatigue  phenomena,  and  its 
appearance  does  sot  stem  from  extreme  value  considerations.  This 
point  will  be  of  interest  later. 
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PRINCIPAL  VABI ABLgS  OF  FATIGUE  FAILURE 


Section  IV  gives  a  detgilgd  regies  e£  the  variables  govern¬ 
ing  a  fatigue  failure  situgtioa.  There,  it  is  deduced  that  tbs 
variables  fill  is  fear  mils  categories:  safaris!  variables,  sur¬ 
face  ai eroggosetry  vsrigbles,  design  variables,  aad  opesrgtiag 
variables. 

Material  variables  are  those  influencing  the  "strength"  of 
the  rolling  system.  Current  fatigue  investigations  (chiefly  of 
the  nos-rollisg  type)  consider  yield  strength  aad  ductility  of 
the  material  as  daaiatat  bulk  (or  matrix)  strength  variables. 
Modifying  these  are  residual  stresses  aad  work  hardening  effects, 
acquired,  in  part,  during  the  course  of  fatigue  life.  1b  rolling 
contact,  the  faaterials  used  are  of  high  hardness,  aad  such 
materials  do  not  react  with  their  matrix  strength.  Rather,  rolling 
contact  life  appears  to  be  determined  by  the  strength  of  the 
material  in  the  vicinity  of  inevitable  eatarial  ieperf ©ctioas 
such  as  inclusions  in  the  matrix  or  microeraeks .  Thus,  the  nature 
of  these  imperfections  is  a  dominant  variable. 

Surface  ai eregeeastry  variables  determine  the  detailed  nature 
of  the  contact  through  mh i eh  loads  are  transmitted  to  the  material. 
The  generalised  roughness  of  the  surface  determines  the  topography 
of  the  eoBiaetiag  surface  eleaeats,  tb©  plastie  behavior  of  the 
material  immediately  adjacent  to  the  surface,  aad  interacts  with 
lubrication,  as  will  be  soon  presently. 

There  are  localised  imperfections  on  surfaces,  mostly  in  the 
fern  of  sharp  depressions  ( "£ arrows” )  which  S ora  stress  raisers 
near  their  §dg©§  aad  are  iafluautial  in  failure.  Of  course,  there 
caa  be  aiay  ether  typsa  e£  surface  variables,  soa@  of  them  arti¬ 
ficial  as  iadueed  by  eoatisg,  special  treatment  of  the  surface,  etc. 

The  cJeiiga  variables  of  the  rolling  contact  ire  dealt  with 
©stescively  by  tie  Laadberg=Palg@r@s  theory,  aad  they  are,  there- 
fore,  quite  familiar.  Track  length,  conformity  between  rolling 
elements,  dimensioBs  sad  nueber  of  these  elements,  contaet  angles, 
parameters  aaflaiag  the  precise  cross  track  geometry  gad  easy 
others  are  influential,  primarily  because  they  determine  the 
macroscopic  (Hg?f*ig@)  stress  field  sad  the  somber  of  cycles  as  a 
function  of  bear leg  ring  rotation.  They  also  determine  the 
magnitude  of  the  highly  stressed  volume. 


Operating  variables  escseptsi  the  external  ssvlrosaent  if? 
which  the  rolling  costaet  system  must  oadure,  iseladisg  lead, 
speed,  lubrication,  tespertture,  sad  atHotpherie  eoaditieai.  The 
influeace  of  load  is  deteraiaiag  the  stress  field  is  obvious,  and 
so  is  that  of  ipeed  in  determining  the  number  of  cycles  per  unit 
time.  However,  these  two  variables  alio  interact  with  the  lubri¬ 
cant  to  determine  the  faydrody  eaai  es  of  the  of  ten-present  pressure- 
bearing  el  as tehydredyoaal e  (EHD)  lubricant  fils  }i  the  eostaet, 
which  redistributes  stresses  sad  has  iaportant  effects  @n  the 
mi crobehavi o r  of  the  contact  area  (asperity  interaction).  Teeper- 
ature  esters  by  influencing  both  the  rasleriel  strength  end  the 
lubrication,  sad  atmospheric  conditions  can  be  of  consequence  if 
they  influence  lubricant  behavior  or  cause  corrosive  effects. 

This  list  of  influential  paraseters  deters  the  theorist  by 
iti  multiplicity.  The  only  practicable  approach  to  the  develop¬ 
ment  of  a  life  prediction  formula  in  the  presence  of  such  a 
multitude  of  variables  is  to  find  a  flexible,  simple  concept 
describing  failure  mechanise,  and  then  solve  the  problem  of  intro¬ 
ducing  each  variable  by  defining  its  impact  on  that  seehaaiss. 
Success  of  such  an  atteapt  depends  on  the  proper  choice  of  the 
mechanism  and  will  necessarily  be  Halted.  There  will  always  be 
variables  that  the  model  cannot  secouaodate,  and  as  time  ptoses, 
the  influence  of  these  will  become  more  and  more  recognized,  lead¬ 
ing  ultimately  to  the  abandonment  of  the  model.  However,  the 
model  will  serve  well  in  the  isteris  if  it  permits  account  to  be 
taken  of  the  most  important  parameters  recognised  to  date.  In  what 
follows,  such  a  model  will  be  outlined.  It  appears  at  the  present 
stage  of  the  study  to  have  the  required  flexibility  sad  to  accom¬ 
modate  many  of  the  most  important  parameters,  including  all  those 
which  the  Luadberg-Palagrea  theory  utilises.  It  will  remain  for 
further  study  to  develop  the  specific  formulations  for  the  incor¬ 
poration  of  new  parameters  into  this  model  and  to  shoe  whether  it 
is  sufficiently  free  free  inherent  contradictions  to  be  practically 
used.  This  further  effort  is  currently  underway. 

3,  THE  FATIGUE  FAILUHE  SOBIL 

Our  failure  model  visualizes  fatigue  damage  as  the  growth  of 
8  crack.  There  are  plestic  flew  occurrence!,  carbon  migratioa, 
and,  of  course,  first  of  all,  dislocation  motions  la  the  notrix  as 
a  result  of  cyclic  stressing,  whieh  precede  or  are  concurrent  with 
crock  formation.  However,  for  modeling  purposes,  these  subtler 
occurrences  are  not  helpful  in  fixing  ideas  of  fatigue  deegge 
because  a@  way  is  known  t©  measure  the  degree  to  which  they  are 
shortening  life  expectancy.  The  effeet  of  a  crack  on  life  is 
intuitively  clear:  when  the  crack  has  become  large  enough  a  piece 


of  material  will  separate  from  the  surface,  forslug  the  fatigue 
spall  defined  is  rolling  contact  technology  as  fatigue  failure. 
Creek  growth  is  visualised  ©s  i rrevers ibl e,  so  that  e  suitable 
measure  of  crack  size  satisfies  the  concept  of  "damage"  as 
irreversible  progress  towards  failure.  The  fatigue  phenomenon 
s tarts , accordingly, wi tb  the  initiation  of  a  crack,  and  proceeds 
through  stages  of  its  prow.h  until  the  craek  has  become  large 
enough  to  form  a  spa'l.  Fatigue  life,  as  determined  by  the  crack 
in  question,  begins  iith  the  onset  of  cycling  and  terminates  when 
the  spall  forms.  Of  course,  it  can  be  argued  that  a  rolling  con¬ 
tact  system  say  be  functional  in  the  presence  of  a  spall  of 
tolerable  siae.  There  is  room  to  accommodate  this  argument  in 
the  sodai.as  will  be  pointed  out  later,  but  the  current  descrip¬ 
tion  of  the  failure  terminates  with  first  spalling. 

It  is  coavgBigal  to  describe  failure  generation  in  three 

phiaes : 


Phase  I  begins  immediately  upon  the  oasetof  cyclic  stressing, 
a  si  i  s  consumed  by  the  formation  and  growth  of  a  aicrocrach.  A 
definition  of  a  sicracraeh  will  be  given  below.  From  the  point 
of  view  of  the  model,  it  is  characterized  by  the  fact  that  it  i§ 
small  enough  aot  to  interact  with  other  microereeks  that  may  be 
present  in  the  rolling  element. 

In  Phase  II,  the  craefe  grows  aaeroscopically  until,  at  the 
end  of  this  pEiase,  it  has  reached  a  critical  sise,  defined  by 
the  fact  that  it  is  now  large  enough  to  cause  the  initiation  of 
precipitous  crack  growth  (in  Phase  III). 

Phase  III  is  occupied  by  precipitous  eraek  growth  at  a  rate 
greatly  is  excess  of  Phase  II  growth.  This  precipitous  eraekisg 
fores  the  spall  itself.  This  Phase  @@y  be  virtually  1 BstSBt§@e©B§ 
or  eeususe  substantial  length  of  time,  depending  on  whether  one 
specifies  a  siniauB  spall  glge,  which  is  aeeepted  as  a  failure, and 
depeading  on  a  variety  of  material  sed  operating  conditions. 

The  objective  of  a  mathematic©!  life  formula  is  to  describe 
crack  growth  through  the  above  three  phases.  In  order  to 
describe  crack  growth,  one  selects  s  measure  of  crack  size  A 
and  formulates  @s  equation  of  the  form: 

A  =  / 1  <N,  X,  )  (2.5) 


where  A  s  creek  §ig§ 

/i  s  fuHftisn  sign 
X£  =  undefined  pgreseterg 
W  s  number  of  stress  cycles 


10 


At  the  present  state  of  the  study,  it  seems  begt  to  select  a 
iorni  of  Equation  (2.5)  common  in  current  fatigue  theory,  vis.  one 
defining  the  first  derivative  of  crack  size  (the  crack  growth 
rate)  in  terrag  of  the  relevant  variables: 

M  =  /S(N,  A.  Xt  )  (2.6) 

dN 

where  /_  =  function  Sign 

Life  prediction  is  then  accomplished  by  determining  from  Equation 

(2.6)  ti»gt  value  of  the  number  of  cycles  which  corresponds  to 

a  critical  crack  size  Ac,  causing  immediate  spalling,  i.e. 

\  =  fa(Ac,  xt  )  (2.7) 

where  =  life  at  failure 

Ac  =  critical  crack  size  at  spalling 

Many  current  theories  of  fatigue  failure  use  the  crack 
growth  Equation  (2.6)  in  the  following  simple  form: 

H  =  A(e e ,  D)  a  A  (ec  (N).  D(N)  )  (2.8) 

where  A  =  function  sign 

ee  =  Bl aa  ti c  strain  at  the  propagating  craek  front 
D  =  ductil i ty 

According  to  Equation  (2.8),  the  only  variables  entering  the 
crack  propagation  equation  area  plastic  strain  (measured  at  the 
propagating  craek  front)  and  a  measure  of  material  duetility. 
Specific  definitions  of  these  two  variables  in  terms  of  measurable 
physical  quantities  are  open  at  this  point,  both  beeau.se  appro¬ 
priate  definitions  for  the  rolling  contact  situation  have  not 
previously  been  determined,  and  also  because  the  plastic  strain 
ec  is  a  mi cropsraraeter  which  is  r.ot  directly  measurable.  Note 
that  Equation  (2.8),  for  all  its  simplicity,  contains  many 
assumptions.  Only  the  firgt  derivative  of  crack  size  appears 
explicitly.  (The  crack  size  itself  enters  by  way  of  its  influence 
on  ec.)  Only  variables  measurable  at  the  location  of  the  propa¬ 
gating  craek  front  appear,  and  these  only  with  their  values 
assumed  @t  the  time  of  the  N-ih  stress  eycle,  (Of  course,  the 
equation  is  compatible  with  a  dependence  ec(N)  and  D(N),  and  the 
specific  form  of  this  dependence  deternines  whether  this  formula 
satisfies  the  Palmgren-Miner  hypothesis.)  There  is  hardly  room 
for  concern  about  res tri c ti venes s  at  this  point,  however,  since 
even  Equation  (2.8)  is  much  too  general  to  be  practically 
@ppl i cabl e . 
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f/Ggr  f,e  add  definition  to  Equstioa  (2.8),  the  concent  of 
or  at  tL  I  ^St?0duced“  h  d©f§et  is  a  location  in  the  aaterigl 

"  *  it  Infr*  i**  t&er@  lf  a  ^a6aaeT  of  crack  genera 

EouaHoJ  (9  a  ll  §3§Ua9a  thSt  tH§  tendeBe?  •nnltaat,  itself  ia 
.(?h8)^bJ  §SB®  property  ef  the  defect  causing  er  to  be 
higher  at  the  defect  than  elsewhere  ia  the  vicinity.  C 

J f®r"  thls  relationship  can  be  writtea  by  aggutB)Ba 
*h! M  ?  r  eC£ /“*  10  UBl§*lal  tension)  ,  critical  Lflaft 
eouid  "e!j,e#jie  strs*n  fi@id  ■■  ** 

yield**  *\d*PSn(il  eRl^  BB  •«.  on  *»tUt  aeJerUy^ad'o/tfr'’ 

yield  strength  ay  of  the  satrix,  i.@. 


sc  =  ?  ' ®»  se#  0  y ) 


(2.9) 


where  e, 


§o  = 


plastie  strain  at  defect 

critical  " undisturbed"  total  (elastic 

plus  pl@g?ic)  strain 

defect  severity  eeagura 

(micro)  yield  strength  of  the  matrix 
fuactionsign 


The  defect  severity  factor  9  is  defined  as  s  „  .  .  „ 

factor  ch  iracteri stic  of  the  defect.  Conveniently  8  i!  dsflJ 
for  all  real  defects  with  strain  veniently,  &  is  defined 

as  "iMfTnTil*  d.flet"  *  esi8iB®  properties,  but  else  for 

a  inutseetive  defect  with  no  severity  at  all.  i  *  nae,  whioh 

doet  sot  rsiae  the  oggnitode  of  the  strata  t  »  & 

®  ,  1  ia&  strain.  For  purposes  of 

ilf!  !  f  ‘  £  ls  ceS9®aie°t  to  define  g„eb  "ineffective" 

defects  the  Halting  egse  of  defects  ^ith  rAn.  e#  iB©^eetiv@ 
properties.  oe^ects  with  real  stress  raising 

letreduei ag  Equation  (2.9)  into  Equation  (2.0),  one  has 


an 


h  (Ot  §0  «  CTy ,  D ) 


(2.10) 


era dl-lLdllVtV  1°  8ep@r§te  tfse  ^riableg  influential  in 
*  ftf  U  inte  the  two  groups;  variable  related  to  defects 

riiTdi:  Vd- @  18  ;*•  varubie 

Tnlldndl  Slmplid7'  the  o^feets^ar^eoasol  idated 

into  a  single  function  y,  i.e. 


dH 


-  A  (i,  y) ; 


Y  =  Y  (e0 ,  of!  D) 


(2.11) 


The  a@n  assumptions  underlying  Equations  (2.9)  through  (2.11) 
are  that  there  is  a  scalar  e0  of  the  macros trai ■  field  which, 
alone,  among  strain  field  character! §ti e§ ,  determines  ec,  sad 
?bs£  all  setrix  parameters  exert  their  iafluene©  ©b  crack  gregrtfa 
via  a  single  quantity  y.  Neither  of  these  assumptions  is  esses- 
tial  in  order  to  arrive  at  a  workable  sedel,  but  ire  made  here 
in  the  absence  of  a  sore  refined  understanding  of  the  actual 
physical  situation,  to  arrive  at  a  relatively  siaple  foraulaiios. 

Further  development  of  Equation  (2.11)  requires  use  of  a 
further  restrictive  concept  such  as  the  hypothesis  of  multipli¬ 
cative  effectg  on  fatigue  life.  This  concept,  also  used  by 
Lundberg  and  Palsgren,  asserts  that  the  rate  of  fatigue  damage 
(crack  growth)  eea  be  expressed  as  a  product  ef  a  aaaber  of 
independent  factors,  i.e. 

=  fl  (  9iai  )  (2.12) 

where  cp  =  unspecified  independent  faetors 
II  =  multiplication  operator 
at  =  constants 

Applying  this  concept  to  Equation  (2.11),  one  say,  by  suit¬ 
able  definition  of  the  functions  @  and  y,  absorb  in  thee  the 
function  A,  and  write 

~r  =  @  •  Y  (2.13) 

uN 

He  will  proceed  now  to  the  examination  of  the  aatrix  factor 
Y  and  the  defect  factor  6. 

4.  THE  MATRIX  FACTOR  Y 

From  the  definition  of  y  given  in  Equation  (2.11),  it  is  a 
function  of  a  critical  total  macrostrain  scalar  at  the  location 
of  the  growing  craek  front,  and  of  a  yield  strength  and  a  ductility 
measure. 

The  yield  strength  measure  ay  will  be  a  mlcroyielri  stress, 
since  small  scale  plastic  occurrence;  are  at  issue  in  rolling 
contact  fatigue.  It  must  be  taken  with  its  value  at  the  time 
of  the  N-tb  stress  cycle,  to  aeeount  for  work  hardening  er  @erk 
softening  of  the  matrix. 
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The  ductility  is  defiled  is  stgtie  teiill©  tests  as  the  reduc¬ 
tion  1b  erogg  sectional  irea  at  fracture.  It  Is  not  obvious  that 
this  simple  definition  will  apply  under  the  conditions  of  ceatact 
fatigue,  bat  ii  is  intuitively  convincing  that  eae  should  isclad© 
is  the  formula  a  material  ductility  property  measuring  the  amount 
of  plastic  strain  the  satris  can  absorb  before  it  erseks.  A 
variety  ef  metallurgical  parameters,  but  also  seae  operatisg 
eonditioag,  sill  determine  ductility.  Tie  sost  Important  eperat- 
lag  condition  is  hydroitstic  coapressive  strass.  It  is  generally 
believed  that  the  high  hydrostatic  compression  eoepeneat  esiitiag 
in  most  of  the  Hertzian  contact  stress  field  retards  crack  ferae- 
tion.  This  fact  #111  aanifest  itself  1b  a  point-wise  varying 
value  of  the  ductility  parameter  D  whea  esesiniag  material 
sleaests  located  at  different  poiats  within  the  Hertzian  stress 
field.  Inasmuch  as  it  gay  depead  on  work  hardening,  ductility  ea@ 
ba  s  function  of  ft.  Thus,  the  ductility  parameter  is  already  Snows 
to  depend  on  saterial  eoastaats,  on  a  parameter  of  the  stress 
field,  and  can  depead  on  ft.  It  say  also  be  related  to  other  oper¬ 
ational  parameters,  e.g.  to  eyeliag  rate.  These  relationships  are 
gyabolized  by  the  following  equation: 


0  ~  °(M,  oh,  N,  q>) 


(2. Id) 


where  U  -  saterial  variables 

=  hydrostatic  compressive  stress 
<p  =  operetieg  factors 

Turning  to  the  critical  aacrastrain  parameter  $0 ,  it  is 
obviously  dependent  os  the  variables  of  load  I,  ceatact 
geoeetry  p,  position  under  the  eoateet  x,  and  elastic  modulus  E, 
deflaiag  between  them  the  elastic  Hertz  stress  field.  With  a 
quasi-elastic  assueptioa,  these  parameters  give  a  relationship 
of  the  fors 


e0  =  Sq  CM,  E.  K,  p) 


(2. IS) 


where  I  =  load 

E  -  elastic  modulus 
1"  =  position  vector 
p  =  contact  geometry  parameters 

The  Bquasi-el§stie”  assusptios  operates  os  the  scale  of  the 
whale  Hertsian  stress  field  and  disregards  the  vicinity  of  defects. 
It  postulates  that  the  aeeroieepic  total  strain  ea  egg  he  calcu¬ 
lated  from  the  elastic  (F.eriilan)  stress  field  of  the  rolling 
coat act.  This  is  the  cese  if  the  leads  are  sufficiently  la#  that, 
at  west,  very  saill  iBeunts  of  aaeroseopic  plastic  flaw  take 
plaea  so  that  the  plastic  coapoeeat  sf  total  strata  ig  negligible 
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and  that  plastic  flow  does  not  result  in  a  significant  redistribu¬ 
tion  of  elastic  stresses.  Except  for  the  generation  of  residual 
stresses  due  to  cyclic  stressing,  this  is  a  reasonable  assumption 
in  all  practical  rolling  contact  fgtigue  situations.  The  question 
of  residual  stresses  will  require  separate  ex aei  n  a  t  i  o  n  .  Generally, 
they  are  handled  by  assuming  that,  after  a  small  ausber  ef  cycles, 
the  residue!  stresses  hare  "shaken  dawn"  to  §  constant  value. 

Then,  they  act  as  a  superimposed  static  stress  field  and  '=  e»b  i  ft  * 
with  the  cyclic  stresses.  The  resulting  tine-variable  stress 
field  is,  of  course,  different  from  that  existing  without  residual 
stresses,  and  assumptions  must  be  made  regarding  the  effect  of 
this  difference  on  crack  propagation. 

A  common  assumption  in  fatigue  theory  is  that  superposition 
of  a  static  stress  field  does  not  alter  the  plastic  strain  e 
influencing  crack  growth  rate.  However,  the  hydrostatic  compres¬ 
sion  coaponent  of  the  residual  stress  field  may  aodify  the 
dueti 1  ity  D. 

Everything  said  above  about  quasi-elastic  behavior  is  festrieed 
to  the  matrix  at  locations  remote  from  defects.  Due  to  the 
stress  raising  effect  of  defects,  it  is,  of  course,  possible  that 
localised  plastic  occurrences  take  place  in  saal 1  volumes  in 
their  vicinity.  Such  ai cropl as t i ci ty  is,  in  fact,  the  condition 
of  cracking  in  the  proposed  model. 

For  any  given  defect  ©,  and  given  matrix  strength  (oy),  it  is 
possible  to  delineate  that  portion  of  the  Hertzian  stress  field 
within  which  the  aacrostrain  e0  is  high  enough  to  cause  plastic 
Bicrostrain  ec.  For  a  given  population  of  defects,  there  will  be 
a  "realistic"  maximum  severity  0.  One  can  delineate  a  highly 
stressed  areo  in  the  Hertz  stress  field  within  which  all  micro- 
plastic  occurrences  occasioned  by  defects  of  "realistic”  severity 
will  be  confined.  This  definition  of  a  "highly  stressed  sons"  will 
be  adopted  in  what  follows,  and  the  cress  sectional  area  of  this 
highly  stressed  seas,  in  a  plane  perpendicular  to  the  rolling 
direction,  will  be  designated  by  s, 

S.  THE  DEFECT  SEVERITY  FACTOR  0 

The  effect  of  §  "defect"  is  generating  plastie  strain  ia  its 
vicinity  is  manifestly  very  complex,  A  simple  relationship  for  S 
will  be  proposed  as  follows: 

9(N)  =  0  (d,  A(K),  $)  (2.16) 
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she  re  d  -  initial  defe-'t  severity 

k  =  A(N),  crack  else  after  N  cycles 
S  =  gise  of  the  highly  stressed  z&ns 


The  variables  of  0  are:  d ,  ths  severity  a£  the  dated  it  the 
onset  efeyelie  stressing;  A(N),  the  crack  else  at  the  tiae  of  the 
N-th  stress  cycle;  and  S,  the  cross  sectional  area  of  ths  highly 
stressed  gone  in  the  Hertz  contact. 

As  mentioned  before,  crack  propagation  can  be  envisioned  as 
occurring  ia  three  phases.  The  first  phase,  aierocrack  propaga¬ 
tion,  can  be  defined  using  Equation  (2.16)  by  postulating  thst 
the  creek  is  so  small  that,  by  comparison,  the  size  of  the  highly 
stressed  jobs  can  be  considered  infinite,  so  that  for  Phis©  I: 


%  =  ©j  Ed,  A(N)3j  A  s  Ap 

where  A  ~  self-propagating  craek  size  at 
p  the  end  of  Phase  2 


(2.17) 


Phase  II,  on  the  other  hand,  can  be  defined  §§  extending 
froa  that  point  in  tins  where  the  craek  has  grown  sufficiently 
large  to  outweigh,  in  its  effect  on  propagation  rate,  the  original 
defect.  Such  a  crack  does  not  require  thB  defect  to  propagate,  it 
is  "self-propagating”.  For  a  c-iek  of  this  size,  or  larger,  the 
size  of  the  highly  stressed  20 .. „  can  no  longer  be  considered 
iafiaitg,  10  that  ono  has: 


0jl  =  0jj(A(N),  S);  A_  &  h  s  Ae 


(2.18) 


where  Ac  s  critical  crack  size 


Here,  kp  is  the  crack  size  at  the  terainstion  point  of  Phase  II. 
This  else  crack  1>  suaed  to  lead  to  a  spall  "instantaneously"  by 
a  precipitous  frastv-.e  seehsaisa.  This  does  not  suggest  that  the 
rolling  contact  systee  bseoaes  inoperative  iaKedintely,  ulthough 
this  egy  be  the  esse.  However,  there  is  a  "estag troph ic"  frewth 
step  betwese  the  crack  of  size  Ag  and  the  cospleisd  spoil,  i.e. 


9  for  A 


(2.19) 


One  c@a  assuee  that  the  critical  crack  size  is  related  to  the  §ise 
of  the  highly  stressed  gone,  or,  in  its  simplest  fora: 


Ae  =  kiS 


(2.20) 


where  fej  =  eeastgat 


. . . . . 


PREDICTION  OF  THE  LIFE  AT  A  DEFECT 


6. 


Substituting  Equations  (2.27)  end  (2.18)  respectl vely,  into 
Equation  (2.13),  oh#  obtgifii  the  following  formulas  for  crack 
propagation  rate  during  Phases  I  and  II: 


'  Y  j  (  ®o  •  a  y 1 


D) 


©J [ d .  A(N)J;  A  *  Ap 


(2.21s) 


en(A(N),  S);  Ap^  A  S  Ac  (2.21b) 


where  subscripts  .1  and  II  apply  to  Phases  I  and  II  respectively 


Using  the  previously  explained  multiplicative  hypothesis  on 
the  function  §,  these  equations  may  be  rewritten  as  follows: 


-  Yj ( gQ  • 


/ 1  (A) 


fCd)  ;  A  £  Ap 


(2.22a) 


Yjj (®e  , 


txyt  D) 


/  8(A/S)  •  /3  ( A ) ;  A  sasa  (2.22b) 

p  c 


In  Equation  (2.22b),  two  functions  fa  and  /3  are  shown,  on© 
representing  the  effect  of  relative  crack  size  by  comparison  to 
the  size  of  the  highly  stressed  zone,  and  the  other  any  remaining 
direct  effect  of  absolute  crack  size  (as  hypothesized  e.g.  by 
Lundberg  and  Palogren  when  introducing  the  effect  of  the  depth 
co-ordinate  z0  of  the  maximum  shear  stress  range.)  Note  that 
in  Equation  (2.22b),  the  original  defect  severity  d  does  not  appear. 
Therefore,  this  equation  contains  only  raaerestraits  and  matrix 
variables,  and  is  independent  of  the  original  defect  population. 

In  Equation  (2.22a),  on  the  other  band,  [Td)  is  different  for  each 
individual  defect,  and  the  equation  is,  therefore,  dependent  on 
the  defect  population. 

Integration  of  the  differential  Equations  (2.22a)  and  (2.22i>) 
leads  to  the  following  forms: 

7 I  <A)  =  N  Y!  r  (d)  (2.23a) 

/jj  <  A,  A/S  )  =  N  Yu  (2.23b) 

Substituting  A  into  Equation  (2.23a)  yields  a  value  Nj ,  the  life 
at  the  end  of  Phase  I.  Substitution  of  Ac  into  the  Equation 
(2.23b)  yields  a  value  Njj,  the  duration  of  Phase  II  life. 


The  life  %  from  ike  beginning  el  cycling  trough  the  ead  e£ 

Phil©  II  life  is  ties  obtained  ai  the  sue  of  these  two  ph a§ 1 i re  a 
=  N|  +  Njx  •  Pine  III  life  say  or  aay  Bofc  be  0,  depending  on 
the  defiaition  of  failure  as  discussed  above. 

Is  principle,  then,  it  is  possible  te  ebtsin  a  prediction  of 
life  te  failure  it  a  particular  defect  with  severity  d,  free 
Equations  (2.23a)  sad  2,23b). 

It  is  noted  that  the  apparent  arbi trgrl ness  in  the  selectlos 
of  the  self-propageti  »g  crick  sise  will  aot  influence  the  total 
life  ~  ^  1 1  the  hypotheses  outliusd  previously  are 

correct,  because  the  Equities!  (2. 23a)  sad  (2,23b)  were  both 
ofeiaieed  from  Ignaties  (2.16)  by  aegleeiing,  for  Equation  (2.23§), 
the  influence  A  ,  usd  for  Equatioe  (2.23b)  tso  isflaeae©  of  d. 
Inasmuch  13  these  spprewiaoiieag  are  volid,  the  two  equatioes 
merely  describe  two  portions  of  the  sas@  function  A(N),  and  their 
domains  of  validity  overlap  §0  that  the  selection  of  Ap  is,  within 
limits,  discretionary. 

Solving  Equations  (2.23a)  and  (2.23b)  for  H,  and  substituting 
§s  described  above,  one  obtains  the  following  expressions  fer  life 
to  failure  ( at  the  esd  of  Phgso  II) 

a  Jl  (Ap)  (2.24a) 

TpTdT 

Nn-  ^XI  Ue»  V*>  (2.24b) 


Nr  =  PL  *  f_ v  (2.24c) 

b  i  II 

One  aast  settle  oa  a  suitable  value  of  Ap  and  determine  that  value 
of  Ag/S  at  which  tbs  eraefe  becomes  critical.  Assuming  fcfe@§© 
deeiiisai  caa  be  reached  geaerally,  Yi  md  f(d)  restis,  ia 
Equation  (2.24a). as  functions  of  ester-sal  parameters  and  S  and  Yg 
is  Equities®  (2.24b).  All  p®ra»@t@rg  of  these  remaining  functions 
are  observable. 
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7,  THE  STATISTICS  OF  LIFE  FOS  AN  ENTIRE  ROLLING  BODY 

In  determining  the  life  of  t  finite-eise  rolling  body,  one 
starts  with  Equations  (2.24)  for  life  in  the  vicinity  of  e  given 
defect,  and  uses  statistical  theory  to  obtain  life  for  e  volume 
of  material  eontaiaisg  a  multitude  of  defects. 

Failure  of  the  rolling  body  will  occur  through  "coapefc i t ioa" 
between  a  sultitude  of  potential  defects  acting  as  failure  nuclei. 
According  to  the  definition  of  Phase  I  fatigue,  gieroeraeks  grow 
at  a  multitude  of  defects,  at  differing  rates,  s§d  independently 
of  each  other.  One  of  those  aicroeraeks,  or  several,  will  reach 
the  beginning  of  Phase  II,  within  the  life  of  the  part.  These 
will  thee  proceed  to  seeuaulate  Phase  12  life  until  such  s  time 
as  one  of  them  has  generated  a  crack  of  critical  size  Ap,  at  which 
a  gpall  forms,  whereupon  the  rolling  body  is  considered  failed. 

All  other  defects  which  have  ©stored  Phase  II  show,  at  the  time 
of  failure,  cracks  of  less  than  critical  size.  In  the  proposed 
model,  variations  ;n  original  matrix  strength  within  a  rolling 
body  are  considered  small  enough  to  be  neglected  (this  position 
can  be  revised  later  if  necessary).  This  leaves  two  sain  sources 
of  variability:  the  systematic  peint=wioe  variability  of  the  macro- 
strain  field  in  the  contact  zone  (and  the  consequent  variability  of 
work  hardening  and  residual  stresses)  end  a  random  variability  of 
defect  severity  and  location  with  reference  to  the  coataet  gone. 

It  is  the  effect  of  these  variables  on  the  life  calculated  from 
Equations  (2.24a)  and  (2.24b)  that  determine  the  outcome  of  the 
competition  among  defects  for  the  generation  of  the  crack  leading 
to  failure. 

A  statistical  treatment  which  can  be  used  to  describe  this 
competition  will  be  illustrated  for  Phase  I. 

Phase  II  life  will  be  considered  a  deterministic  quantity, 
calculable  for  each  poiat  in  the  rolling  body,  from  the  knowledge 
of  macroscopic  gtrsin  sad  matrix  parameters  alone. 

To  express  the  statistics  of  Phase  I  life,  consider  the 
highly  stressed  sons  to  be  composed  of  elementary  "cells”  of 
uniform  size,  selected  small  enough  to  contain  only  one  defect, 
but  large  enough  for  a  crgefe  of  size  Ap  to  be  wholly  confined 
within  the  cell.  Then,  Phase  I  fatigue  damage,  originated  witbiB 
a  cell,  will  remaia  confined  within  it.  Fatigue  damage  existing 
in  one  cell  will  not  influence  the  behavior  of  adjacent  cells. 

On  this  assumption,  Phase  I  life  of  each  cell  is  independent  of 
the  life  of  all  other  cells. 


Assume  sew  that  there  is  a  known  probability  dig tribut ion  of 
defects  of  varying  severities  d  for  e§eb  cell,  i.e.  there  is  a 
beoea  canal  a t i re  distribution  function  Fid)  suefa  that  is  gay  cell 

Prob  (dj  s  d)  =  P(d)  (2.25) 

where  F  =  cumulative  distribution  function  of  d 

Ties,  Equgtioa  (2.24a),  establishing  a  funetloasl  relationship 
between  d  and  Nj,  peraiti  del ersi sit i os  of  a  transformed 
probability  distribution  6(N)  §ueh  that 

Prsb  (  Nj  s  Nj  )  s  €  (  Sj  |  f  )  (2.26) 

where  G  =  cumulative  dig t ribatios  fuaetiea  af  Nj 

S'  =  position  vector 

la  the  genera!  case  shore  the  £uaetisa§  ,  Vi ,  sad  f  is 
Equation  (2.24s)  vary  from  point  to  poist  ia  tie  rolliag  body 
because  of  the  asa^uaiforsi fcy  of  the  aserostrgis  field,  or  for 
any  other  reason,  the  distributioa  S(Nj.)  sill  depend  os  the 
polities  vectors  shown  ia  Equation  (2.26). 

EquatioB  (2.26)  states  the  (cumulative)  probability  that 
the  sell  with  position  eo-ordiasie  s  sill  reach  the  end  of 
Phase  1  life  ia  Nj  eyelet  er  less. 

Frog  Equation  (2.26),  stgteaents  ess  be  sad©  regarding  the 
probability  ef  failure  of  the  satire  rolling  body.  The  rolling 
body  sill  fail  if  exactly  ©a©  of  its  cells  fail.  It  is,  there¬ 
fore,  required  to  express  the  probability  distribution  of  the 
life  of  that  cell  ia  tie  rolling  body  which  fails  first  of  all 
eel  1 g . 

If  tio  fatigue  ph©@0B@Ba  ia  each  cell  are  iadepeadeat  gg 
assuaed,  the  prsbebility  that  the  rollisg  body  i»r wiwee  is  tig 
product  of  tfce  survival  probabilities  of  all  cells  is  it,  i.e. 
ths  life  distribution  of  as  entire  rolling  body  is: 

H  (K2)  =  1  -  a  [  1  -  G,  (Nj)  ]  (2.27) 

where  H  =  cumulative  distribution  function  of  for  the 
entire  rolling  body 
0:  eultiplieat ion  operator 

Equation  (2.27)  follows  froa  Equation  (2,26)  by  observing  that  the 
probability  ef  survival  is  obtained  by  subtracting  the  probability 
@£  failure  frog  uaity. 
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Equation  (2.27)  is  general.  A  specie!  case  is  that  for  shieh 
all  are  identical.  This  is  a  simple  approximation  ox  thg  sags 
ef  a  thrust  loaded  bearing  ring  (the  Herts  stress  field  la 
iadepeBdeat  of  contact  position)  is  which  there  is  a  single  defect 
dis  tribatlos  throughout  the  ring  §ad  the  defects  are  infreqaeBi  so 
that  a  "cell"  can  be  represented  by  a  short  "siiee"  of  riag 
between  two  closely  spaced  crose-sectioBS .  In  this  simple  eise, 
Equation  (2.27)  reduces  to  the  following: 


H  (Nj )  =  1  ~  [  1  -  G  (Nj )  j 


(2-20) 


where  m  =  number  of  cells  in  the  rolling  body 


As  is  shown  in  a  later  section  of  this  report  for  certain 
general  classes  of  distributions  G(Nj),gsd  for  increasing  ra. 
Equation  (2.28)  approaches  the  form  of  a  Weibuli  distribution: 


( .  )  ) 


(2.29) 


H  (Nj)  =  1  -  exp  V  -  - J  i 

where  N*  m-l/k  =  a  "characteristic  life"  or 

scale  parameter 

k  =  constant  dispersion  exponent 
N0  =  minimum  life,  (N0  0  ) 

The  constants  are  determined  by  the  specifics  of  the  distri¬ 
bution  function  G(Nj).  N0 ,  the  minimus  Phase  I  life,  can  be  con¬ 
sidered  sero,  since  a  crack  of  sige  A  may  pre-exist  in  the  matrix. 
Equation  (2.29)  emerges  from  general  theorems  on  the  asymptotic 
properties  ox  extreme  value  distributions  and  is  not  a  separate 
hypothesis. 

In  the  model  presented,  the  (Phase  I)  life  distribution  of  a 
rolling  body  is  not  sorely  observed  ss  a  phenomenological  fact,  it 
is  related  to  the  physically  meaningful  distribution  of  defects. 

The  relationship  between  the  distributions  F(d)  in  Equation  (2.25) 
and  S(Nj)  in  Equation  (2.26)  permits  inferences  from  one  distri¬ 
bution  to  the  other,  thereby  identifying  suitable  measures  of 
defect  severity.  As  shown  ia  subsequent  sections,  it  is  possible 
to  conjecture  appropriate  defect  severity  distribution  and  test 
these  conjectures  by  the  effect  they  have  on  the  life  distribution 
of  the  part. 

Equation  (2.29)  contains  a  volume  effect  on  Phase  I  life. 

For  cells  of  fixed  size,  their  number  m  is  proportional  to  the 
stressed  volume.  Therefore 


N*  ~  V' 


(2.30) 
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i.e.  the  scale  parameter  of  the  Seibull  life  di stribut ion  is  an 
iBfgrgs  function  of  the  stressed  volume  as  in  the  LuBtiberg  = 

Pgiagrea  theory. 

The  observable  Hie  of  @s  eetire  rolling  body  is,  accord¬ 
ing  to  Equation  (2.24c),  the  sub  of  Ph@§©  I  and  Phase  II  lives  §£ 
the  defect  shere  Hill  be  shortest.  This  esa  be  calculated 
§ lapl y.  only  cade?- one  ox  TWe  JoTToaing  tee  conditions  t or  Phase  II 

(  §)  Nj  j  «  Nj  ,  or 

())  Njjjy  '’OBstant 

For  either  condition,  the  shortest  Nl  will  coincide  with  the 
shortest  Nj  ,  so  that  Nj  can  be  taken  from  Equation  (2.27),  (2.28) 
er  (2.29).  Njj  is  then  calculated  from  Equation  (2.24b)  for  the 
jf  at  which  the  shortest  Nj  occurred  and  the  two  are  added. 

If  neither  of  the  above  conditions  for  Njj  is  satisfied,  then 
the  (deterministic)  variation  of  with  sT  soy  cause  the  shortest 
N^.  For  this  case,  the  statistical  treatment  Bust  include  Nj^and 
this  say  lead  to  difficulties  in  satisfying  the  assumption  of 
independence  of  §11  failures  required  for  Equation  (2.27). 

0O  APPLICATION  OF  TBE  LIFE  EQUATIONS 

Equations  (2.24a),  (2.24b),  and  (2.24c),  and  the  statistical 
interpretation  contained  in  Equations  (2.25)  to  (2,29),  represent 
the  framework  of  the  proposed  fatigue  life  prediction  Bedel.  la 
order  to  apply  this  sodel,  it  is  necessary  to  sake  specific  assump¬ 
tions  for  §11  functions.  This  work  will  bo  pursued  in  the  next 
project  year.  St  is  instructive,  bosever,  to  illustrate  the 
approaches  that  eaa  b@  tifeoa  towards  application  of  the  aed@l  by 
eeagiderisg  g  few  special  examples  at  this  time. 

a.  The  Luadherg-PilBgren  Cage 

The  following  is  as  example  of  one  of  several  possible  methods 
by  which  the  Lasdber g-Palagreu  formulas  caa  be  obtained  as  a 
special  ease  of  the  proposed  model. 

If,  in  Equations  (2.24?),  (2.24b),  and  (2.24c),  oa@ 
accuses  that  T(d)  is  a  material  constant,  yj  depends  only  os 
the  tsgxieue  shear  stress  rang©  t0  ,  Ap  is  proportional  to  a 
Hues?  dimension  of  the  highly  stressed  eross  section,  say, 


22 


is  depth  ee^ordi nats  zB  of  the  maximum  shear  stress  range  sad 


f!,t  =  Os  then  one  obtains 


(^°  }  or  Nl  y  (  t0)  /I"1(a0  ) 


coast, 


V  (T0) 


(2.31) 


If  one  uses  power  functions  for  y 


c/e 


Vi 


and  / j ,  then: 


=  const. 


(2.32) 


The  designation  of  the  constant  exponents  is  that  used  in 
Equations  (2.1)  end  (2.2),  Equation  (2.32)  states  that  the  life 
distribution  of  eaeb  defect  aad  consequently  that  of  the  entire 
rolling  body,  is  sealed  by  the  maximum  alternating  shear  stress  To 
and  its  depth  co-ordinate  zQ  exactly  as  specified  by  the 
Lundberg-Palagren  theory.  Equation  (2.30)  states  that  the 
stressed  ?ol ume  is  another  scale  factor  for  life.  Equations 
(2.29)  and  (2.32)  are  equivalent  to  Equation#  (2.1),  (2.2),  ai.d 
(2.4).  Thug,  the  Lundberg-P almgren  formula  appears  as  a  special 
ease  of  the  new  model. 

b.  Deviations  from  the  Lundberg-Palmgren  Type 
Seibull  Distribution 

It  has  been  found  experimentally  that  there  is  a  non-zero 
minimum  life  prior  to  which  there  is  no  finite  probability  of 
bearing  failure.  This  fact  is  not  explicable  by  the  Lundborg- 
Palmgren  theory.  However,  it  is  immediately  obvious  if  Equation 
(2.24b)  is  assumed  to  give  a  non-zero  value  for  NTj.  In  this 
case,  Nt  may  still  be  distributed  according  to  a  ieihuli  distri¬ 
bution  with  No  =  0,  Equation  (2.29),  but  the  total  life  to  failure 
will  have  a  positive  minimum  value. 

c.  Effect  of  Material  Cleanness  (inclusion  Content) 


Assume  that  there  are  two  materials,  one  "cleaaer",  the 
other  of  lesser  elemaess,  i.e.  possessing  different  inclusion 
severity  di § tribui loss  F(d)  par  Equation  (2.25),  but  otherwise 
identical.  The  probability  of  encountering  inclusions  of  grout 
severity  is  higher  for  the  material  of  lesser  cleanness.  One  can 
then  expect  e  larger  number  of  effective  stress  raisers  in  esy 
gives  rolling  body  made  of  the  less  clean  material.  According 
to  Equation  (2.24a),  the  relationship  between  Phase  Z  life  Nj 
and  defect  severity  is  inverse.  Consequently,  the  distribution 
of  lives  G(N)  in  Equation  (2.26)  will  show  shorter  lives  with 
higher  probability  for  the  steel  of  lower  cleanness. 
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The  distribution  of  cell  life  sill  be  soiled  to  low®,  values.  The 
scale  paraaete?  is  Equatiea  (2.29)  will  b©  a  ssalle?  number, 
i.e.  typical  rolling  body  life  sill  be  reduced. 

d.  Two  Cgapetiag  Failure  lodes 

Agsuee  there  is  sot  one  family  of  defects  (e.g.  inclusions) 
but  two  families  (e.g.  laelutioes  aaci  surface  defects).  Assuee 
that  these  two  fasilies  of  defects  operate  1 sdepesdea t ly  e£  each 
ether  gsd  escfe  bag  a  severity  di  stributioj?.  It  is  possible  to 
define  cells  sack  that  they  have  either  one  er  tfea  ether  type  of 
defect,  but  sot  both.  This  state  of  affairs  is  realistic: 
failure  in  rolliay  coattet  has  been  shown  to  occur  either  sub¬ 
surface  (e.g.  from  iaeluiiaas)  or  to  start  at  the  rolling  surface 
(due  te  surface  defects).  Cells  with  ieclusioai  art  volue© 
eleseats  sot  ©steading  te  the  surface,  cells  with  surface  defects 
ire  surface  areas  underlaid  by  §  thin  "ilia"  of  seferisl.  The 
distribution  of  the  two  typos  of  defects  is  ladepeedeat , 

Ohs  obtains  tm  cell  life  distributions  of  the  type  of 
Equation  (2.26):  ose  for  inclusions  aad  the  other  for  surface 
defects: 


Prob  (M,  as  N1<v)  =  Gy  (Njl  «  •  s>  (2.33a) 

Prob  (PJ  s  W | ,g )  -  Gg  (N|  |  s.y)  (2.33b) 

Equation  (2.33a)  applies  to  subsurface  boluses  (three 
e©°©rdinat@g)s  iitd  Equation  (2.33b)  applies  to  surface  areas 
(two  ee-ordi aatee ) . 

The  probability  of  relliag  body  failure  can  be  obtained 
by  considering  that  survival  of  the  rolling  body  necessitates 
survival  ©f  all  calls  from  both  pspulatioas.  Aeeordiagly,  is 
iialogy  te  Equation  (2,27) 

H  ,Nj)  a  l-  (  II  Cl  -  Gv  (Nj)  J  n  [1  -  Gs  ( N, )  ])  (2.34) 

V  S 

or  substituting  Ires  Equities  (2.20) 

H  (Kj)  s  I  -  (  [I  -  Gv  (Nj)  n""  [1  ~  G#  (Nj)  3"*  )  (2.35) 

share  av  =  ausbtr  of  cells  with  inclusions 

m$  s  number  of  cells  with  surface  defects 
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One  can  determine  whether  the  resulting  asymptotic  distri¬ 
bution  for  §  large  number  of  cells  is  a  Weibull  dig tributi on  by 
comparing  the  h arae t eri s li e §  ox  the  cell  life  distributions  6^ 
and  Og.  will  be  Seibull  only  if  and  Gg  are  suf f i ci ea t i y 

similar. 

e.  The  Effect,  of  Hesidual  Compressive  Stresses  in 

the  Material 

If  there  ere  residual  stresses  in  the  highly  stressed 
zone,  their  hydrostatic  pressure  component  can  he  expected  to 
influence  ductility.  Accordingly,  yr  and/or  Yjr  will  change. 

This  may  have  as  effect  on  crack  initiation  (Phase  I),  or  only 
on  crack  propagation  (Phase  II).  Depending  on  experimental 
evidence,  appropriate  modifications  can  be  introduced  in  the  first 
or  the  second  of  Equations  (2.24). 

It  is  also  possible  to  account  for  the  effect  of  the  non- 
hydrostatic  component  of  residual  stresses  if  a  reasonable 
assumption  can  be  made  regarding  the  effect  of  a  static  stress 
component  on  the  mi cropl as t i c  strain  §0  .  As  said  earlier,  many 

fatigue  theories  assume  that  this  effect  is  nil. 

f.  Effect  of  Hardness 

Three  parameters:  e0,  Gy  and  D  enter  the  matrix  strength 
function  V  .Equation  (2.11).  Clearly,  the  assumption  made  by 
Lundberg  and  Palagren  that  the  only  s t ress»v ari abl e§  influencing 
life  are  the  maximum  alternating  elastic  shear  stress  and  its  depth 
eo-ordinste,  is  tenable  only  if  the  material  is  kept  a  constant 
so  that  ciy  and  D  do  not  vary.  For  different  materials,  y 
will  depeed  on  the  excess  of  elastic  stress  above  the  mieroyield 
stress  Gy  ,  which  is  known  to  be  related,  although  not  equivalent, 
to  indentation  hardness.  Thus,  a  general  expression  for  y  will 
contain  a  material  strength  parameter  of  the  type  of  hardness. 

g.  NsG=Hertai an  Contacts 

Practical  contacts  are  nen-Hertzian  in  two  respects: 

Roller  to  race  contacts  are  macroseepi  cal  ly  non=-Hertzian 
because  roller  profiles  are  ao t  correctly  approximated  by  second 
order  surfaces  (roller  crowaiag,  edge  effects.) 

All  contacts  are  ei eroi eep i cal  1 y  non-Hertsian  because  of 
the  presence  of  surface  asperities,  whieh  will  be  discussed  latc-r. 
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No  Hertsiifi  asiaaptioa  la  inharnnt  is  any  e£  the  deriv§~ 
fines  gives  abe?s.  Provided  that  the  total  straits  eaa  be  defined 
peiat-wii©  ia  lie  e@at§et,  it  if  passible  to  calculate  eeli  life 
aad  lease  rolliag  bed?  life,  ugisg  the  prepared  serial. 

h .  Lubricities  Effects 

The  most  tberseghlj  explored  lubricities  effect©  t@  date, 
are  those  ef  el  as tssydredf naet e  1 ubri cat i e b ,  i.e.  tbs  trsasaitta* 
of  eoa-tact  pressure  via  §  pressurised  viscous  file  of  lubricant. 
Tie  effects  ef  such  §  file  ea  lif©  e§n  be  treated  oa  two  levels: 
tie  saereseepic  red! itributiea  of  pressure  brought  about  by  the 
eia§  tohydrodynasi  c  film  caa  b®  used  to  correct  the  saeresoopic 
fsireia  field  ia  the  sstris  said  defiees  the  fusctloB  y . 
iiersseepieilly,  it  is  possible  to  describe  the  influence  of  a 
sarli ai  ©1  as  tehydrodyBasii  c  fils  ©a  asperity  i steree ti obs .  The 
result  i§  a  sodificitios  ef  the  eierostress  field  ia  tie  viciaity 
ef  tk©  §«rfae©,  aad  s  ses  distribution  of  cell  life  is  the  pepula- 
tiea  ef  cells  eoatsislsg  surfge©  defects.  In  priasiple,  it  is 
ting  possible  e@  cgleulat®  these  el astohydrodysaaic  effects  on 
relligg  body  life, 

i.  Sig©  Effects 

Because  the  present  study  is  ultiaately  §iaed  at  improve- 
seats  ia  radar  auteagg  heerisgs,  aad  radar  an  tew  as  bearings  are 
aaeng  the  largest  aid©,  gig©  effects  are  ef  greet  isperttaee, 

Th©  aedel  presented  offers  several  eluss  to  giae  affects. 

The  self-propagetisg  ersefe  gig©  A  is  defined  eg  the 
smallest  crash  which  progresses  at  a  rgte  Isdepegdeut  af  the 
original  defect.  Clearly,  the  sagsitude  of  this  eraefc  aust  depend 
ss  tfe©  ^effective  radius”  ef  the  most  severe  ©rlgieal  defect,  i.e. 
it  east  be  related  to  tie  volume  wit&is  which  the  stress  raigiag 
effect  ef  the  defect  is  felt.  This  velum©  almost  eertaialy  is 
propsrtieaal  t@  defset  ditm©t©r.  Aecordiagly,  depends  ©a  the 
siae  distribution  @£  the  erigisBl  defects.  Since  eaters  isto 
II©  deteralasliea  of  16©  Phase  I  lif©  ef  e  defect,  It  predue©s  a 
siae  effect  oa  this  life,  Absolut©  beariag  si§©  will  eater  inte 
this  effect  iaaseueh  is  it  isfluewees  steal  processing  practices, 
sad  thereby  the  gig@  distribution  of  defects. 

Tfe ©  other  Phase  I  gig©  effect  coaceras  the  cell  siae 
istredueed  la  eesaectioa  with  Eqattiaa  (2.25),  Since  cells  should 
eeatgie  only  ea©  defeet,  their  §ig@  s§  related  t©  the  defect 
spicing.  This  gpieiag  is,  agile,  iaflueaeed  by  absolute  fceariag 


gig©  through  the  correlating  between  reductiea  Is  srsss  iectioasl 
are®  uadsrgasa  by  the  steel  durlay  rolliag  ©r  fergiag  ana  defect 
dig Eributioa.  Cell  iige  indirectly  eaters  Equation  (2.38) 
beciuse  tie  aueber  of  ceils  present  is  a  given  relliag  body  depends 
oh  their  gige.  Of  eeune,  the  distribution  F  (d)  iu  Equation  (2.25) 
also  dapesdg  ea  cell  §iie  so  that  the  effect  is  set  a  simple  oh©. 

Gives  a  defect  sige  distribution  and  a  ceil  sise,  the 
stressed  volume  of  the  rolliag  body,  in  terpss  of  aultiples  of 
unit  cell  volume,  introduces  the  volume  effect  of  life,  repre¬ 
sented  by  Equation  (2.30). 

Turning  now  to  sis--  effect  on  Phase  II  life,  Equstioa 
(2.24b)  shows  two  siss  effects.  Oae  is  represeated  by  Ae 
It  is  tie  (hypothetical )  effect  of  absolute  crack  size  on 
propagation  rate.  The  ether,  aore  obvious  effect,  is  represeatsd 
by  the  ratio  As/8  tad  represeats  the  fact  tast  8  eraefe  east  grow 
to  §  certain  sigs  with  reference  to  the  size  of  the  highly 
stressed  soao  before  a  spall  can  for®  (the  crack  sust  propagate 
from  the  depth  of  maxisus  shear  "tress  to  the  surface,  er  vice 
versa,  since  spalls  are  typicall  of  a  depth  comparable  to  that 
of  the  location  of  eaxisue  shear  str-sses).  This  effeet  is 
related  to  absolute  bearing  sise  through  5. 

There  are,  af  course,  auasrous  size  effects  lap! led  in 
Equities  (2.15)  relating  total  sacrostrain  to  external  load  and 
contact  geometry,  lost  of  these  effects  are  accounted  for  la 
the  Luodberg-Palsgrea  theory.  Steel  processing  effects  are  also 
implied  in  Equities  (2.15)  through  the  feet  that  the  yield 
strength  ©y  say  sell  depend  oa  the  absolute  size  of  the  part. 

The  same  applies  to  tbo  duelility  parameter  D  in  Equation  (2.14), 
Clearly,  the  proposed  model  offers  ample  rooe  for  the  exploration 
of  gigs  effects,  providing,  ef  course,  that  the  necessary  expert- 
seats!  evidoace  eats  ba  secured. 

9.  THE  STRESS-ST8AIN  RELATIONSHIP 

The  relationships  between  esteraal  load,  resulting  elastic 
stregaes,  total  strain  s0  defined  in  Equation  (2.15)  and  aicr©- 
plastic  strain  sc  defined  in  Equation  (2.9)  are  key  elements  of 
the  model.  In  the  present  state  of  the  work,  a©  choice  has  been 
mode  among  the  possible  approaches  by  which  ®n  eaa  be  related  to 
a  ealeulateu  stress  field.  However,  it  seems  obvious  that  high 
plastic  strains  sill  be  associated  with  high  shear  stresses,  and 
considerable  effort  was  spent  is  identifying  sources  of  high 
sheer  stress  uader  acre  general  conditions  this  previously  avail¬ 
able,  Three  stress  analyses,  described  in  subsequent  sections, 
are  of  Importance. 
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g.  Coat ours  of  Equal  Shear  Stress  Bunge  is  Elliptical 

Harts  Contact 

The  proposed  sedal  requires  g  definition  of  5,  the  highly 
stressed  gene.  A  logical  definition  may  be  the  gone  bounded  ia 
each  crass  section  by  §  suitably  selected  line  of  equal  shear 
stress  range,  and  extending  around  the  rolling  body.  Since  calcu¬ 
lations  for  such  linos  are  aot  available  free  the  literature, 
they  hose  been  computed  and  gre  shown  in  this  report.  Numerical 
equations  and  graphical  relationships  have  beea  developed  fra® 
which  areas  &,  bounded  by  selected  centaurs  of  equal  shear 
stress  ru..ge,  can  be  obtained.  One  say  in  which  these  contours 
day  bp  used  is  by  observing  that  the  distorted  tubular  arnulus 
between  any  two  contours  is  a  zone  of  roughly  equal  shear  stress 
range,  sad  therefore  presumably  equal  total  macrostrain  *o • 

It  is,  therefore,  as  area  ia  which  the  functions  Yy  and  yjj 
of  Equations  (2.24a)  aad  (2.24b)  esn  be  considered  emstant  (with, 
perhaps,  s  secondary  influence  of  the  variable  hydrostatic  com¬ 
pression  on  ductility.)  A  correct  way  of  evaluating  Equation 

(2.27)  to  obtain  rolling  body  life  might  be  to  consider  the  dis¬ 
tribution  functions  G  |  ( N j) i doa % i e a  1  for  cells  within  each  of 
these  aaauli  of  equal  shear  stress  range; to  apply  Equation 

(2.28)  to  one  annular  set  of  eells  at  a  time,  and  then  to  coapute 
H(N^'  from  Equtuioa  (2.27)  by  multiplyiag  the  survive!  probability 
functions  for  each  annulus. 

b.  Near-Surface  Stresses  Due  to  Asperities 

A  major  result  of  research  into  rolling  contact  fatigue, 
referenced  later  is  this  report,  Has  been  the  recognition  that 
destructive  spalling  fatigue  failure  of  rolling  bearings  is 
often  preceded  aad  precipitated  by  "surface  fatigue",  $  sequence 
sf  starting  with  plastic  flow  in  the  Jasadiate  subsurface 

layers  of  material t  followed  by  profuse  ei erocraoki ag  at  the 
surface,  and  eventually  leading  te  spell  feraatios  for  surface 
engijsted  cracks.  This  series  of  phenomena  was  found  to  be 
) ..  brl  cut  io  o  related  and  has  been  associated  with  the  thickness  of 
tse  el  as  t  >- ;\  •*}  rodyn  esi  e  lubricant  fils,  compared  to  tee  roughness 
si  the  cse  tl  ag  surfaces.  Until  recently,  it  was  thought  that 
el  as  tshydr?r:Tn§rai  e  films  prevent  surface  fatigue  by  reducing 
tractive  (tangential)  forces  at  the  surfaces  which,  in  the  gb,.©ne© 
of  n  full  .1  agtF’Byiirodynami  c  file,  are  transmitted  between  con- 
tectliu;  ■.<•) -c  ,*i  1 1  «•«  •  Because  the  magnitude  of  these  hypothetical 
tract)  liives  cannot  be  determined  generally,  aad  because 
everei  pfoeos  w  /-cot  rarii  c tory  evideace  eould  net  be  gceofasedaiod, 
ssetfeesat i  •  ►  '  KO«’.@ling  of  aurfgee  originated  fatigue  was  at  g  standstill. 
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1§  tfee  last  sasils,  it  Mag  beds  proposed  that  lurfecs  fatigue 
©eeur-r  encse  fee  explained  by  purely  sorgsl  iertgita  preacur®, 
acting  is  agperity  dimessisHg.  Tfa®  role  of  tie  el  tg  tohydredysafii  e 
film  i  a »  Sa  ill®  vi@w,  that  it  preveatt  ghs?p  pressure  gradleatg 
fros  gr^ciag  it  each  §ap§r§i@  asperity,  wigr@n  th©§©  pressure 
g  radl  @a ..  s  do  occur  whenever  twe  asperities  cose  laio  contact 
through  a  boundary  lubrieaat  film.  This  hypothesis,  if  proven 
correct,  can  serve  ai  a  basis  sf  taatheiaatieal  modeling  for  sur¬ 
faces,  To  explore  its  workability,  stress  analysis  was  conducted 
on  a  raodel  asperity,  represented  by  a  prisssstie  ridge  topped  by  a 
small  radius.  It  is  shown  in  subsequent  sections  of  this  report 
that  high  near-surf  gee  shear  stresses  arise  under  the  aerial 
asperity  whan  it  is  pressed  against  a  flat  surface  (or  opposing 
asperity),  and  that  the  magnitude  of  this  shear  stress  depends  on 
two  parameters:  the  degree  of  depression  of  the  asperity  (in 
partial  elastshydredyaaBsic  coetoct  this  is  determined  by  the  film 
thickness  to  roughness  ratio)  and  by  the  typical  slope  of  the 
asperity  side.  Experimentally,  it  was  shown  that  asperities  on 
relatively  rough,  e.g.  li-ground,  surfaces  have  steep  slopes  of 
the  order  of  30°,  asperities  on  finely  honed  surfaces  have  only 
slopes  of  the  order  of  4°,  whereas  the  finest  achievable  lapped 
surfaces  of  hearing  bails  have  asperity  slopes  of  less  than  1°. 
Calculation  shows  that,  the  shear  stresses  under  ground  asperities 
reach  the  yield  strength  of  hard  steel  for  relatively  little 
asperity  depression,  those  under  asperities  of  honed  surfaces 
break  into  the  plastic  range  at  substantial  depression  values 
(low  film  thicknesses)  only,  whereas  those  on  lapped  surfaces 
should  not  besose  plastic  under  any  conditions. 

In  principle,  it  should  bo  possible,  based  on  this 
analysis,  to  calculate  the  life  of  calls  Immediately  adjacent  to 
a  rough  surface,  based  an  Equations  (2.24a)  to  (2,24*). 

The  sequence  sf  ©vents  leading  to  surface  originated 
failure  is  Isis;  Near-surfae©  plasticity  occurs  under  contacting 
asperities,  and  eventually  causes  slcrocrackg  at  the  surface  to  g 
shallow  depth  of  several  hundred  si  ere! aches .  The  life  prior  to 
the  formation  sf  these  creeks  can  b@  calculated  free  Eqastlsas 
(2.24a)  ts  (2.24b).  The  fersatios  of  surface  fatigue  cracks  (not 
to  b©  confused  with  deep  cracks  which  say  originate  at  the  surface 
sad  will  be  discussed  a nbsequeatly)  represents  a  failure  ph§g@ 
prior  to  Phase  I.  Experi aestall y,  it  is  a  distinct  phage  both 
s@r©  widespread  and  hot©  rapid  of  progression  than  the  subsequent 
spelling.  The  surface  ©quipped  with  Biereerseks  can  be  treated 
as  a  source  ©£  Pfe§§®  I  crack  generation  because  the  aieroeraefes 
are  so  genii  sad  shallow  that  they  do  sat  shoe  §s  apparent  tendency 


I©  xtts£§@se  i|§  sarf of  settriii  ting  requires 

@  Isrg©  population  ©£  gts  defects,  vis.  tie  gurfteg  Igtigyg  ( rsei 
Om©  eai  tlas  gaits?  Igaatjeai  (2„24a)  to  (2.24a)  again  §ad  ealeu- 
Igl©  ipil 1 i gf  f it I y a©  life  £©?  1 1  i  §  population  af  defects . 

8>  Stresses  Dade?  a  Surfgc©  leperfeetlea 

Aside  £fo@  the  geugraliaad  rsugasagg  of  the  §  urf  ece, 
tierg  ©r©  loeillggd  lapit I aeti g s§ .  generally  of  ifes  scratch  ©? 
fis-rsfeii  type,  §a  §f® ry  practical  bearing  surface.  It  has  bean 
recognised  recently  that  sucM  furrows  or©  points  of  origin  £tr 
spells.  They  ?@preg©st  §  iteoad  population  e£  defects  (alluded 
£©  previously)  which  eisptteg  with  the  subsurface  defects  of 
the  iaslisgits  type  to  geserata  spalling  failures.  Saerees  the 
Laa4b©j’f»p §Ia§re a  theory  is  eiptbl©  sf  predicting  life  icr 
failures  erigiBitiig  at  iaeluiloss  by  using  the  macroscopic  shear 
stre&i  rang©  la  lieu  if  the  total  strain  parameter  es  ,  end  by 
assusisg  that  tbs  str©§§  raising  properties  ef  the  inclusions 
(their  severity)  §?©  a  material  §oa§ig@t(  th©  approach  is 

set  feasible  for  surface  originated  failures  because  there  is  bo 
kifb  shear  ot?@§8  at  the  surfae©  of  s  Hertzian  contact  between 
ideally  isootfe  surfaces  ia  th©  abitae©  ef  traetiaa.  This  diffi¬ 
culty  diggppggrg,  heweytr,  If  it  can  be  shows  that  high  shear 
sifts/iee  arise  in  the  vieialty  of  localised  surface  i ape? factions 
^h§i  they  eattr  the  Htrtaiga  contact  area.  The  stress  analysis 
gifts  is  §  later  section  of  this  r©p©rt  has  shows  that  this  is 
tfte  ©@s@.  rhe  iarfasi  stress  field  was  ealeul gt@d  under  a  furrow 
type  inpogfoctioa,  represested  by  §  lesg  prissatle  dtpressioa  in 
a  pliM.  ‘feo  depression  has  two  rouaded  edges  gad  sufficient 
depth  t@  prevent  contact  at  the  bottea  btlow  the  reuaded  edggg. 
High  sheer  strggggi  were  found  usd©?  the  rousded  edges.  §atl  their 
fiiatioigfaip  to  the  geoagtrle  pi?ae©terf  of  tie  defect  were 

esleulati©!  permits  the  assignment  of  §  severity 


itrsis  0 tiers! i all oi 

^eeeriiiig  to  Equalise  (2.8),  the  plsstie  strain  physically 
relevant  to  eraefe  propiifiEies  r§E©  ig  e08  the  atreia  at  the  site 
ef  ih@  erieh  frest.  it  ig  eesvesieot  gad  glspi©  to  visualise  this 
sit  of  s  pre-esistisg  total  strata  $s  |g  tte 
h  e©  whieh  epefitgg  the  strain  rgisiig  ©ffaet  a# 


Tii©  actual  strain  field  la  the  vicinity  ef  an  i  super  fee- 
ties  ii  quite  complex  and,  even  if  the  macros irai u  in  the  vicinity 
of  tfie  inclusion  egg  be  considered  ueifnrs,  depends  oa  all  three 
pri-.i^sl  styaisi.  Thus,  geaernlly 

§c  -  /  («i «o  <?y«  ©i e0  >  (2 . 36) 

with  different  defect  severity  factors!  applicable  along  each 

principal  slrairs  axig. 

The  preceding  paragraphs  contain  two  examples  for  th© 
calculation  of  a  do? eet-induced  maximum  shear  stress,  bgsed  on 
defect  geometry  (i.e.  severity)  and  assuming  n  simple  macrostress 
field.  Given  this  maxi  bub  shear  stress  at  the  defect,  it  may  be 
possible  to  put  forth  a  simple  model  of  plastie  flow  at  the 
defect,  yielding  ec. 

Bbetlier  it  will  be  necessary  for  the  application  of  the 
model  to  proceed  to  the  actual  determination  of  the  plastic  strain 
magnitude  at  defects,  or  whether  calculation  of  a  maximum  shear 
stress,  attached  to  known  defects,  is  a  sufficient  refinement  of 
the  Lundberg^P almgren  aaeroshegr  stress  criterion  to  accommodate 
currently  available  experimental  evidence  is  a  subject  for 
future  investigation. 

10.  OUTLOOK 

k  mathematical  model  for  the  prediction  of  the  life  of 
rolling  contacts  h©s  been  proposed,  based  on  a  concept  of  crack 
propagation  from  pre-existing  defects.  Numerous  parameters 
characterising  the  material  and  geometry  of  the  contact  and  the 
operating  conditions  are  incorporated  in  the  model,  including 
lead,  lubrication,  matrix  strength,  defect  population,  and 
defect  severity. 

Further  work  now  underway  is  aimed  at  closer  definition  of 
as  many  of  the  influential  variables  us  feasible,  considering 
currently  available  information.  Starting  with  the  Lundberg- 
Palmgrea  predicting  of  bearing  life,  the  incorporation  of 
experimentally  documented  effects  not  previously  accounted  for, 
will  be  attempted,  @@d  prediction  will  be  correlated  with  life 
test  results. 


SECTION  III 


SYNOPSIS  OF  LUNDBERG-P ALMOST  fHEOHY 


The  first  iacee&sfui  syitiaatie  attempt  to  treat  rolling 
bearing  fatigue  life  analytically  was  made  by  G.  LHndber@  and 
A.  Palwgrea  is  1947,  under  sponsorship  of  AB  G §  lP,  Gothenburg, 

Swedes  (1)*.  This  work  was  further  pursued  with  e  special  vies 
towards  roller  bearing  fatigue  life  end  reported  In  (2). 

This  wo*-R  of  Luodbsrg  and  Polsigren  is  the  basis  for  the 
life  prediction  method  standardized  by  the  Anti-Friction  Bearing 
Manufacturers  Association  (AFBMA)  (2),  the  AS \  < 4)  and  ISO  (5) 
for  computing  rolling  bearing  life. 

1.  FAIkUBE  PROBABILITY  DISTRIBUTION 

The  fallowing  is  Lundfasrg -Paling  ran  !s  developm/sit  of  the 
probability  §(L)  that  a  bearing  ring  will  survive  to  life  L. 

S0  is  today's  parlance  is  terasd  a  reliability  function.  Its 
arithmetic  complement  F(L)  =  1  -  S(L)  is  the  cumulative  form 
of  a  failure  probability  function. 

Let  X (N )  be  a  hypothetical  function  which  describes  the 
fatigue  ”008411108”  of  a  differential  volume  &V  of  a  ring  or 
rolling  element  material  at  a  depth  z  below  the  rolling  surface 
after  N  cycles  of  stress  are  endured  by  that  volume.  Let  AX(M) 
represent  the  change  in  the  material  fatigue  condition  within  u 
small  number  of  additional  cycles  AN. 

The  probability  that  the  volume  develops  a  "failure”  (i.e.,  a 
crack)  in  the  interval  (N,  N+AN  )  is  tuken  to  be 


/(  \<N)  ).JU<N)-aV 


The  probability  of  not  developing  a  failure  (creek)  in  this 
interval  is  the  arithmetic  complement,  of  Equation  (3.1),  i.o. 


1  -  /(  x(N)  ).A\<N>'ov 


^Numbers  in  paranthese s  refer  to  the  References  at  the  end  of 
this  report 
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Changes  in  the  material  condition  at  a  depth z  are  token  to 
depend  on; 

a)  maerostress  t(N)  which  is  most  dangerous  from  the  point 
of  view  of  material  fatigue.  This  stress  is,  hosed  on 
observations ,  taken  to  be  the  alternating  shear  stress 
in  the  direction  parallel  to  rolling, 

fc)  tin  material  condition  T(N), 

c)  the  depth  ^ 

Thus  the  change  in  material  condition  di.(N)  for  a  small 
number  of  loed  cycles  UN  may  be  written: 

=  j  (  KN),  T(N),  2  )  (3.11) 


The  reliability  function  SCL)  is  uniquely  determined  from 
Equations  (3.10)  and  (3.11)  if  the  functions  G  and  J  are  known. 

It  is  postulated  that  these  functions  are  power  functions 
«s  fol lows  : 


•g(x(N>)  ~  k(\(N)) 

-) 

j(\(N),  T(N).z)  =  (\<N))  -k(t(N)#3) 
ghere  g  uni  j  are  constants. 


(3.  12) 

(3. 13) 


i 

i 


5 

i 

\ 

y 

\ 


35 


From  Equations  (3.11)  and  (3.13) 

i 

(UN))  d X ( N )  -  N  ) e  2 ^dN 


(3. 14) 


From  Equations  (3,10)  and  (3.12),  the  condition  S(N  =  0)  =  1 
requires  that  X(N  =  0)  =  0 

Integration  of  Equation  (3.14)  between  the  limits  of  0  and  IS 
therefore  givess 


-  f  K(T(N),z)dN 

J+1  0 


(3. 15) 


introducing  j  +  l  -  gives 


./  K( T<N)(*)dN 
0 


8 


(3,16) 


\<N)  - 

Using  Equations  (3.12)  and  (3.16)  in  Equation  (3.10)  gives.* 

log  -±—  =  k(«  ].j  |  (  r  tt(  Tin-;,  a  ;  an  j  f  av  (3.17) 


S  k(a)*f  (  fi(  T(W),  a)  dN) 
S(N)  '®  JU  '’o  ;  J 


If  the  amplitude  ^(N)  is  independent  of  N,  equation  (3.17) 
becomes j 


log  ^  dV 


(3.10) 


la  u  Hertz  stress  field  of  given  geometry „  th©  value  of  t(g,>c) 
at  depth  s  and  coordinato  position  K  measured  perpendicular 
to  the  rolling  direction  from  the  midpoint  of  the  contact 
ellipse  is  given  by; 


T 


?g  ■/(£*  } 


8  « 


(3, 19) 
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where  "0  and  z0  represent  the  maximum  shear  stress  arapli' 
tuda  sad  a  Is  th©  half  axil  of  the  contact  ellipse  is  the 
direction  aeroia  the  raceway. 

Using  Equation  (3. 19)  in  Equation  (3.18)  gives; 

“«  OTT  ~  K"('T»* ’<?,  f.) 


where  10  is  the  length  of  the  raceway 
Introducing  the  change  of  variables  P  =  ■&•  V  = 


gives  : 


where 


J  Ke<  T.  2)  dV 
v 


Z0-l0-«p(  T0, 20  } 


«P(to«*o'  =  K°  T0f  (u.v),  Z0V  du  dv 


(3.20) 


(3.21) 


ing  Equation  (3.20)  in  (3.18)  gives  • 


log-™  N ®a  •Zo'lo'P^  To*zo) 

S(N) 

The  relationship  log  or; 

£ 

F(N )  =  1-.»S(N )  *  l-exp-fflj) 


defines  the  probability  distributios  of  bearing  lives  when  many 
identical  bearings  are  operated  under  identical  conditions.  This 
distribution  is  today  known  a§  a  (two  parameter)  Weibuil  distri¬ 
bution. 


It  is  postulated  that; 


fP  (T0, 2s)  “  T, 


c  ~h 
0  ’  so 


(3.22) 


whore  c  and  h  are  unknown  positive  constants  which  satisfy 
the  intuitively  reasonable  relationship  that: 


S(N)  -  1  if 


0  or  z«  =  @  • 


3T 


. . 


Us  lag  Equation  (3,22)  ia  iquatioa  (3  0  2 1 )  fiftss 


(3.23) 


At  this  point  ia  the  developsoRti  the  quantities  fo  ant*  g9 
are  related  through  Herts's  equations  for  the  contact  of  elastic 
bodies,  to  the  bearing  geometry  (rolling  body  diameter  08,  raceway 
dieeeter  On  and  raceway  curvature  ia  a  pleas  perpsadicul ar  tc  the 
rolling  direction)  she  elastic  constants  of  the  material  (Young's 
modulus  and  Poisson's  ratio)  and  th©  contact  force  §.  Th®  contact 
force  ia  here  assumed  constant  end  indepeadsat  as  position  oh  th© 
ring.  The  notion  of  as  equivalent  load  is  later  introduced  to 
acoouat  for  those  coses  (@,g,  the  stationary  outer  ring  of  a 
radially  loaded  bearing)  where  load  varies  with  sing  position. 


Also  introduced  ia  th©  number  of  contact  cycles  per  revolu 
tins  u  dsfiaed  through  the  rei at ionshi pt 

N  =  uL  (3 

whore  L  is  bearing  life  in  Billions  of  revolutions. 


(3.21) 


Fcr  point  contacts  the  hslf^width  of  the  contact  ellipse  is 
replaced  by  Its  expression  ia  tores  of  load  from  the  Hortaian 
equations. 

la  th©  1947  tsoutmosto  the  half=width  s  for  line  contact 
is  lakes  equal  to  three  quarters  of  £h©  roller  length,  the 
19 52  ©stesuioB  of  the  wash  deals  specifically  with  roller  beariags 
and  Introduoes  a  now  expression  relating  stress  and  ©ffeetive 
roller  leagth. 

For  point  contacts  chef©  gubst itaiioos  lead  (for  th©  eosttet 
at  either  bsariag  ting)  tos 


log 


ss+k-S  c  ~h-f i 

#.0a  0=6+*“.  Q  § 


(3.25) 


whore  I  eoataiai  bearing  geoeeferiesl  parameters  to  powers 
which  ar©  linear  oosbigatioBe  of  th©  exposenis  e8  h  and  e.  Th© 
rolling  body  diaetter  app@ars9ls  th©  function  only  Is  the 
fors  Ua/ds,  where  da  is  th©  pitch  diastler  of  the  rolling  ©leaents. 
The  survival  probability  gives  by  Equation  (3,25)  is  that  of  the 
risg=rolIiag  eleaeaf  eestset.  Both  th©  ring  esd  the  rolling  element 
are  goasldesed  equally  likely  to  fail. 
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Qe  li  de  f i as  a  to  b©  the  ring  load  Q  for  which  the  life 
for  tile  fined  vslue  S  -  0,9  is  oee  million  revolutions. 

The  life  L^q  for  the  s age  survival  probability  uatigr  a 
load  Q  Ls  found  fro@; 


e~h-H 


s  ”  h +i 


0  9  “  Lio 


L  _(£s.)  ■’  =(V)P 

LIO  ~'Q  '  '0 


(3.26) 


where  p  -  £±jii£. 

8  © 


The  rolling  element  load  Q  is  proportional  to  the  bearing 
load  F,  and  hence,  Ci  or  Co  the  dynamic  capacity  of  the  (inner  or 
outer  ring)  contact  .defined  as  the  bearing  load  for  which  the 
ring  will  endure  one  million  revolutions  with  a  survival  probabi¬ 
lity  of  0,90a  is  proportional  to  Qc.  Accordingly,  from  Equation 
(3.26); 


(3.27) 


2.  EQUIVALENT  LOAD 

From  Equation  (3.25)  the  logarithm  of  the  reciprocal  of  the 
survival  probability  is  proper! ionQl  to  whore  6J  = 

end  Q  may  vary  with  angular  position  %  on  the  bearing  ring,  i.@. 
Q  =  Q(f>. 


In  view  of  Equation  (3.23),  the  summing  over  the  complete  ring, 
of  probabilities  that  individual  ring  segments  of  angular  length 
will  survive,  gives  rise  to; 


leg  4  ~  J*  Q*<t>I?  JO.  d*  =  Qe  ttD„L 


(3.20) 


where i n 


(3.29) 
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Oe  is  thus  the  equivalent  constant  ring  load  for  which  the 
ring  contact  survival  probability  Is  the  same  os  it  is  under  the 
spatially  variable  loading. 


a 

S 


I 


A  continuous  load  distribution,  the  components  of  which 
integrate  to  equal  the  applied  redial  and  axial  load,  is  intro¬ 
duced  as  an  appvoximot ion  of  the  discontinuous  rolling  element 
loads  and  the  integral  of  Equation  (3.29)  is  evaluated  as  a 
function  of  the  ratio  of  radial  and  axial  load. 

3.  CAPACITY  OF  A  COMPLETE  BEARING 

The  probabilities  Si  and  30  that  the  bearing  inner  and  outer 
ring  contacts,  respectively,  will  endure  beyond  a  life  L  under  a 
benring  load  F  are  given  by  • 


i 

t 

j 

1 


log  J.  =  lrf  j?wLe 


(3.30) 


log  _i  =  k0F*Lc 
0 


(3. 31 ) 


where  k j  and  kQ  are  constants  of  proportionality. 

The  probability  S  that  the  complete  bearing  survives  to 
life  L  is  the  product  of  S0  and  Sj  hence  : 


log 


1  =  (  k1+  ke>-  Fw.  L 


(3.32) 


By  definition  whan  in  Equations  (3.30)  to  (3.32)  the 
survival  probability  is  taken  *.q»oi  to  0.9  and  the  life  L  1 
the  loads  will  be  equal  to  the  respective  dynamic  capacities. 
Thus, 

i»«oS  =  k*'c"  =  ktfco  =  (  k«  +k°  )>cW 


from  which  it  is  found  that  ; 


(3.33) 
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4  ^ 


DET‘-;BB  I  NATION  Of  COf'STANTS  IN  [HE  LIKE  FOHHUl.A 


a.  Point  Contact 

Taking  logarithms  of  Equation  (3.25)  gives  log  log  g-  elog  L 
plus  terms  not  containing  L»  e  is  the  parameter  characterizing 
the  dispersion  of  lives  L  of  a  sample  of  identical  hearings  oper¬ 
ated  under  identical  conditions.  e  is  today  generally  designated 
as  the  "We! bul l  slope”  of  the  life  distribution.  Bearing  life 
test  results  yield  an  estimate  of  e  as  the  slope  of  the  straight 
line  which  is  obtained  when  the  percentage  of  unfailed  bearings  is 
plotted  against  life  on  paper  so  ruled  that  the  ordinate  is  propor¬ 
tional  to  log  log  |  and  the  abscissa  to  log  L  (such  a  diagram  is 
called  a  hoi  bull  plot). 

From  Equation  (3.27)  it  is  seen  that  the  exponent  p  -  - '  may 
be  determined  from  tests  conducted  under  various  loads  F  as  the 
slope  of  the  line  obtained  when  L  and  F  are  plotted  on  log-log 
paper. 

3 c+ h— a 

Finallyj  from  Equation  (3.28)  the  exponent  m0y 

as  the  slope  of  the  line  on  logarithmic  coordinates  of  the  value  of 
Qc  (or  C)  plotted  against  roller  diameter  Da.  The  tests  must,  in 
view  of  the  fact  that  Da/dm  appears  in  the  function  #,  be  run  for 

constant  values  of  Oa/dm. 

The  results  of  these  three  test  series  are  then  solved 
simultaneously  to  give  c,  h  and  e. 

t.  Line  Contact 

The  treatment  of  the  line  contact  problem  given  in  the  1947 
work  (l)  is  amplified  and  revised  in  the  1952  publication,  (2) 
which  is  devoted  exclusively  to  roller  bearings.  In  this  treat¬ 
ment-  «he  contact  load  is  foksn  to  be  propo rt 1  - " ” '  to  the  1.1 
pow-r  r  ■  ‘y.j  del ormat ior. ;  in  tn@  earlier  worn  ..  linear  relation¬ 
ship  was  assumed  for  line  contacts. 

It  was  found  that,  for  line  contact  the  exponent  p  in  the  load- 
life  relationship  is  4  rather  than  3  as  in  the  point  contact  case. 

It  is  possible  for  some  roller  bearings  to  have  point  contact 
within  one  range  of  loads  and  line  contact  within  another.  It  is 
©ven  possible,  under  some  conditions  in  a  roller  bearing,  for  the 
roller  to  make  line  coatset  with  one  raceway  and  point  contact 
with  the  other. 
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As  a  compromise  made  for  simplicity's  sake  a  conservai  ye 
compromise  exponent  of  p  --  ]0/3  is  assumed.  The  capacity  is 
thereby  modified  by  a  factor  which  is  calculated  so  that 
when  line  contact  prevails  at  both  contacts  the  error  in 
using  the  exponent  10/3  rather  *hnn  4  is  made  as  small  as 
possible  oyer  the  most  frequently  used  life  range. 

In  computing  the  capacity  of  roller  bearings!  a  reduction 
f actor^ ) i s  introduced  as  an  attempt  to  account  fcr  the 
stress  coscentrat i on  which  may  occur  at  roller  ends  as  well 
as  the  effect  of  inexactly  aligned  rollers. 

Another  reducti'  i  factor  ^<1  is  introduced  for  thrust 
loaded  bearings  to  account  for  the  effect  of  the  greater 
sliding  undergone  by  the  rolling  bodies.  In  thrust  leaded 
bearings,,  roller  loading  Is  not  cyclic,  but  remains  virtually 
constant,  exacerbating  the  problem  of  rolling  element  sliding 
and  ^action  forces. 

5.  FACTORS  OMITTED  FROM  LUNDBERG-PALliGREN  THEORV 

Lundberg  and  Palmgrenfin  their  pref ace,  acknowledge  the 
absence< from  their  development,  of  several  factors  known  to  affect 
endurance  life.  Specifically  cited  as  areas  for  future  investiga¬ 
tion  are: 

a)  Effect  of  stress  history 

b)  Work  hardening 

e)  Lubrleant  effect  on  stress  distribution 

d)  Effect  of  residual  stresses  (set  up  in  the  rings  by 
interference  fits) 

s)  Effect  of  edge  loading  in  line  contact. (2) 

f)  Effect  of  radial  load  on  contact  angle  ia  boll  bearings. 

g)  Effect  of  indexing  of  the  ball  rotational  axis  which, 
since  any  given  point  on  the  bell  is  cyclically  stressed 
for  only  part  of  the  bearing  life,  results  in  a  lesser 
number  of  failures  initiating  in  the  boll  material  than 
in  the  ring  material 

h)  Effect  of  surface  traction  (in  thrust  loaded  bearings) 

i)  Effect  of  geometrical  imp© rfeet i oni  on  load  distribution. 

Many  other  life  factors  are,  of  course,  absent  frem  the 
Lundberg  and  Paimgren  treatment  without  having  been  specifically 
enumerated  by  these  authors. 
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SECTION  IV 


7AEIABLES  AND  MECHANISM  OF  ROLLING  FATXGtJI 


1,  FAILURE  MECHANISMS  (GENERAL) 


.  B 


The  definition  of  functional  failure  in  rolling  bearings 
depends  on  the  application.  Except  in  instrument  bearings  in 
which  torque  is  often  the  main  criterion  in  defining  “failure8, 
the  most  common  definitions  of  functional  failure  involve  visible 
damage  to  the  rolling  tract.  As  an  example  of  specific  interest, 
the  failure  of  many  large  radar  antenna  pedestal  bearings  has 
been  found  to  fail  under  this  definition  involving  smearing, 
spalling  and  surface  distress  in  the  rolling  track.  From  these 
broad  concepts  of  failure,  several  groups  of  reasonably  well 
defined  changes  in  the  rolling  bearings  can  be  identified  which 
represent  failure  modes  (6). 


:/  1 
m  r 


The  failure  of  a  particular  bearing  is  a  consequence  of 
several  competing  modes  of  failures  classified  in  Table  1,  Aisong 
the  modes  of  bearing  failure,  the  present  study  pertains  to  the 
prediction  of  contact  f  at  lone  life  (mode  3  of  Table  1).  This  is 
justified  since  wear  and  plastic  flow  (failure  modes  1  and  2  of 
Table  l)can  he  eliminated  in  most  rolling  bearings  by  suitable 
design  and  operation  controls  whereas  cracking  (failure  raede 
No.  4)  ig  not,  specifically,  a  rolling  contact  failure.  On  the 
other  hand,  ail  loaded  rolling  contacts  are  eventually  subject 
to  fatigue  failure. 


ft?  • 


The  functional  failure  due  to  contact  fatigue  is  character¬ 
ized  by  the  emergence  of  a  fatigue  crack  causing  removal  of  a 
sizable  piece  f  metal  from  the  rolling  surface  (spalling). 

The  spalling  itself  is  preceded  by  the  initiation  and  propaga¬ 
tion  of  one  crack  out  of  possibly  several,  which  arise  in  the 
rolling  element  through  stress  cycling. 


The  present  study  considers  that  the  fatigue  process  in 
rolling  contact  is  associated  with  the  generation  and  propaga¬ 
tion  of  fatigue  cracks  in  the  bearing  material  and  the  fatigue 
life  is  taken  as  the  number  of  cycles  at  which  fracture  (spalling) 
occurs.  This  is  a  concept  taken  from  numerous  fatigue  studies 
most  of  which  are  not  for  rolling  contacts.  Fatigue  theories 
related  to  fatigue  crack  behavior  in  materials  can  be  divided 
into  two  distinct  approaches;  one  is  the  so-called  ^ngi  neeri  ng** 
(or  "phenosaeno  1  ogi  c  a  1 approach  (e.g,,  Manson  (7),  Coffin  (8) 

@nd  Morrow  (9)  on  low  cycle  fatigue  and  Dugdale  (10),  Paris  (11) 
on  sheet  specimens)  whieh  is  interested  in  quantitative  treatment 
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Table  1 


FAILURE  MODES  OF  ROLLING  CONTACTS 


1 „  Bear  type  fai  lures 


2 .  Plastic  flow 


3.  Contact  fatigue 


4.  Bulk  f oi lures 


1.1  Surface  removal 

1.1.1  P c ..I o v a  1  of  loose 
particles  l  "Wear  ’) 

1.1.2  Ciiealcsl  or  electric., 
surface  removal 

1.2  Cumulative  material  transfer 
between  surfaces  ("Smearing" 

2.1  Loss  of  contact  geometry 
due  to  cold  flow 

2.2  Destruction  by  material  soft 
ening  due  to  unstable  over¬ 
heating 

3.1  Spalling 

3.2  Surface  distress 

1.1  Overload  cracking 

4.2  Overheat  cracking 

4.3  Bulk  fatigue 

4.4  Fretting  of  fitted  surfaces 

4.5  Permanent  dimensional 
changes 
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of  macroscopii,  '/arnfiblei,  and  the  other  i  $  the  octal!  or gicai 
( r  r  microscopic)  approach  (e.g.  Laird  02),  Bond  (13)  and 
f.  .asekrtutz  (14))  which  deals  with  tho  basic  microscopic 
ph eiionenn  such  as  dislocations,  Kilobauds  and  microcrack 
initiation  in  materia]  under  cyclic  strain. 

Both  approaches  teach  that  plastic  occurrences,  either 
-macroscopic  (St«<  h  as  in  low  cycle  fatigue)  or  raicro.-copic  (in 
high  cycle  fatigue  (19))  are  a  source  of  crack  generation  sad 
propagation  under  fatigue  loading. 

Han y  attempts  Have  been  ,nade  to  relate  quantitatively  the 
fatigue  life  to  the  magnitude  of  cyclic  plastic  strain.  In 
published  literature,  quantitive  1 rsatment  has  been  made  by 
Hanson  st,  al.  regarding  crack  generation  and  propagation  in 
notched  specimens  (15)  without  taking  into  accou  t  the  micro¬ 
scopic  plastic  behavior.  The  prediction  of  the  life  for  a 
complete  machine  component  has  been  found  very  difficult  be¬ 
cause  of  the  complexities  in  geometry  and  loading.  Peterson 
in  (16)  has  studied  quantitatively  the  fatigue  problem  in 
turbine  components;  associated  with  cracking  due  to  thermal 
cycling  by  determining  cyclic  plastic  strain  cox  ,  «=  sp  Ouu  i  ng  n 
the  condition  of  operation,  taking  account  of  stress  concen¬ 
trations.  Sice  snd  Brown  ( 4 i >  have  attempted  to  interpret  the 
fatigue  fracture  of  loaded  structural  elements  in  terras  of 
fatigue  crack  propagation  and  to  forrnuiate  the  fundamentals 
of  a  statistical  theory  in  these  terms. 

Quantitative  treatment  as  described  above  has  to  rely 
on  many  assumptions  and  phenomenological  conclusions  from 
specimen  testing.  The  basic  physical  mechanisms  of  crack 
generation  and  propagation  in  materials  under  fatigue  loading 
is  of  course  a  subject  of  much  fundamental  research  in 
metallurgy,  Sore  specifically,  cracking  st  non-sseta  i  1  i  c 
inclusions  has  been  g  subject  of  considerable  interest  in 
basie  raetal  lurgieal  studies  but  quantitative  treatment  correla¬ 
ting  various  macroscopic  quantities  has  not  been  available 
for  engineering  design  application. 

Prom  the  reviewed  literature,  the  following  can  be 
gleaned: 

Iu  spite  of  the  great  physical  complexity  involved  in 
the  process  of  crack  initiation  from  inclusions,  it  is  well 
recognised  that  non-set  si  1 ic  inclusions  are  stress  (or  strain) 
rai’sers  and  potential  sites  of  plastic  deformation.  The  above 
described  literature  pointsto  a  close  association  between 
fatigue  life  and  the  magnitude  of  the  cyclic  plistic  strain. 

It  ©ay  be  poiiiblg  to  build  a  fatigue  life  -lodel  by  predicting 
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the  number  of  strege  cycles  needed  to  generate  a  self- 
propagating  crack  ss  d  function  of  the  magnitude  of  (cyclic) 
plastic  strain  at  the  stress  raisers.  The  plastic  strain  is, 
in  turn,  dependent  on  the  macro-strain  field  and  the  inherent 
characteristics  of  the  stress  raiser.  Such  an  engineering 
approach  has  been  used  by  Benson  £-  Hirechbcrg  (15)  and  by 
Peterson  (16)  to  relate  crack  initiation  life  to  the  stress 
concentration  factor  at  a  notch. 

The  quantitative  dependence  of  fatigue  life  on  fatigue 
ductility  and  plastic  strain  ccn  be  described  as  follows: 

Coffin  (0)  and  Manson  ( 2t )  has  investigated  the  fatigue 
behavior  in  a  wide  range  of  materials  by  strain  cycling  at. 
constant  strain  amplitude  at  out  zero  mean  strain.  In  the  low 
cycle  region  where  failure  occurs  in  10®  cycles  or  less,  the 
fatigue  life  Is  found  to  be  related  to  the  plastic  strain 
amplitude  by  an  equation  of  the  fol. owing  form  for  ail  materials 


where  A§r 


B  and  g 


As  =  a  (  2N  ) 

P  f 

-  total  plastic  strain  amplitude 
=  number  of  cycles  to  failure 
=  constants 


(4.1) 


iaason  (15)  has  ©steaded  the  treatment  of  strain  cycling 
data  to  cover  the  ©at ire  range  of  fatigue  lives,  from  the  low 
cycle  region  where  strains  ©re  predomi nantly  plastic  to  the 
high  cycle  reqions  where  stains  « 'c  predominantly  elastic. 
Maason8§  equation  relating  total  strain  amplitude  sad  cycles 
to  failure  is  given  by: 


Agf  =  +  AS  = 


(2N/)‘ 


L  (2N/) 


(4.2) 


where  a® 


T 


-  total  strain  amplitude 
s  plastic  strain  amplitude 
Ai  =  elastic  strain  amplitude 


B  and  z  =  the  sgrae  constants  given  in  Equation  (4.1) 
L  s  a  constant,  (see  Hgure  03) 
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Totol  Strain  Amplitude ,  A  €  T080! 


where  E  -Young’s  modulus 

and  y  =  a  constant,  (kgp  Figure  .tii 

This  general  equation  covering  the  entire  range  of  cycling 
lives  is  shown  in  Figure  .13  and  represents  a  curve  which  can  be 
expressed  as  the  sum  of  the  two  straight  line  components,  the 
steeper  component  to  the  left  represents  the  plastic  strain  component 
while  the  elastic  strain  component  is  represented  by  the  line  with 
shallow  slope. 

The  fatigue  behavior  of  a  material  can  therefore  be  charac¬ 
terised  by  the  four  constants  of  Equation  (4.2)  i.e.  H,  z,  L,  ■ 
and  w  . 


Tli e  plastic  portion  of  the  curve  reflects  the  matrix 
properties  of  the  material  while  the  fact  that  the  elastic 
portion  a£  the  curve  slopes  at  all  is  due  to  the  presence  of 
defeets.  This  latter  conclusion  follows  when  one  considers 
that,  in  purely  elastic  reversals  of  strain,  nothing  changes 
in  the  metal  to  cause  failures.  Since  failure  does  occur  in  the 
nominally  elastic  stress  range,  it  must  be  due  to  localized 
yielding,  probably  in  the  vicinity  of  defects. 

Thus,  failures  in  the  elastic  region  are  defect  dominated. 
The  plastic  region  of  the  curve,  however,  is  derived  from  tests 
where  large  plastic  strains  develop  throughout  the  bulk  material. 
While  it  is  true  that  defects  will  exert  some  influence  on  this 
portion  of  the  curve,  properties  of  the  matrix,  e,g„  strain 
hardening  and  ductility,  will  tend  to  dominate  the  results.  The 
predominance  of  matrix  properties  in  the  plastic  region  has  been 
confirmed  by  matrix  properties  in  the  plastic  region  and  has 
been  confirmed  by  Morrow  (49).  He  states  that  little,  if  any, 
difference  can  be  found  between  clean  and  dirty  steels  of 
similar  processing  and  composition  when  tested  in  the  plastic 
strain  region. 

Tnus  the  plastic  portion  o£  the  curve  is  suited  to  describe 
the  ductility  properties  of  candidate  materials.  The  intercept 
of  the  plastic  strain  line  at  2Nj.=  i  (M  0f  Equation  (4,2))  nnd 
the  slope  of  the  line  (z)  are  considered  fundamental  fatigue 
properties  by  Sanson  and  Morrow.  They  have  been  given  the  names 
’’fatigue  ductility  coefficient”  and  "fatigue  ductility  exponent" 
by  Morrow  (16). 

It  le  apparent  from  inspection  of  Figure33  that  an  increase 
in  the  intercept  value  (M),  lor  a  fixed  slope  (z),  will  trans¬ 
late  the  line  upwards,  giving  increased  fatigue  lives.  Similarly 
a  decrease  in  the  slope  at  a  fixed  intercept  will  rotate  the  line 
upwards  with  the  §arae  effect  on  life. 
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Many  variables  have  h  e  e  n  found  to  affect  fatigue  life  of 
rolling  hearings.  Based  on  the  current  knowledge  of  the  opera¬ 
tion,  lubrication,  metallurgy  and  theory  of  fatigue  fracture 
of  rolling  bearings.  Table  2  presents  nn  ex  ten  si ve  list  of 
Hexfcrnai”  variable'  whiub  cor.tio1  the  fatigue  life  of  a 
roiling  contact.  These  variables  can  be  grouped  into  categories 
as  shown  below; 

a.  Material  Variable." 

1)  Factors  that  effect  the  material  yield  strength 
and  ductility;  material  analysis,  hardness,  soft  con¬ 
stituents  (retained  austenite,  ferrite,  bain lie),  grain 
size,  alloy  segregation. 

2)  Factors  that  modify  the  applied  stress  field; 
residual  stresses  originating  from  beat  treatment, 
grinding  and  plastic  flow  during  operations. 

3)  Material  imperfections  acting  as  stress  raisers 
such  os  non-metal lie  inclusions  or  imbedded  micro-cracks. 
(Lenticular  carbides  which  develop  during  bearing  opera¬ 
tion  may  also  serve  as  localized  stress  raisers.) 

4)  Modulus  of  elasticity  and  Poisson  ratio,  as  they 
affect  materi"l  rigidity  and  control  the  maximum  stress 
level , 

b.  Surface  Si c ro-Gcometry  Variables 

1)  Surface  imperfections  such  as  grinding  furrows, 
scratches  and  dents  which  serve  as  surface  stress  raisers. 

2)  Genera.1  surface  roughness  as  induced  by  methods 
sueh  as  grinding,  honing  and  lappiBg,  and  os  characterized 
by  a)  amplitude,  e,g.t  the  composite  r.m.g.  surface  rough¬ 
ness  defined  as  the  square  root  of  the  sum  of  the  squares 
of  the  r.m.s.  roughnesses  of  two  surfaces  rolling  together 
and  b)  a  plasticity  parameter  of  the  asperities,  such  as 
their  typieel  slope. 

3)  Compositional  and  hardness  properties  of  the 
surface » 

4)  Surface  coatings. 
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c.  Benign  Variables  Meiated  to  Design  Dimensions 

1)  Moiling  body  design  such  as  rolling  track  length 
(pitch  diameter),  groove  conformity,  ball  (roller)  diameter 
end  number  of  rolling  elements,  contact  angle,  roller 
length  and  crowning  and  1  lending  radii,  accuracy  parameters 
influencing  load  distribution,  dimensions  controlling  sliding,, 
etc.  These  variables  affect  stress  distribution,  stressed 
volume  and  number  of  stress  cycles, 

2)  Cage  and  auxiliary  par  s  design  as  it  affeets 
rolling  element  forces, 

d.  Operating  variabi  s  s»:h  s  s; 

1)  Logo  magnitude  and  direction  as  it  affects  stress 
level  flrd  stressed  volume, 

2)  Speed,  as  it  affeets  lubrication  (e.g„,  EHD  film 
thickness), 

3)  Temperature,  as  it  affects  lubrication  and  aaterial 
strength . 

4)  Lubricant  properties,  a)  rheological  properties, 

e. g.  (1)  viscosity  and  (2)  boundary  lubricating  ability 
(chemistry,  additives,  etc.) 

5)  Atmospheric  conditicns,  contaminants,  etc. 

The  above  listed  variables  are  believed  to  be  applicable 
for  rolling  bearings  of  all  types  and  sizes,  including  large 
radar  antenna  pedestal  bearings. 


3,  MECHANISES  OF  FAILURE  IN  ROLLING  CONTACT 


The  present  study  is  based  on  the  concept  that  there  are 
two  competing  fotigue  mechanisms  operating  to  promote  spelling 
failure  in  rolling  contact,  namely,  sub-surface  initiated 
failure  and  surface  initiated  failure.  In  the  former 
the  cracks  are  generated  from  highly  stressed  sub-surface  regions, 
whereas  in  the  latter,  fatigue  cracks  are  generated  from  the 
rolling  surface.  Beeent  findings  of  metallurgical  investigations 
(6,  17,  16)  gupport  this  classification  of  failure  mechanisms. 
Structural  changes  have  been  observed  in  bearing  material  under 
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repealed  loading  117,  10).  These  structural  changes  encompas- 
Cl)  plastic  deformation  around  sub-surface  weak  points  (e.g., 
'"butterflies”),  which  ind.cstes  that  the  mi  c  ro-de  f  cc  t  s  act  as 
stress  raisers,  (2)  the  formation  of  plastic  deformation  bonds 
(Hwhite  etching  areas'5)  in  the  sub-surfnee  high  shear  stress 
zone  signifying  the  existence  of  a  region  within  which  the  yield 
limit  is  exceeded  for  some  volumes  and  £3)  near-surface  micro- 
plastic  occurrences  oi  asperities  and  at  surface  defects. 

Among  the  structural  changes  described,  the  formation  of 
'butterflies"  ernund  inclusions  has  long  been  closely  associated, 
with  sub-surface  crack  generation  (1,  16).  A  similar  associa¬ 
tion  between  structural  changes  and  cracking  has  recently  been 
demonstrated  for  the  near  surface  changes  (3  above)  and  with 
some  indirection,  for  the  generalized  changes  (2  above). 

The  following  summarizes  current  knowledge  of  the  two 
fstigue  failure  mechanisms  described? 

a.  Sub-surface  Failure 

The  concept  of  sub-surface  failure  is  well  covered  in 
the  Lundherg-Pa  ling  ran  theory  (1)  in  which  the  fatigue 
crack  is  assumed  to  start  from  week  points,  e.g.,  slag 
inclusions,  in  current  terminology,  these  weak  points 
give  rise  to  the  local  stress  concentrations  and  plastic 
flow  in  the  surrounding  matrix  material.  According  to 
tundberg-f alragren  the  site  of  crack  generation  is  the 
zone  of  high  shear  stress  in  the  sub-surface  region  of 
a  roiling  element.  Fatigue  cracks  are  found  to  start  at 
weak  volumes  and  to  grow  under  repeated  loads  until  an 
advanced  stage  of  fatigue  cracking  is  reached.  This  stage 
is  characterized  by  the  distortion  of  the  macroscopic 
stress  field  due  to  the  cracks.  Eventually  the  destruc¬ 
tive  process  of  spalling  sets  in  at  one  (or  more)  loca¬ 
tions  causing  removal  of  a  sizeable  piece  of  metal  from 
the  surface. 

It  is  recognized  that  the  formation  of  localized 
plastic  deformation  around  inclusions  (i.e.,  the 
"Butterfly55  structure)  is  stress  dependert.  In  (1), 
Lundberg  and  Palmgren  hypotha s i zed  that  the  actual 
stresses  o£  nub-surface  '’weak"  points  of  the  material 
are  proportional  to  the  magnitude  of  the  macroscopic 
stress  and  are  modified  by  factors  dependent  on  the  size 
ind  microscopic  ships  of  each  of  th©  "weak  pointo" 
(inclusions),  lecent  met  all urgi e@ i  1 nvestigeti on  has 
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supported  this  hypothesis,  g.g,,  Liftman  end  Midner  CIS) 
pointed  out  that  each  n  on-met  si  1  i  c.  inclusion  has  its  own 
stress  concentration  factor  which  depends  on  the  size, 
shape,  physical  end  mechnnict!  properties  o£  the  inclusion; 
Martin  and  Eberhardt  (IT)  have  shown  that  '’Butterfly" 
structure  can  occur  at  locations  where  tho  calculated 
"macro”  stross'3  are  below  the  threshold  level  of  plastic 
flow  indicating  chat  inclusions  are  stress  raisers. 

b„  Surface  Distress  and  Surface  Initiated  Fatigue 

The  above  described  sub-surface  failure  will  occur  for 
any  lubrication  and  for  smooth  surf sees  (thick  elasto- 
hydrody Bamic  films).  Under  such  conditions,  life  to  failure 
appears  to  be  independent  of  both  lubricant  end  surface 
texture.  There  is,  however,  an  altogether  different  mode 
of  rolling  contact  fatigue  referred  to  as  surface  distress 
(6,  20).  Although  surface  distress  does  not  imply  raceway 
destruction,  it  is  a  precursor  and  often  a  precipitator  of 
spoiling  failure,  apparently  by  generating  severe  surface- 
adjacent  defects  (about  100  pin.  depth)  which  then  serve 
as  crack  initiation  points  (6).  Metallurgical  evidence 
has  shown  that  this  s„urface  distress''  involves  a  type  of 
plastic  working  end  subsequent  fatigue  microcracking  of 
the  immediate  surface-adjacent  layers  of  the  metal.  The 
occurrences  of  near- surface  plastic  deformation  are  of  the 
following  two  kinds: 

a)  Hide-spread  microp last xc  flow  at  surface  asperities 
caused  by  asperity  interaction; 

b)  localized  plastic  flow  at  surface  imperfections, 
such  as  grinding  furrows,  scratches  and  debris  dents. 

Begarding  the  occurrence  of  near-surface  plastic  defor¬ 
mation,  lallian  (6)  has  suggested  that  it  is  controlled  by 
the  elastohydrodynaaie  lubricant  film  and  by  surface  rough- 
ness  parameters,  e.g.,  the  composite  surface  roughness  r.m.s. 
end  a  typical  asperity  slope  angle.  In  (6)  it  is  hypothe¬ 
sized  that  the  near  surface  plastic  occurrences  are  a  conse¬ 
quence  of  severe  interaction  of  asperities.  The  approach  of 
the  contacting  asperities  depends  on  the  mean  EUD  film  thick¬ 
ness  defined  as  the  distance  between  the  mean  line  of  the  two 
asperity  profiles.  The  degree  of  approach  determines  the 
severity  of  the  plasticity  in  the  asperity. 

It  has  also  been  reported,  based  on  recent  investiga¬ 
tions  of  failed  bearings  (IT),  that  for  rolling  bearings 
mode  of  improved  clean  steel  (obtainable  by  the  vacuum 
melting  process),  a  larger  fraction  of  failures  is  found  to 
be  tsaocioted  with  the  surface  defects.  This  if  because  of 


fewer  and  smaller  subsurface  non-set ni 1 ic  inclusions  being 
available  to  initiate  fttigue  cracks  in  clean  steel;  thus 
the  role  of  surface  defeets  becomes  more  important. 


4„  FAILUEE  PIOCESS  OIAGlAi  (BOLLING  CONTACT) 


Taking  into  account  the  previously  described  failure 
variables  and  failure  mechanisms,  it  is  possible  to  dras  a 
diagram  of  the  roiling  contact  failure  process,  Figure  1 
shows  this  diagram  which  depicts  tiie  i nterdependence  of  the 
v  ar tables  and  their  effects  on  subsurface  and  surface  initiated 
spalling  occurrences.  This  chart  is  a  brief,  self-explanatory, 
illustration  of  the  previously  described  failure  process. 

The  interacting  effeets  of  the  variables  are  shown  in  captions 
enclosed  in  the  rectangular  blocks  of  Figure  1.  The  arrows  counted 
en  the  lines  signify  an  effect. 

The  "terminations"  of  this  flow  chart,  shown  on  the  right 
hand  side  of  Figure  1,  represent  1)  subsurface  spalling  occur¬ 
rences  controlled  by  subsurface  crack  initiation  within  a 
stressed  volume  and  2)  subsurface  spalling  occurrences  con¬ 
trolled  by  the  surface  crack  initiation  within  a  stressed  area. 

The  stressed  area  end  volume  are  determined  from  the 
quasi-elastic  subsurface  and  surface  stress  field.  Subsurface 
(or  surface)  crack  initiation  is  affeeted  by  subsurface  (or 
surface)  weak  points  and  the  magnitude  of  subsurface  Cor 
surface)  plastic  atrainjhis,  in  turn,  dependents  on  the 
quasi-elastic  stress  field  and  the  material  yield  strength 
which  have  the  effect  of  controlling  the  amount  of  plastic 
deformation  for  a  given  stress  field. 

Dactility  is  believed  to  have  the  ability  to  suppress 
crack  initiation  and  propagation  in  a  material  matrix  under 
a  given  cyclic  raacroslraiu  field.  It  is,  of  course,  a  matrix 
parameter. 

In  Figure  1,  the  variables  listed  on  the  left  side  consist 
of  two  kinds,  i,e,,  1)  material  impurities  and  surface  imper¬ 

fections  acting  as  weak  points  and  2)  the  variables  affecting 
the  macros  tress  field  in  the  material. 

The  following  describes  the  external  variables  which  affect 
the  aocrostress  field  in  bearing  materials 
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a.  The  effect  of  design  variables  ( diaens 1 ons  ) ,  elastic 
modulus,  and  load  on  the  elastic  isscrostress  field  has  been 
sell  covered  in  the  Lundberg-Palagren  theory  (1>  based  on  the 
application  of  Hertzian  theory  (32). 

b.  Residual  stress  serves  as  a  stress  modifier, 

c.  The  Lubricant  properties,  speed  and  temperature  at 
the  rolling  contacts  are  parameters  which  have  been  found 
in  recently  developed  e 1  a s t ohydr ody nami e  theory  to  affect 
EHD  film  thickness.  Soiling  contact  experiments  have 
shown  {48)  that  the  EHD  film  thickness  and  surface  finish 
(characterised  by  surface  roughness  iSS  and  asperity  slope 
angle)  have  significant  influence  on  the  occurrence  of 
surface  distress  or  n@§r~aurf ace  micro-plasticity  (6,44). 

d.  The  effect  of  speed  i.e.  strain  rate  (or  the  time 

a  material  element  is  under  stress)  and  temperature  in  the 
highly  stressed  volume  (or  surface)  on  material  strength 
is  hypothesized  on  the  basis  of  non-rolling  contact  experience. 


SECTXQK  V 
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Both  sub-surface  and  surf  gee-  failure  aeck  ani e@§  require  a 
crack  initiation  phage  followed  ty  §  crack  propagation  phase  to 
reach  the  point  of  a  functional  failure  (spelling)  of  a  rolling 
element  (6). 

The  crack  initiation  in  material  under  cycling  stressing  J.g0 
in  general  (except  la  §herply  notched  specimens),  a  process  re¬ 
quiring  sany  cycles  of  stressing.  The  necessary  conditions  for 
cr*ck  initiation  are  generated  by  plastic  deformation  ground 
material  ishemogeaeities  or  stress  raisers  on  a  result  of  cyclic 
stressing.  The  crack  initiation  phase  can  bo  defined  as  that 
isvolviag  th©  formation  of  cracks  oa  the  microscopic  scale.  It 
is  usual  to  consider  that  this  phase  is  ended  when  a  self-propa¬ 
gative  crack  emerges  at  the  site  of  crack  initiation,  i.e.,  one 
which  propagates  without  further  influence  of  the  original  defect. 

An  analysis  by  Mangos  and  Hirschberg  (15)  postulates  that 
the  tine  to  the  initiation  of  a  fatigue  crack  depends  on  the 
magnitude  of  the  cyclic  plastic  strain  at  the  strain  raiser. 

The  strain  raiser  considered  in  (15)  is  a  macroscopic  one,  i.e., 
a  notch  in  a  specimen.  As  usual,  the  fatigue  process  is  divided 
into  two  stages,  i.e.,  crack  initiation  and  crack  propagation, 
k  certain  crack  size,  called  "engineering  site",  is  used  to 
demarcate  the  two  stages.  It  wag  found  that  when  an  ^engineering 
size55  ersck  emerges,  at  the  notch  the  controlling  strain  for 
further  crack  growth  no  longer  depends  on  the  strain  eoBeeatra- 
tlos  factor  at  the  notch  but  os  the  nominal  strain  in  the  sp©ei= 
B@ne  Because  the  strain  range  is  relatively  high  at  the  notch, 
the  role  of  aiere-defects  can  be  neglected.  The  concept  ia 
this  analysis  is  useful  in  understanding  crack  initiation  in 
rolling  contact,  although  it  i§  realized  that  ia  rolling  contact 
oat© rial  the  strain  raisers  sra  actually  all  microscopic. 

In  rolling  contact  the  aiero-defeets  can  be  treated,  for 
convenience,  as  strain  raisers,  Little  is  known  about  the 
cyclic  strain  concentration  near  a  sderostreas  raiser,  that  Is, 
for  a  given  real  inclusion,  there  is  no  adequate  means  to  know 
either  by  observation  or  by  analysis  the  strain  raising  factor. 
Metallurgical  experience  (3)  has  shown  that  the  size  of  an  in¬ 
clusion  has  a  large  effect  oa  ersck  initiation  life;  other 
variable!  sueh  as  shape  and  content  aay  also  have  an  effect  oh 
crack  initiation.  The  size  and  shspe  of  the  defects  presest  in 
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is  any  volume  of  metel  are  t-ioehastic  quantities.  This  calls 
for  a  statistics!  treatment  in  whieh  the  cyclic  strain  rsising 
factor  or  severity  (designated  as  d)  of  a  defect,  §§  a  function 
of  its  size,  shape  and  contest,  is  introduced  as  a  random  variable. 

In  formulating  mathesati ea 1  models  of  rolling  contact  fatigue 
for  the  surface  and  sub-surface  failure  mechanisms,  it  is  hypo¬ 
thesized  that  the  fatigue  cracks  in  both  mechanisms  initiate  from 
micro-defect  %  or  nt?eak  points".  In  the  sub— surface  failure  mecha¬ 
nisms,  the  micro-defects  are  embedded  i  the  matrix  material 
eiereag  ia  surfsee  failures  pre-exi sti ng  surface  micro-defects 
and  surface  fatigue  cracks  are  considered  as  the  weak  points  of 
§  surface  eleeeni.  This  spre=e*ci  iti  n@  micro-defect  hypothesis'5 
is  postulated,  of  eoune9  only  for  hard  steels  as  used  in  bearings. 

Current  concepts  of  metal  fatigue  assert  (22)  that  the 
controlling  pargeeters  of  fatigue  life  are  the  applied  eyclic 
strain  §ad  material  ductility  g§  they  influence  crack  growth. 

TM©  rate  of  growth  ©f  a  crack  is  denoted  by  dA/dM  where  A  is 
the  i Rstaataaeoua  crack  area  and  ft  is  the  number  of  stress  cycles, 
it  is  hypothesized  thatlhe  growth  rate  dA/dM  of  g  crack  Is  depen¬ 
dent  os  the  plastic  mierostraia  <0  aad  ductility,  D,  prevailing 
la  tke  vicinity  of  a  “defect53  consisting  of  a  “pre-exi  st  i  ag 
origlsal”  defect  aad  of  the  crack  initiated  therefrom. 

Tfeug  a 

34 

it  =  A  <  ,,  ,  0  ) 

The  ductility  0  Mas  been  defined  is  static 

0  =  l©g(A©/Ap)  =  lag  y_ g _ a / 

sfeere  A e  §ad  A#  are  the  initial  aid  fis§I  area  of  the  fracture 
erssi  section  ia  the  tengile  test,  aad  I„A.  is  the  convent! onal 
redaetisa  iu  ares. 

However,,  there  exist  other  ductility  related  quantities 
which  cau  be  defiacd  from  fatigue  tests.  It  was  shown  by  Coffin 
(8)  that  a  straight  lino  results  when  plotting  the  logarithm  of 
cycles  to  failure  against  the  logarithm  of  eyelic  reversing 
plastic  strain  ia  specimen  fatigue  tests,  i.e., 

&  sp  =  M  (2  fif)  (5.3) 


(5.1) 

test  as  (21); 

(5.2) 
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tsfeeee  the  reversed  cyclic  plastic  atraia 

N /  •  number  of  cycles  to  failure 

S  :  a  material  constant  called  the  fatigue  ductility 
coefficient  (9) 

2  ;  s  constant 

Thuss 8  and  z  are  ductility  related  constants  which  can  be 
determined  from  specimen  testa  of  material.  The  above  iatro=- 
duced  quantity  0  can  be  considered  as  a  function  of  I  and  z. 

It  is  keosa  (23)  that  increasing  hydrostatic  compression 
increases  material  ductility.  in  the  above  formula,  (Equa¬ 
tion  5.11j  the  beneficial  effect  of  compressive  hydrostatic 
stress  can  be  taken  into  account  by  the  use  of  the  ductility 
0  existing  under  the  given  stress  conditions. 

The  hypothesis  expressed  in  Equation  (5.1)  takes  account 
of  the  facts  that  £1)  microplastie  deformation  is  a  criterion 
for  fatigue  crack  generation;  and  for  a  given  ductility,  the 
crack  generation  rate  increases  with  the  magnitude  of  cyclic 
plastic  strain  amplitude  at  the  stress  raiser  and  (2)  the 
generation  or  the  groath  of  a  crack  is  controlled  by  the 
material  ductility,  i.e.„  highly  ductile  material  craeks 
less  rapidly. 

It  is  assumed  that  there  exists  a  macroscopic  strain  eQ  at 
the  location  of  the  defeet  if  the  defect  were  absent  so 
that  8  quantity  called  '’defect  severity”  can  be  defined  as 
the  strain  raising  factor  operating  oa  eQ  to  give  the  strain  at 
the  (subsurface  or  surface)  defeet.  8,  the  severity  of  a 
combined  defect  and  craek  is  hypothesized  as  a  function  of  the 
original  defect  severity,  d,  the  instantaneous  crack  area 
h  and  the  size  of  the  stressed  volume  (for  sub-surface  fatigue) 
or  area  (for  surface  fatigue)  §„  i.e.s 


©  =  0(  d,  A(N),  S  )  (5,4) 

lb©  severity  of  localised  plastic  siroia  concentration 
around  the  defect  (as  demonstrated  e.g.  by  KiatterflyB 
structure )  is  dependent  ons  1)  the  strain  raising  properties 
of  the  defect,  ©  2)  the  macrostrein  s  assuming  the  strain 

raiser  is  absent  and  3)  the  yield  strength  <?y  of  the  matrix 
material „  i ,  e, 

§e  =  §c  (gg, a  Q) 

\MS  5  ) 

®  i@  is  defined  as  total  (elastic  plus  plastic)  mac ros t r §i a . 
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Introducing  Equations  (5.4)  and  (5.5)  into  Equation  (5,1)  one 
has: 


tlA 

dN 


(5.6) 


It  is  convenient  to  separate  the  variables  influential  in  crack 
growth  into  two  groups:  variables  related  to  defects  and 
matrix  variables.  In  Equation  (5.6).  ©is  the  variable  related 
to  defects  wher@ggs0,ay  and  D  are  related  to  the  matrix.  For 
simplicity,  the  matrix  effects  are  c on g ol i ri ated  into  a  eingle 
function  Y  i.e, ; 

M 

dN 


=  A  (©,  y). 


V  -  ( e 


O'  °y*  D) 


(5.7) 


Assume  that  the  crack  growth  Tate  dA/dS  can  be  expressed  as  a 
product  of  he  severity  function  and  the  matrix  function,  i.e. 


M  =  @.Y 

dN 


(5.82 


Using  this  designation,,  the  defect  function  plays  the  role  of 
a  Hpreporti onal 1 fcy  constant*  between  the  crack  rate  and  the 
matrix  function  of  the  raaeroetrain  e0  the  ductility  D  and  the 
micro-yield  strength  oy  .  This  proportionality  constant  is, 
of  course,  sa  inherent  quality  of  §  given  defect. 

la  the  eraek  iaitistlon  phase,  the  effect  of  stressed 
area  S  on  ©  can  be  neglected.  The  fuaetion  §j  can  be 
written  as  a  product  of  a  fuaetion  of  d  designated  byF(d) 
i  function  of  A„  designated  byfj(A).  Thus  one  has: 


©I  =  fi(A)- r(d) 


(5.9) 


and 


dA 


-  fj(A)-r(d).  Yj  (oy,D,  e0) 


meters  ey , D 


In  the  crack  initiation  period,  the  matrix  material  para= 
and  macro-strain  s0  can  be  considered  to  be 
constants  with  cycling.  Letting  y^  be  independent  of  N, 
Equation  (5.9)  eae  be  further  written  as  follows  after  inte¬ 
gration  with  respoet  to  Ms 


VA) 


Nj'  Vj-r^d) 


(5. 10) 
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. . . . . . mill . . 


NUMBER  OF  CYCLES,  N 


etjematic  lepreseat ati oa  of  @rot?fh  of  liero 


A  quantity  ft  y 
required  to  reach  a 


representing  the  iiumbei  of  stress  cycles 
crack  size  may  be  define*.'  fromr 


f!<V  :  NI  V(d) 


(5.11) 


where  fT  represents  a  function. 

Since  Ap,  the  self  propagating  crack  size,  is  a 
constant,  Equation  (5.10)  implies  that  for  &  given  defect  with 
severity  Hda  existing  in  the  highly  stressed  volume,  there  is 
a  correspond! ng  life  Nj  to  crack  generation,  in  numbers  of 
cycles,  associated  with  this  defect.  Thus,  Equation  (5,6)  can 
be  solved  for  Nr*. 


f I ( Ap  ) 
1  ~  Vi  TCd) 


(5.12) 


The  function  r ( d )  can  be  assumed  to  be  &n  increasing 
function  of  d,  thus  increasing  d  corresponds  te  decreasing 
N j .  Equation  (5.12)  implies  that  every  defect  has  a  life 
associated  with  its  defect  severity  for  8  given  matrix 
factor  Yj.  Figure  2  shows  a  schematic  representation  of 
growth  of  micro-cracks  from  defects  of  various  severities 
and  the  dependence  of  cycles  to  self -propagating  crack  genera¬ 
tion,  on  original  defect  severity  d  8S  defined  in  Equation 
(5.12),  provided  tha,t  the  matrix  parameter  is  the  same  for 
each  case.  The  graph  shows  that  increasing  d  corresponds 
to  increasing  slopes  of  the  curves  plotting  A  against  N 
(i.e,,  the  crack  rate)  and  shorter  life  to  crack  generation. 

The  degree  of  concentration  of  the  plastic  strain  around 
an  original  defect  diminishes  the  distance  from  the  original 
defect.  One  can  define  a  macroscopic  size  of  the  crack  beyond 
srhieh  the  crack  propagation  rote  is  dependent  only  on  the 
plastic  strain  at  the  crack  tip  and  no  longer  on  the  plastic 
strain  around  the  original  defect.  (An  analogy  of  this  case 
is  given  in  ianson9s  fatigue  life  analysis  of  notched  speci¬ 
mens  in  US)).  Thug,  for  large  Ae  it  can  be  assumed  that  the 
severity  of  a  "defect  after  N  cycles”  is  dej endent  on  the 
instantaneous  crack  area  A  and  the  applied  strain  e0 but  is 
independent  of  the  original  defect  severity  d.  The  propaga- 
tion  of  a  crack  of  this  size  constitutes  Phase  1J  of  the 
fatigue  process.  For  a  large  crack  area,  an  increasing  ratio 
A/S  corresponds  to  an  increase  of  stress  concentration  at  the 
crack  tip.  For  such  craeks  dA/iK  increases  with  Vs  .  Taking 
into  account  this  size  effect,  kb  assume  that  the  severity 
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function  in  Phage  II  is  &  product 

and  a  function  of  A  g@  follows; 


g  inaction  ef  4 


-  6  j  ^  ( A  ,  S  )  .yt  j  (  ^ ,  D,  Oy  ) 


(5,13) 


drr  =  fg  (A/S)  •  f.  (A) 


(5, 14) 


where  the  use  of  the  function  >tt  denotes  the  fact  shit  this 
enaction  may  be  different  for  the  propagation  and  initiation 
pha  geg. 

Integrating  Equation  (5.13)  gives: 


fC  dA 

“A-,  ?rrn  f. 


(5.15) 


where  nos  represents  the  number  of  cycle*  measured  from 
the  starting  fitae  when  A  =  Ap,  Equation  (5.15)  becomes  upon 
i ntegrati on  . 


fn(A,  A/S)  =  Yn 


(5.16) 


Tho  crack  arcs  reaches  critical  six®  Ac  after  *31  cycles 
in  Phase  II  where,  from  Equation  (5.16)  2 


Nil  =  f  1 1  <A,Ac/S)/yu 


(5.  17) 


The  total  number  of  cycles  required  to  produce  a 

creek  of  size  Ac  is  ; 


NL  =  Hi  +  Nfi 


(5.18) 


Neglecting  any  additional  cycling  require!  in  Phase  III  for  a 
crack  to  grog  from  size  Ac  to  produce  failure,  represents 

the  life  in  cyeles  associated  with  a  defect  of  severity  d.  It 
is  not  yet  known  whether  the  amount  of  life  thus  neglected  may 
be  appreciable. 
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Oil  n  g  Equations 


C S e 1 2 > s  («i.i7)e  iquitien  (5,16)  becomes 


f 

!  • 

n 

y 

i 


fjUp)  L  f  n(Ac,  Ac/Ag) 
Yj ' f( d )  Yl 


(5. 19) 


Equation  (5.19)  in  purely  deterministic  in  natures  that  is, 
if  the  functions  are  all  known  and  the  ductility  and  stressed 
area  are  prescribed,  then  Equation  (5.1§)  will  jioid  the 
number  of  cycles  that  elapse  until  a  crack  emanating  from  a 
defect  of  designated  severity  d  situated  in  the  stress  field 
to  as  to  be  subjected  to  a  strain  «0(  achieves  the  area  Ac  , 


It  has  been  shown  in  non-rolling  fatigue  testing  (22) 
that  the  erack  initiation  phase  has  higher  dependence  on  applied 
stress  (or  strain)  than  the  craek  propagation  phase,  which  implies 
that  the  increase  of  stress  level  will  shorten  the  crack  initia- 
tion  phase  sore  rapidly  than  the  crack  propagation  phase  <  (24), 
(25)),  Thus  the  Phase  I/Phase  II  ratio  decreases  with  increasing 
stress  levei. 


The  formation  of  a  spall  from  a  critical  sized  crack  is 
considered  a  rapid  process  associated  with  cleavage  and  dimple 
rupture  (27),  It  is  recognized  (26)  that  final  fracture  is  a 
complex  subject  in  itself,  especially,  when  large  plastic  strain 
is  involved.  It  is  assumed  that  the  critical  crack  size  Ac  is 
related  to  the  size  of  the  stressed  area  S  in  a  rolling  element 
cross  section; 


A„  =  ktg.  &i  '  const.  (5.20) 

It  ig  assumed  that  after  reaching  a  critical  crack  size 
the  fracture  enters  a  new  phase,  i.e.,  Phase  III,  which  involves 
a  relatively  sapid  cracking  (spalling)  process,  resulting  in  a 
visible  spall  on  rolling  surface. 
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SECTION  VI 


STATISTICAL  TfiSOI?  OF  HQLlfNG 
ELEMENT  FA  I LUBE 


1.  GENEiAL 

Figure  29  illustrate,  sere  tally  the  meaning  that  attaches  te 
same  of  the  terse  introduced  previously  and  required  In  the  subse¬ 
quent  development. 

Viet?  A  in  Figure  29  is  an  arbitrary  cross  section  through  a 
rolling  body  shewing  the  area  S  of  the  highly  stressed  voluae 
defined  to  be  the  volume  within  which  the  stress  critical  for 
fatigue  exceeds  a  designated  magnitude. 

The  highly  stressed  20ne  is  divided  into  8  large  nuaber  of 
cubical  cells  ©f  small  volume.  The  volume  is  so  selected  that 
within  it  can  be  contained  an  "engineering  slue"  crack  as  it  exists 
at  the  end  of  failure  Phase  I .  Each  c  f  these  cells  is  considered  to 
contain  a  defect  from  which  a  crack  aay  initiate. 

Defect  severity,  as  defined  in  Section  V,  is  a  function  of 
a  defect’s  size,  shape  and  compos i t i on , 

In  effect,  under  this  model,  the  ring  Is  assumed  to  be  built 
solely  of  cells  containing  defects,  each  having  different  severity. 
Naturally  for  clean  stools  a  large  proportion  of  the  cells  will  be 
occupied  by  defects  having  a  severity  close  to  unity,  (i.e.  with  no 
effect  on  the  strain  at  the  cell). 

View  B  of  Figure  20  is  a  side  view  of  the  rolling  body 
stretched  out  into  a  rectangle.  The  length  denoted  by  le  is  the 
circumference  at  the  neutral  axis  of  the  rolling  body  cross 
section.  View  B  also  shews  how  the  stressed  area  S  varies  with 
the  coordinate  f  for  a  general  loading  condition. 

Equation  (5.12)  of  Section  V  states  that  each  cell  in  the 
highly  stressed  zone  of  the  rolling  body  has  associated  with  it  a 
number  ef  stress  cycles  required  for  a  crack  originating  in  that 
cell  to  become  of  critical  size.  Equation  (5.17)  defines  an  addi¬ 
tional  life  increment  Njj  for  a  cell’s  crack  to  grow  to  critical 
size. 
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Tha  lives  ftj  in  any  get  ef  cells  are  statistically  independent. 

The  Phase  II  cell  lives  N|j  are  not  necessarily  independent 
ef  each  ether.  It  is  highly  likely  that  the  Njj  life  ef  eoe  cell 
aill  be  altered  if  a  crack  in  a  nearby  eel!  completes  its  Pfesse  I 
life  and  the  crack  extends  beyond  the  initial  eell  boundaries. 

Without  the  independence  assumption  the  probigs  of  finding  the 
probability  distribution  of  the  railing  body  life  is  quite  complex. 

The  treatment  which  follows  will  focus  on  Phase  I  life  under 
the  assumption  that  Phase  II  life  is  either  i)  negligible  compared 
to  Phase  I  life  or  ii)  constant  for  every  cell. 

The  Phase  I  life  predicted  by  equation  (S.l2)  is  a  determin¬ 
istic  function  of  the  defect  severity  and  the  operating  variables. 

Stochastic  considerations  enter  when  one  considers  that  the 
defects  in  a  steel  matrix  vary  in  their  severity. 

Every  defect  in  a  rolling  body  has  a  potential  Phase  I  life. 

The  rolling  body  life  is  the  smallest  ef  these.  The  defects  in  any 
one  rolling  body  constitute  a  sample  from  a  population  of  defeets 
with  varying  severities.  The  defect  producing  the  smallest  Phase  I 
life  will  thus  vary  randomly  between  rolling  bodies  so  that  the 
life  of  a  randomly  selected  rolling  body  is  a  random  variable. 

We  seek  the  probability  distribution  of  rolling  body  Phase  I 
life  ever  the  population  of  rolling  bodies  which  are  subject  to 
identical  material,  geometrical  and  environmental  conditions.  At. 
this  stage  in  theory,  w©  ignore  the  (relatively  small)  variability 
between  and  within  components,  of  the  matrix  properties  such  as 
yield  strength,  and  ascribe  the  variation  in  life  between  rolling 
bodies  te  the  variability  in  defect  severity, 

2 .  DEFECT  SEVERITY  DISTBXBUTION 

The  severity  of  the  defect,  situated  at  coordinates  x,  y  and 
z,ls  assumed  to  be  a  continuous  random  variable  independent 
of  the  defect  location  coordinates.  That  is,  a  defect  having  given 
severity  is  not  predisposed  to  be  located  near  any  particular 
set  of  x,  y  and  z  coordinates  (group  ef  vol::ae  elements  in.  the 
rolling  body). 
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i.  ’intauuwer' mumnnminuniuitn  i  m  itnrm  01  u  m  mTuiii  m  fo  qi  „ 


Let  F(d)  denote  the  cumulative  distribution  function  fee 
defect  severity,  the  probability  that  the  severity  d,  of  the 

defect  located  at  coordinates  sc.  y  and  z  is  less  than  a  value  tis  1  s 
given  by  the  continuous  function  F  ( d  ) . 


P  rob 


<  d|  (k, 


f(d) 


;  d>  i 


(6.1) 


=0  ;  d<  1 

Shere  ds  represents  a  realized  yalue  of  severity  and  where  the 
left  hand  side  in  the  above  statement  is  read:  "the  probability 
of  dz<d  given  g  get  of  values  k,  y,  and  z  is  . ”. 


The  fact  that  the  given  coordinates  x,  y  and  z  do  not  appear 
as  par  aae  ters  in  the  expression  for  F(  d )  states  what  was  said  above 
about  the  independence  of  the  severity  distribution  from  the  defect 
coordinates. 


Surface,  as  opposed  to  subsurface,  defects  pose  no  special 
problem  at  thig  point.  They  simply  represent  the  special  case  z-0, 
where  the  cells  on  the  surface  degenerate  into  square  platelets 
of  negligible  depth.  No  information  presently  exists 
concerning  the  functional  fora  for  F(d).  Two  forms  sill  be  advanced 
in  Appendix  VI  which  satisfy  soae  plausible  general  conditions  that 
any  defeet  severity  distribution  must  possess.  Other  distributions 
just  as  plausible  say  be  found  h©t?ever»  and  experimental  inforaation 
is  ultimately  required  to  choose  between  candidate  severity  distribu¬ 
tions.  Nonetheless,  valuable  insight  can  be  gained  regarding  the  sts“- 
tistics  of  possible  failure,  without  precise  knowledge  of  this  distri¬ 
bution,  The  diocuggion  will  therefore  proceed  in  terms  of  a  general, 
unspecified,  severity  distribution  F(d). 


3.  DISTRIBUTION  OF  "DEFECT  LIFE 


From  Equation  (5.12)  one  may  find  the  relationship  between 
a  defect.'g  severity  d  and  its  phase  1  life  Nj  as  follows: 


d 


- 1  .  B 

r  ^  N 


(6.2) 


where 

_t 

f  denotes 


the  inverse  ef  tho  function  T  (d)  defined  in 
is  a  non-randoa  quantity  depending  @n  A^, 


Section  V, 

{?.,  and  ?>„ 

y 


Equation  (6,2)  way  be  used  to  tranifsrs  tie  dtitribstien  of 
defect  severity  into  §  d  i  g  t  r  1  bn  t  i  on  sf  life  associated  with  th- 
de  f ec  t  pop u 1 stion. 

~l 

Assuming  that  the  function  T  lea  single  valued  and  aonotsni- 
cally  Increasing  function  sf  d,  one  say  writ©? 

G(Njl*.y.*)  =  ProbfN^Njj  =  Proto  f^J  (£>  3) 

=  1-Prob  |  *  (  f  )] 

Pros  equation  (6,1)  ^  *  -* 

G(Njj  x.y.z)  =  I-F[r  jfj]  ]  (6,4) 

Equation  (6.4)  is  the  distribution  of  life  at  coordinate 
position  x,  y  end  z  associated  with  the  population  of  defects 
having  the  severity  distribution  of  Equation  (6.1). 

Of  course  in  any  given  ro  1 1  i  ng' body, the  volume  element  at 
x,  y,  z  Hill  have  a  specific  defect  severity  and  hence  a  determined 
life  calculable  from  Equation  (5.12).  It  is  when  one  considers  the 
varying  values  of  the  severity  of  the  defect  that  can  occupy 
the  cell  located  at  coordinates  x,  y  and  z  in  different  rolling 
bodies  from  the  population  of  such  bodies,  that  a  distribution  of 
cell  life  results. 

Every  cell  has,  associated  with  It,  a  distribution  function  of 
the  fersi  of  Equation  (6.4),  therein  the  parameters  (since  they  de¬ 
pend  on  0,  eg,  and  S)  will  depend  upon  the  location  of  the  cell. 

The  life  of  a  complete  rolling  element  Is  identical  to  the 
life  of  that  cell  which,  of  all  the  cells,  has  the  shortest  Phase 
I  life  according  te  Equation  (5.12). 

Fros  probability  theory,  the  probability  that  at  least  one 
cell  produces  failure  before  Nj  cyeieg  is  the  complement  of  the 
probability  that  all  cells  survive  beyond  PJj  cycles. 

Assuming  independence  of  the  eell  lives,  (@§  one  say, 
for  Phase  I  life  considered  here,  based  on  the  choice  of  cell 
size  to  include  completely  a  croek  of  "engineering  size"),  the 
probability  that  all  cells  survive  is  the  product  of  the 
probability  that  each  survives.  Thus  H  (Nj)  is  given  by: 

H(Nj)  =  1  -ft  [l-G,  <V  ]  (6.5) 
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where  G{  denotes  the  distriiiution  of  Equation  (6 .4)  it  tie  i  - 1  ft 
eel!  in  the  structure.  The  §  y  safari  ■rri  denotes  that  the  t  eras 

evaluated  at  eaeh  cell  in  the  structure  including  tie  sarfic 
eel  Is  are  to  be  multiplied  together. 

Equation  (6.5)  is  of  great  generality  and  could  be  evaluated 
numerically  for  any  general  ehspe  and  leading  condition  if  the 
functiong  in  Equation  <h.  5)  and  the  severity  distribution  §  f  Equation 
(6.1)  its  re  k  seen  , 

These  functions  are  net  presently  known  however,  and  say 
not  be  practically  ds terminable  in  the  near  future.  One  has  as  a 
recourse,  the  option  of  asking  plsusible  assumptions  of  tie  sakaawn 
function  farms. 

In  view  of  the  uncertainty  with  which  any  function  can  be 
assumed,  it  is  preferable  in  this  regard  te  limit  the  number  and 
re s t r i c t i ve nes s  of  such  assumptions.  In  particular  it  preferred 
to  assume,  if  passible,  only  the  asyeptotie  behavior  of  a 
given  function  rather  than  its  specific  fore  since  this  merely 
limits  the  possible  functiong  te  an  adsissible  class. 

At  the  present  state  of  theory,  as  will  restrict  the  discussion 
te  applications  ©herein  the  distributions  Gj(ftj)  say  be  considered 
identical  for  all  values  of  i. 

This  corresponds  to  the  Lundberg-Palsgren  assuaptien  on  failure 
locations  for  the  case  of  thrust  leaded  bearings.  The  Lundberg- 
Palagren  assueptlen  is,  broadly,  that  all  points  inside  a  highly 
stressed  gone  ar@  equally  likely  t©  fail  and  no  point 
eutiide  this  gene  sill  fail.  For  a  thrust  loaded  bearing,  the 
highly  stressed  area  is  the  game  in  all  cress  sections,  and  thus 
the  above  stated  cage  applies. 

4.  ASYMPTOTIC  OISTBIBLiTIOH  0!*  SMALLEST  DEFECT  LIFE  IN  A  ROLLING 
BODY 

If  a  vgfidea  gaepia  of  size  m  is  drawn  from  a  distribution 
baring  a  distribution  function  F(x)  the  ordered  values  ,  x3 , 

.  .  .  Ke  (where  x%  ,  <  gs ,  <.  „  .  xB)  are  themselves  randea 
variables  hgving  s  di s t r ! but  1  on ,  over  repeated  samples,  of  a  fora 
which  depends  upon  order  nueber  i,  sample  sige  a  and  parent  dis¬ 
tribution  F  ( X  )  . 
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T.«  particular  £@ar  1  =  1B  that  is  for  tfe©  smallest  aesbe?  si 
tli©  saapie,  the  41  s&ribufci on  F](Sji  is  gives  bf6  C45>| 


Fs  ( s  ) 
l  i 


1- 


1-F(s  ) 


16. 6) 


As  the  sssple  sise  a  increases,  tbe  distribution  FjCSj) 
wader  ; ether  general  conditions  converges  to  §  specific  fora*  la 
tbe  ea««i  t?her>?  it  does  converge,  it  will  eoavarge  to  sae  of  three 
di  striiJutioasl  liras,  depeadiag  upc  the  behavior  of  the  pareat 
distribution  F(x )  in  the  vicinity  of  the  smallest  adsissibls 
value  fer  ss  (See  Gu&bel,  (4?)  p.  162). 


The  Weibull  distribution  is  one  such  Halting  fora  far  the 
di §  t r I batten  of  smallest  values  and  it  is  applleeble  te  parent 
di  i t f i b« t i o n s  which  have  finite  admissible  values.  In  order  for 
the  distribution  ef  X  j ,  t@  converge  te  the  W-elbull  di s t r 1  but  lea, 
the  parent  distribution  F(x)  erst  behave  like  a  power  function  in 
in  the  vicinity  of  the  ainiaua  admissible  value  .  that  i&i 

,  k 


F(x ) 
x— 


0(x~xar 


(6.7) 


Where  $  and  k  are  positive  constants  ana  x0  Is  the  siaallest 
admissible  value.  If  Equation  (6.7)  is  satisfied,  it  has  been  shown 
(Epstein  (2))  that  for  x<xg,  Fj  (x.)  has  the  following  fora; 


Fl(x1t  -  1-exp  [ -03(xl-xo)k 


(6.0) 


Thug  the  distribution  ef  rolling  body  lives  (shortest  defect 
lives)  sill  be  0 ei bull  distributed  if  the  patent  distribution  of 
eell  Phase  I  life  CS(Nj)  satisfies  the  condition  of  Equation  (6„7) 
wherein  the  ainisus  Phase  I  life  corresponding  to  infinite 
defect  severity  i  s  so  o. 

Vis.  S>]]=  8  B>'‘  «*■« 

If  Equation  (6.9)  Is  satisfied  for  the  actually  applicable  functions 
F  and  r,  the  di stri but ion  H  (Nj)  becoaes,  by  Equation  (6.8)  s 

Ni  '* 


H(f¥ 


)  =  1-exp  ~Bp(Nj)k  -  1-expjJ- 


where  N*  =  l/(m&) 


/k 


1  N® 


(6.10) 


Saepie  §i§e  s?  is  the  total  nusber  of  defects  (cells  is  a  rolling 
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Beaeting  the  number  ©f  defect?  per  ©nit  velarne  by  T),  ©ne  l§s; 


M  =  TlSlg  =  1|V 


Sabi t i t  a  t i ng 


(6.11)  in 
H(N  ) 


(6.10)  gives 


Equetien  (6.12)  §h@wa  the  stressed  velum®  appearing  a 
scale  parameter  ©f  the  Seibull  distribution.  This  is  the 
effect  of  the  Lundberg-Palmgren  formula,  which  is  believed 
in  geed  agreement  with  existing  ball  bearing  fatigue  data. 


Me  have  new  ehawn  that  the  Lundberg-Palmgren  formula 
compatible  with  the  fatigue  failure  model  developed  ab©ve 
the  restrictions  pn  the  tw®  functions  T  and  F  eabsdied  in 
(6.9)  abeve. 


(6.13) 


(6.12) 

s  % 

val use 
to  be 

is 

under 

Equ  ex  ion 


SECTION  VII 


IACSOSCOPIC  ST1ESS  CALCULATION  FOB  THE  HIGHLI  5THE-5SE 0  VOLOBSS 

IN  A  HEiTZI AN  CONTACT 


la  rolliag  contact,  there  euists  a  '‘reversing  Bhesr  stress” 

(1)  which  has  iosg  been  regarded  as  the  best  criterion  of  rolliag 
contact  fatigue  strength;  ueuiliy  on  the  ground  that  its  range  of 
variation  is  always  equal  to  or  greater  than  that  of  the  maximum 
ghear  stress  and  that  there  is  some  evidence  suggesting  that  total 
amplitude  of  shear  stress  correlates  with  life  (1,29,  43),  Tkis 
assamptioB  was  used  by  Lundberg  and  Palmgrar.  (1)  who  stated  that 
most  subsurface  fatigue  cracks  are  generated  at  a  depth  equal  to 
that  of  the  reversing  shear  Cat  the  central  cireuaferentisl  plana). 
Latest  Greenert  {28)a  ia  testing  toroids  of  various  curvatures, 
found  that  the  reversing  shear  range  gives  good  correlation  with 
fatigae  life.  In  available  literature,  the  shear  stress  has 
been  determined  only  for  the  central  ci rcuraferenttai  plane  of  a 
rolliag  element  oa  which  the  maximum  shear  range  is  equal  to  that 
of  orthogonal  shear  stress  ( ? yz ).  Since  the  fatigue  life  i§ 
closely  relates  to  the  highly  stressed  volume,  it  is  necessary  to 
know  the  volume  of  material  subject  to  a  certain  level  of  shear 
rang©  for  p  given  applied  Hertzian  load.  The  determination  of 
this  volume  requires  a  knowledge  of  the  stress  distribution 
under  a  Hertzian  load. 

A  computation  of  (reversing)  shear  stress  range  in  a 
Hertzian  elliptical  stress  field  has  been  performed.  Previous 
workers  have  studied  the  maximum  alternating  shear  stregs 
on  the  caster  piano  of  a  contact  x  -  0,  (see  Figure  3).  At 

x  =  0,  the  shear  stress  is  maximum  in  a  plane  at  a  keown  distance 
belew  the  surface  and  parallel  to  the  XY  plane.  As  ■  shown  below,, 
the  alternating  shear  stress  for  other  planes,  (a^O)  is  Hasiauro 
on  planes  inclined  to  the  XY -  plane,  and  therefor©  is  not  ortho¬ 
gonal  „ 


The  computation  presented  here  ig  a  generaligatioa  @i 
Laadberg's  solution  which  was  confined  to  a  central  ci reuaf ere  a- 
tial  plan©  only.  The  present  solution  provides  a  means  to  deter¬ 
mine  the  stressed  volume  at  various  stress  levels,  eoasidering 
the  shear  range  as  the  critical  stress. 

The  baaie  aathsraati c al  formulation  for  the  subsurface 
stresses  at  all  locations  in  a  Hertzian  stress  field  is  presented 
In  Appssiix  I,  These  formulas  were  developed  by  Lundberg  and 
Sjevall  ia  19§l  (34). 


T  h  ©  eurreat  computation  involves  the  rie termi net io  B  B  for  a  ay 
point  Ip  e  ero§s  aectioa  (perpendicular  to  ro 1 J i ag  ri i rec t i ea\  of 
thereage  ofsigiBue  shear  stress  (  t  acting  is  the  direction 
of  rolling  at  &  certain  value  of  y  (where  y  is  the  distance  from 
the  central  plane  of  symmetry  is  the  direction  of  rolling)  but 
on  any  plane  inclined  to  the  surface  (with  an  angle  at  ) .  The 
determination  of  the  values  of  y  and  e  st  which  tg  is  maximum, 
requires  an  iterative maxirai gat i on  process. 


The  following  stress  analysis  was  conducted: 

a)  Based  on  the  equilibrium  of  stresses  acting  on  a  sub¬ 
surface  element  shown  in  Figure  3e  the  shear  stress 
acting  on  an  arbitrarily  inclined  plane  can  be  expressed 
as  the  combination  of  two  orthogonal  sheer  stresses 

Tzy  antl  Txy  ^  ° e ' 


T0  =  Tzy  0030  +Txy  sin  B 


(7.1) 


where  T2Jr  and  TKy  can  be  computed  from  formulae  available  in 
(34)  for  gives  values  of  x,  y,  z  and  a/b  (the  ratio  of  major  axis 
and  minor  axis  of  a  contact  ellipse). 

Since  both  TKy  and  TZy  are  odd  functions  of  y9  "Tq  is  also 

an  odd  function  of  y.  Doth  vxy  and  tz„  vanish  when  y  s  0. 

Furthermore,  glace  ?8y  =  0  for  all  y  when  %  =0S  the  amplitude 
of  i0ia  equal  to  that  of  t  at  the  central  plane  in  the  direction 

of  rolling  (where  x  =  0).  This  justifies  Lundberggs  analysis 

(1)  of  setting  the  amplitude  of  Tzy  equal  to  the  muxisuaj  shear 
amplitude  for  x  =  0. 

b)  Fop  given  values  of  x  and  z  the  maximum  value  of 
eaa  be  found  by  setting: 


Ha  =  T 

ae 


xy 


cos  0  -  Tyg.sin  6  =  0 


(7.2) 


and 


Ha  =  d 


( Tgy)  COS0+  i  ( Vxy).S 


in  6  =  0 


9  y  d  y 

These  two  differential  equations  can  be  reduced  to: 


(7.3) 


“  -  (  + 


d  y 


xy 


zy 


)  =  o 


75 


,11  mm,  null  . . . I 


which  isapllgg  that  when  T  |y  +■  i  ly  is  raaxinum,  7  g  is  also  a 
msxlsuB.  Therefor© 9  a  value  of  y|  ,  corres ponding  to  maximum 
Tiy  *  ^ly  eSH  b©  found  fro®  the  stress  calculation  and  using 
the  tallies  of  and  ?sy  at  y  =  y^t  the  value  0  =  €h  for  which 

is  BisisuB!s  eaa  be  derived  fro®  Equation  (7-2) 


l.S8; 


tan  8j  =  <  'rny/^y)y=7l 


(7.4) 


Tbs  rnguisua  shear  range  7 g  is  equal  to; 


2(’» 


) 

max . 


2  (  t@) 


y=y  i 


Using  the  values  of  ®i  aad  thus  obtained,  the  maximum 
value  of  7q  eaa  fa©  obtaiaed  from  Equation  (7.1)  after  substi- 
tut ioB. 

Using  &  digital  computers  numerical  computation  has  been 
made  for  values  of  vgy  aad  tq  (in  terms  of  the  maximum 

contaet  pressure  p0  )  as  functions  of  y/a  and  ®. 

Figure  4  plotse  from  these  computations  g  the  variation  of 
Tg  as  a  fuaefcioH  of  k  aed  %  Sot  a  contaet  with  a/fa  =  10. 

The  variation  of  9j  as  a  function  of  %  and  z  is  shows  by  dotted 
curves  i  a  Figure  4# 

At  k  =  a»  all  of  thetg  aad  91  curvos  very  nearly  inter- 
aoet-  This  iedicotos  that  at  x/a  ~  1.0  the  iaasiauia  stress 
range  is  equal  to  0.04pe  independent  of  depth  s  for  the  z  values 
shows  in  Figure  4e  Further,  the  plan©  os  which  the  maxima®  occurs 
is  defined  bj  ©j  =  35°  independent  of  z. 


Figure  5  shows  the  variation  of  v  and  7  gjith  y  at  a 

fix©d  laeatioB  in  the  ring  cross  sectiona  i.e.  II  x  =  8b  aad 
a  -  0.6b.  Figure  ^ shows  that  for  a  given  depth  g8  the  shear 
stress  range  fe-tneheg  a  maxisus  at  s  =  0  aad  decreases  nsono- 
toaieally  with  g. 

Using  the  sstlod  described  above,  th@  contours  of  equal  shear 
rgsige  on  a  ring  c  ros  s=see6i  sa  for  e  gives  major  axis/  minor  axis 
ratio  (characterising  the  contaet  ellipse)  can  be  drawn,  based  oh  the 
values  of  shear  range  computed  for  a  large  number  of  grid  points 
equally  spaced  is  the  highly  stressed  sea©  of  a  ring  eross  section. 

Figure  6  shows  contours  of  equal  shear  range  for  the  @eo©n~ 
ferisity  ratio  e/fa-10  representing  a  typicgl  ball-race  contact  of 
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i  ©  deep  grease  ball  basting.  These  cot  Lours  ur&  drawn  based 

i  0t?  0  0$  grid  point  ejsies  in  a  quadrant  o£  the  plane 

;  containing  the  ssajor  axis  of  the  contact  ellipse.  It  is  noticed 

■  that,  the  shear  range  along  i\  —  0  i  s  very  close  to  tbit  in  §  ese= 

!  tact  of  two  lifiBitgiy  long  cylinders, 

'  ^  limiting  ease  for  s/b  -  »  has  been  de<  uced  from  formulas  in 

£34)  which  corresponds  to  an  infinitely  narrow  cosiest  ellipie 
:  (see  Appeedix  II,)  la  this  cess,  the  stress  distribution  ia  the 

.  central  ei  rcusafereatisi  plage  (x  =  o)  is  idealist!  to  the  plaae 

j  a  t  v  o  i  h  solution  far  a  Heraiaa  two-dime  a  s  1  c  n  a  ?  coat  act.  Using  the 

i  results  ©f  the  lieitiag  easee  it  ia  possible  ta  plat  tfee  vals©§ 

;  shegr  rings  ea  a  pl§a@  eiife  eeordiaates  x/a  and  g/b, 

j  Figure  T  shows  tfee  coafcours  of  equal  sheer  raag®  for  s/b  s® 

j  i.e,  aa  infinitely  narrow  coatact  ellipse  with  its  aajor  axis 

i  lying  on  the  X~axlss  perpesdicular  to  tfee  direction  ©f  rolling, 

\  By  comparison  silk  figure  6,  it  is  seen  that  the  ceaiears  of 

;  equal  shear  rasge  for  tfee  tee  cases  are  nearly  identical  except 

■  tfegt  |a  fcfe©  eis@  a/b  = © #  tfee  level  of  Tg  lies  slightly  deeper 

|  uade?  ths  surface  tfeaa  for  a/b  =  10. 

i 

j  Colours  of  equal  shear  range  ?g  for  th©  ease  a/b  =  7.5 

•  have  ©1st  best?  obtained  and  these  are  shoes  ia  Figure 

1 

|  Figure  9  plots  tfee  area  S  expressed  as-  multiples  of  ab  aad  ea- 

|  closed  by  contours  of  equal  Tg  as  g  function  of  Tg/PBas.  for  values 

i.  ®f  %/b  ~  7c5„  10,  and  ®  o  Tfe©  pcleis  whore  these  curves  intersect 

r  tbe  abscissa  correspond  to  the  eaxisus  values  of  v  tferoaghost 

i  stressed  region,  i.e.  2  ?0„  It  is  noted  that  the  curves  lie 

’  very  close  to  each  ether.  I&e  thro©  carves  coincide  even  a©re 

closely  if  one  plots  S/as0  vs  rB/2rQ,  by  auliiplylsg  the  vertical 
coordinate  of  Figure  9  by  b/zo  and  the  horizontal  coordinate  by 
I  t©>  nMss@  3 g  i §  the  depth  of  tfee  point  below  th©  surface 

|  ebich  is  subject  to  the  maximum  shear  stress  raag<  2v0.  The 

values  of  b/ze  aad  Pga„/To  8g  functions  of  a/b  ore  obtained 
frea  (34)  aad  are  tabulated  ia  fable  3.  Figure  10  is  a  plot  of 
S/aZg  vs  7g/2ve  and  shows  that  the  curves  having  parameters 
s/b  -  7,5,  10  and®  nearly  coincide.  This  means  that  tie  stressed 
§ree  S  at  any  given  value  is  proportional  to  the  product  of 
j  0  *o  wEiere  ®  is  fb©  sesi-eajor  asis  and  2  is  the  depth 

j  leoetioe  coordinate  of  This  approximate  relationship  appears 

j  to  be  v@.lid  for  7,5  4  a/b<  «  ;  a  rsage  eaccepossing  tfee  usual 

!  diaofsgioes  of  a  ball-raceway  coataet  ellipse,  including  line 

j  eooiaet.  This  fiadiag  supports  La  ndber  g=*Pa  lagree  s  §  analysis 

i'l)  ia  wfeiefe  tfee  stressed  area  in  tfe©  ring  cress-section  is  sg« 

■  sueed  t,o  be  proportional  to  a?g. 
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A  selection  of  tit©  values  to  he  aged  in  defining  tie 
stressed  area  can  be  eade  by  setting  Tg  equal  to  a  threshold 
stress  level  for  plastic  deformation  for  th#  gi?@a  aatsrial 
composition  and  hardness,,  If  ^elected  in  this  b^risse,  Tg 
is  a  constant  for  a  given  material.  For  future  use8  it  is 
more  convenient  to  plot  S/az  as  @  function  of  2tg/tg  shewn 
os  dotted  line  in  Figure  10= 

This  eur?9  shews  that  plasiicslly  stressed  (aacroscopie) 
area  exists  only  for  2t0/tg>l-  The  size  of  the  plastically 
stressed  area  eio  be  appr oxiraated  by  i 


B.  2, 


2  T„ 


Ta 


nO.75 


T8 


(7.5) 


This  formula  shows  that  the  @rea  enclosed  by  th©  eoatoar 
of  equal  '?g  can  be  expressed  approximatsly  as  a  fuaotioa  of 
Tg.  It  is  also  of  interest  to  know  the  average  shear  range, 
*TJ^ev@  »  ia  stressed  area,  S,  enclosed  by  the  contour  of 

equal  Tg6  Appendix  I?  gives  ife©  details  of  the  computation 
using  Equation  C7.5).  She  result  yields  th©  following  spprcxi^ 
«8tion; 


S 


830 


-6.45 


4-  12.9 


(7.6) 


where  tfR^ave  s  average  segaitud©  of  a  laaxiisuB  reversing  shear 

rang©  in  a  closed  contour  of  equal 

3  ~  area  enclosed  by  the  closed  eoatour  of  equal  Tg 

a  =  s©  lajor  axis  of  ths  contact  ellipse 

^  =  maximum  reversing  shear  stress 

z0  -  depth  where  ?0  exists 
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BETEBM1NATION  OF  SHEAS  STRESS  NEAR  ASPERITIES 


To  demonstrate  s  possible  mechanism  for  the  generation  of 
plaafic  dsforssiion  In  asperities,  an  elastic  analysis  is  present¬ 
ed  below  for  the  stress  distribution  in  a  contacting  gcpsrity 
using  an  idealized  asperity  profile.  The  result  of  this  analysis 
may  enable  one  to  predict  the  location  and  severity  of  a  plastic 
occurrence  and  take  into  account  the  relevant  variables  such  as 
lubricant  film  thickness  and  surface  roughness.  It  should  also  he 
noted  that  in  the  foil  owing  analysis  the  friction  at  the 
contact  surface  is  neglected. 

Figure  11  shows  the  distribution  of  slopes  for  threq  V. 
typical,,  abrasively  finished  surfaces.  From  those  distribu¬ 
tions,  it  is  seen  that  the  typical  slope  of  abrasively  finished 
hard  steel  surfaces  varies  widely.  For  8  ground  surface,  the 
95th  percentile  slope  is  ©95  ~  29°;  for  a  mush  smoother  honed 
surface,  the  95th  percentile  is  ©95  =  3.8°,  while  for  the  still 
smoother,  lapped  surface,  it  is  only  Ogg  -  0.79.  The  80S 
surface  roughness  values  of  these  three  surface  finishes  were 
found  to  be  o  =  13e  1.0  and  0.4  pin,  respectively. 

Quantitative  insight  into  the  plastic  occurrences  at 
asperities,  as  a  function  of  parameters  listed  above,  can  be 
achieved  by  introducing  a  simple  mathematical  model  in  which 
the  contact  of  0  single  asperity  having  an  idealized  profile 
with  an  elastic  half  plane  is  considered.  This  model  contains 
no  assumption  regarding  the  height  distribution  of  a  population 
of  asperities. 

Figure  12  shows  an  idealized  two  dimensional  asperity  shape, 
similar  to  what  one  would  expect  to  find  on  ground,  honed,  or 
lapped  surfaces,  formed  by  0  multitude  of  cuts  by  sharp  and 
straight-edged  abrasive  grains.  The  asperity  is  a  plane-sided 
ridge  with  o  curved  tip  of  radius  g»  The  slope  angle  of  the 
sides  is  6  which  varies  with  the  process  of  surface  finishing 
and  is  accessible  to  experimental  determination  (6). 

The  plane  contact  problem  with  the  above  described  profile 
and  a  straight-edged  half  plane  can  be  solved  based  on 
Muskheli shvili fs  method  of  singular  integral  equations  (36), 
provided  that  there  exist  no  sharp  curves  or  corners  along  the 
entire  profile.  The  derivation  of  the  solution  for  this  eon° 
tact  problem  is  given  in  detail  in  Appendix  III. 
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h  iaa©r3eil  eeeputatien  af  surface  preaears  seas  per  farmed 
gad  is  glows  ia  Figure  1 2»  bssgi  on  tfes  eese  a  -  2b  wherein  the 
width  of  tfes  contact  region  is  ©qual  6s  iwic©  the  width  of  the 
earned  base.  Computatios  for  ether  v§iue§  of  a/b  is  equally 
£a@§lbl@e 


The  surface  pressure  distribution  plotted  in  Figure  12 
wea  coaputed  based  oa  §  closed  form  solution  derived  £roa  ( 36 ) e 
H o-^e?er 0  for  £fe©  sub-gurfooe  stress  distributioa  there  is  no 
closed  fora  solution  available  for  eoaputatioa  sed  i  aueserieal 
teefesiqu©  is  required  to  obtain  tfes  sub-surfact  stresses^ 
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HALF  WIDTH  OF  CQtiTAGT 


Figure  13  shows  the  eoateayg  ef  equal  oetii@dral  shear 
§tr©§§«  ia  Issfia  o I  ifes  product  el  ¥ousgss  Modulus  aad  sspsrity 
slept,  §’§  i a  tie  crass  seetioa  at  tie  plant  of  tae  eoetaet  d©~ 
pie  ted  ia  Figure  12,  It  eis  be  §@@a  that  II®  d©ptfe  wfe@re  the 
misieua  value  of  octahedral  si§ir  stress  ©sours  is  shout  O.S3 
ef  tbs  half=»idtb  of  the  asperity  contact  regies  00  the  surfec®, 
For  a  typical  asperity  eoatget  width  oa  a  fiaely  finished  surface 
(1  p  in,  BMS)  (6)^  this  depth  i?  of  the  order  of  teas  to 
huadrsds  of  aieroinehes.  Thus,  high  shear  stresses  are  iadeed 
gsaersted  close  to  the  surface, 

fis©  following  relationship  has  been  derived  for  the  eagsi-> 
tad©  ef  the  slssr  stress  along  the  asig  of  a  symmetry s (Appendix 
III  ) 


45° 

where 


V 

u 


""Z-9'!"  1  •[  v  "  taB  ( Y  •  ? sn  v) 
rr  b  *  t  2 

=  *  /a 

=  %/  (1  +  fa)  1/2 


y (rr/2  -u)j  (8„l) 


=  cos“* (b/a) 


Figure  14  shows  a  plot  of  the  saxisus  value  of 


'45° 


computed  from  Equation  (6.1)  against  b/a.  The  maximum  value 
of  ?4gdg  seen  to  be  proportional  to  ©  and  is  a  function  of  the 
ratio  a/b„  lies  b/a  •“❖(),  i.e.  when  the  contact  area  is  wills 
is  comparison  with  th©  asperity  tip  width,  the  following 
limiting  case  exists; 


W- ’E*e  /^s 

45° 


SS8JJ 


0,9O4'/TT 


(0.2) 

(  8.  3) 


*Tha  octahedral  aha§r  stress  is  defined  as  tocj  =  *\f(  ol  =  ea  )  s* 
“  Os T®  +  ( Og  -  Cj  )s  /3,  where  ffj ,  c8 ,  o8  ar©  the  thro© 
principal  stresses,  whereas  the  Yen  Mises  yield  stress  is 
OB  -  3/1  fget/2  Note  that  Teej  acts  on  one  of  the  faces  of  a 
recular  octahedron  with  vertices  en  the  axes  X,  x  and  Z„ 
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and  tie  eoataet  fores  I g i 

p  — e  '  e  §/  tt 

Tiis  se  relationship!  isply  tie  following  £ Btaresti eg 
facts  ^les  b/a  — £  08  i,#e  tie  ©leetie  deforaitioa  of  aa  asperity 
is  large  in  eospgrisoE  with  its  tip  width  (for  §xasple9  is  tfe© 
cage  of  a  sharply  tipped  asperity  or  heavy  siastic  depression 
of  ae  asperity),  This  is  the  case  for  low  k/a  (tie  eiaisua 
EH$  film  thieknags/eesposit®  surface  roughaess  115  ratio) 
v  a  1  a  ©  s  s 

1.  The  saxiau®  shear  stress, ?a8X  ,  approaches  a  eoagtaat 
value  proportional  to  the  asperity  slop©  sEgl©B  bat  itsdepeadeat 
of  the  load, 

2.  Tie  depth,  feas  s  of  the  point  at  which  %gg  occurs 
increases  with  the  l/2~tfe  power  of  th©  gesi-=@ldtl  s®9  of  the 
asperity'  coniicte 

3.  Tie  load  carried  by  tie  individual  asperity  contact  is 
proportional  to  tie  contact  width  sara  and  the  eagle  ©.  This 
iiragl©  asperity  model  yields,  asymptotically,  a  proportionality 
between  P  and  a,  in  agreement  with  the  &rcharde§  postulate  (35), 
elicit  states  that  in  dry  contact  the  real  asperity  coatacf  area 
iacreasei  linearly  with  th©  load. 

Of  course,  there  must  be  cases  in  which  the  degree  of 
indentation  or  a/b  is  not  large.  This  occurs,  geaer^lly,  when 
g)  these  is  a  thick  lubricant  film  separating  th©  two  rough 
surfaces,  i8e,  t/ais  large  or,  specifically,  when 
b)  tie  tip  radius  of  th©  asperity  on  real  surfaces  is  large 
in  comparison  eith  the  asperity  spacing  30 

The  following  presents  a  simple  sod©l  for  asperity  inters 
action  for  the  purpose  of  relating  nesr^gur f ac©  plastic  occur* 
rsnaes  to  the  ainiau®  BOB  file  thickness  and  surface  roughness 
partasterg. 

Figure  15  (&)  slows  the  two  diaeasictaal  eontaet  of  a  siagl© 
asperity  assumed  to  be  rigid,  with  an  elastic  straight  edged  half 
plane?  th©  degree  of  approach  of  th©  two  bodies  is  controlled 
by  the  EQ3  minimus  fils  thickness  h.  For  cos?enieac@,  th© 
verticil  distaaee  b©twe@i'  the  outer  line  of  tie  asperity  profile 
end  the  asperity  tip  is  set  equal  to  3ewh@r@  q  is  th©  IMS  valu© 
of  the  asperity  profile  and  is  of  tn  acceptable  order  of  Hagai- 
feud©  approximation  for  real  profiles.  For  th©  hypothetical 
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pr@£ilee  felts  esisaptiea  |g  arbitrary  feeeaas©  II©  fslitfi  @f 
the  asperity  nr©  safe  deflagd,. 

The  vertical  distane©  beineea  16©  asperity  tip  and  the 
eeataet  ©dg©  §  is  ©blal.s©4  f?ss  tie  geesetry  ef  Figure  iS9 
B&iel  is 


8  =  6e  +  ''a-’b)  tin  9 

where  =  S  (1-c  s®),  the  depth  of  the  ssperitj  tip  aa«S 

b  —  M  g i »t  9  C  B  ©  5  ^ 

Oat  stay  consider  §  eisJUitad©  of  rough  surfaces  with 
asperities  of  the  shape  shows  is  Figure  IS  with  different 
values  of  or  .  A  reasoaabl©  asguaptiea  for  the  degree  @£ 

Shair  rounding  at  tie  tip  is  tint  §0-  \a  ghee©  Kiss  ©©uaiaut, 
proportional  to  1,  i3e,  that  fefte  rounding  occupies  e  ©©aslant 
iraetisa  of  tfeo  esptrity  height  between  tip  aed  prsfil©  e©at©r~ 
li  a©» 

Frea  iqaatieg  C9.5)  oae  obtains,  using  60  =  A  a 


g  =  X.<?-(l-ees  S )“l  (a. 6) 

fr  =  X-c-sla  ©’(I-ces  9>°*  C80D 

The  deforced  profile  ©f  the  elastic  half  plana  is  eeataet 
will  tha  rigid  asperity  is  gi?ea  by  the  following  forguSa  C4) 


v(s)  =  1  ^"^ss  -  a3 

.  Pa  f  ,.  1  ^ 

"a  ( t-s)^  a®-!"8 

iffl  m 

-  i  f  siua 

dg  dt 

n  -a 

4  t-s 

sh@r©  £'(t>  =  t/i 

!tl  <  b 

=  b/i,  b 

<  1 1 1  <  a 

0i”  ?Ck)/s  =  @/B/(x/a 

,  b/a) 

«i=9) 

The  ab§xi=,£3£’SHj.§  ean  be  integrated  after  espasdiBg  the 
iatagraad  \fa*~x3/{  t-x)  into  @n  ascending  p©w©r  series  in  a/jc.  Th© 
©eepat©d  dleensi oaless  deforaalioa  outside  the  soataet  regies 
is  plotted  i@@is8t  x/e  ia  Figure  16  using  the  ratio  b/g  @§ 
th©  p§ras©t©r8 
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la  e§§©  tf  a  real  agperiiy  eosiact  t-ro  stetl 

ssrfeceg,  th©  ©gfstrlty  is  ist  rigid.  It  lag  the  fta©  ©Isatic 
prs®©rti§s  ai  tla  half  plus©,  la  addition,,  both  surfaces  fee?© 
asperities*  Per  (39),  t  seateaiaat  approximation  ef  the  ees<= 
tact  possesses  is  arrived  at  if  considering  ©ae  surface  hfriag 
a  reughuata  prefils  ©quel  te  the  composite  of  tit  two  setae! 
serfage  r  eugh'ne  s  se  s  ,  while  the  ether  surface  is  flat,  C37)  sh§w§ 
feos  tb#  composite  roughness  sag  be  calculated  for  egdeferised 
asperities  (among  other  formulas,  one  finds  aa  =  ux s  +  gBa 
@§  quoted  before).  Assuming  that  the  same  composite  feraisg 
procedure  is  acceptable  frr  deformities  e al cal st i oas9  Figure 
15  (A)  reprgseats  aa  asps  ity  of  the  composite  reugfegifi.  la 
order  to  esleulase  deflection,  this  composite  asperity  mast 
be  eeasidered  elastic  and  the  smooth  half  plane  rigid. 

Assuming  lew  slept!  §si  large  radii  for  the  asperity,  the 
deflections  shews  it  Figaro  IS  for  the  elastic  hsli-plsa© 
will  still  apply  bat  mass  sow  fee  measured  separately  from 
its©  aadefesaed  asperity  profils.  In  order  to  be  able  to  do 
so,  horigental  reference  lines  must  be  ©itablisiad  oh  1. 1»©  de= 
served  half  plasie  (tie  satire  surface  of  the  infinite  half 
plaea  in  castect  is,  strictly  sp©akiage  deforest  ahd  the  ©agii= 
f.sie  ef  ike  deflection  does  sot  converge).  A  reasaaebl©  refer¬ 
ees©  liue  is  the  heriseatal  lino  drawn  through  two  symmetrical 
paints  ©a  the  profile  at  a  iisiaaee  of  B/a  from  the  line  of 
sjss?aetry  of  the  asperity,  where  Sis  the  asperity  spacing. 
Scferwatioia  ef  the  surface  beyond  these  points  is  thoa  to  b© 
seglected,  Fer  a  symmetrical  asperity  profile,  tha  profile 
peinti  at  n  distaae©  of  S/s  free  th©  asperity  tip  lit  ©a  th© 
profile  eeatcrliae.  fhs  use  ef  those  points  as  a  iefarastlea 
rofereae©  implies  that  all  asperity  deflection  occurs  is  the 
pertiea  area  the  ceaterliaa  t©  the  tip. 

Hith  these  rtfertiee  points  established,  one  obtiias  the 
iefesaed  shape  ©f  the  asperity  by  subiruetiag  the  shaded  area 
frea  its  aadefarsed  shape.  <§ee  Figure  15  (IK) 

fh©  average  E10  film  thickness,  h,  betwtes  two  rough 
sarf.isei  Its  beta  iasiied  <375  @s  tbo  distaac©  between  their 
prefil#  of sterli ses.  If  the  composite  roughness  profile  is 
used,  tisea  w  is  the  diseases  between  its  centerline  sad  the 
sa©ftl  r@£©r©ae©  surface.  By  this  definition  on©  obtains 
free  Fife?©  1 0 : 

h+6  ^v(S/2)  -3a  (9.10) 

where  3e  is  the  undefarstd  §§pt?i  ty  height  from  the  profile 
eealerliae  te  the  tip. 


Is  diaeaglenleaB  for®,,  this  tquation  ogi  hs  written  b f' 
seats  a  oi  Iqueiieni  <  8 » 4  >  and  (0.73  as 

4=3  -b/a  -  v(S/2 )/g 

t) 

=  3  -  60/ a  -  (fi-fo)tan  8 /''j  -  v  C  §/2 )  /s 

••=  3  =  ?.  -  K  (  a/b=i  )  •  <  1  +  cogO)  /cose  -  v(S/2 )/a  (§,H) 

The  knows  function  v(x)in  Equation  (8.9)  can  be  re&rraig- 

eri  in  the  foilcsiu^  fora; 


By  Bttsns  of  Equation  (8=5)  sad  (8.6)  one  has 

a 

*(S/2>/a  =  \<1+C8.9)  •  ({)  /(|b,i,i)  «8-13) 

For  given  values  of  X.  9  and  a  ,  one  can  date  rain-'  h/o 
ue  u  function  of  a/b  (or  vice-versa)  by  using  Squat,  ma 
(8.11)  and  (8.13). 

Computation  oih/e  has  a’so  been  performed  for  the  case 
that  a/B  <1,  l.s.  the  region  of  contact  of  the  asperity  falls 
within  the  curved  tip  without  touching  the  straight  sides  of 
the  asperity  profile.  This  caB©  occurs  only  when  h/o  is 
relatively  large  but  less  than  3.  (It  is  noted  that  in  the 
overeat  sod©!,  there  is  no  asperity  contact  for  h/<y>3). 

Figures  IT,  18  anti  19  plot,  for  the  ground,,  hoaed  aad 
lapsed  surfer ‘t  finishes  discussed  above,  a/b  as  a  function  of 
h/.v  .  Sbi  iKraaeffE  for  eseh  curve  is  the  relative  tip 
width,  sjlieh  is  a  function  of  the  tip  radius  E. 

The  curves  show  the  general  trend  that  b/a  decreases 
Kith  sg  h/V  .Since  sbe  is  a  constant,  deeressiag 

b/a  n,’v8s;isj  the  contact  width  ea8.  Since  the  aaxi- 

(gas  ake,.i  Ptrpsi  increases  with  increasing  9&e  (or  ieereaeiBg 
b/a),  o ae  else  rs>  e  test  for  a  given  asperity  profile  tis© 
eigj;irnu@  sh-.t  ssr-jfg  ia  the  asperity  increases  eith  decreasing 


It  is  also  9 boss  is  Figures  I?  -  19  that  for  a  gifts  b/o 
v@lugB  decrtiiisg  value®  of  X  ar  1  (sharper  tipped  asperities) 
esrrespoid  ts  iieressisg  e/b  or  the  eeataet  width  ’ a 1  ssS 
ehaar  stress  faaH  .  Mows?©?,  it  is  recalled  that  for  large 
vsiaes  of  ’a*  (  or  b/s=# 0  ) # Taas  sill  reach  as  asymptotic 
value  isdepeadest  of  I,  The  relationship  between  b/a 
and  ?@a5,  for  the  three  surface  finishes  ei aeasssd  Mere  is 
plotted  ia  Figure  20,  which  permits  tbs  scaling  of  the  ordi= 
nt@s  of  Figures  17-19  la  terms  of  tbix,  ia  addition  to  b/a. 

Ibis  feei  been  done  by  showing  a  eoalinssr  ordisgi©  scale  on 
these  plots, 

From  the  above  rsiulia  based  os  a  simple  plane  asperity 
interaction  model,  it  is  possible  to  relate  the  eaxisua  near 
surface  shear  stress  to  three  paraseters  explicitly  gives  for 
a  gives  lubrication  condition  end  surface  finish.  Basely  &/©.  8 
and  E/a  where  h  =  sinisue  fils  thickness,  g  »  combined  surface 
roughness  (rss),  9  s  asperity  slope  angle  and  I  s  asperity  tip 
x adi ui. 

It  can  be  seen  froa  Figure  17  that  for  ground  surfaces 
with  a  typical  asperity  slope  angle  9  =  29®e  th®  eaaieus  shear 
stress  level  is  quite  high,  i.e.  at  h/g  =2  and  k^Q.5e  the 
maximal  shear  stress  is  ~  106  pel  which  is  considerably  higher 
than  the  magnitude  of  maximum  Hertsisn  shear  stress  usual  In 
rolling  contacts.  It  is  expected  that  plastic  floe  sill  occur 
causing  a  ^blasting”  of  the  tsp&rliy  tip.  For  a  smoother  sar= 
foe©  finish*  ®.g.f  e  honed  surface  with  6  s3.8°  and  a  1.8  pin., 
it  can  be  seen  from  Figure  1®  that  the  maximum  shear  ifregi  is 
considerably  lower  than  that  of  i  ground  surface.  For  ©saspltB 
it  U/o  =  2  @ad  k  =  0,5*  the  maximum  shear  stsss§  «.  l.S  x  10® 
psi*  which  is  of  the  same  order  of  magnitude  of  maximum  ligrtsian 
shear  gtr©§sa  Thus  severe  plastic  deformation  is  not  expected 
in  the  asperity* 

It  is  i©@u  (by  using  Figure  19)  that  still  lower  maximum 
sheer  stress  sill  occur  ia  lapped  surface,  @.gB  at  h/er  =  2,0, 
the  mtxi'sum  sheer  stress  in  this  case  is  only  3  x  10^  p@i 
(which  is  lower  than  the  yield  strength  of  hard  steel)  and  no 
plastic  deformation  (or  surface  distress)  should  occur. 
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Figure  21,  Scheaatie  Hepre sent ati on  of  an 
Idealized  Surface  Defect 
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SECTION  IK 


DETERMINATION  OF  SHEAR  STRESS  BENEATH  A  FUI10S 


A  t.  «o  dieensienal,  exploratory  §F>al  ysl  s  ax  the  cent  act 
8 tress  of  bodies  containing  an  iaeeliged  surface  mi cr ©-defect 
S3  show)  in  Figure  21  has  bean  conducted. 

The  dsfect  in  Figure  2l  is  a  two  dimensional  "depre  ss  1  on ", 

1,0,  infinitely  long  and  Figure  2i  shows  its  cross-section.  the 
rias  ef  which  are  formed  by  two  r&dii,  tangential  to  the  surround¬ 
ing  (original)  surface  at  points  2c-  apart.  The  Herta  area  is  *4a 
wide  where  a  >  c,  and  the  contact  does  not  extend  closer  t©  the 
bet  tea  of  the  defect  but  rather,  there  is  a  free  surface  of 
width  2b  in  the  defect,  whereby  e  >  b  >  o,  i.o.  the  free  sur¬ 
face  et  the  bottom  of  the  defect  is  free,  its  shape  ia  irrelevant 
provided  that  It  is  sufficiently  depressed  not  to  contact  the 
opposite  body. 

In  addition,  the  profile  is  assumed  rraooth,  having  radii 
of  curvature  at  all  points  considerably  greater  than  the  char¬ 
acteristic  dimensions  of  the  defect.  The  contact  region  consists 
of  two  portions  due  to  the  presence  of  the  open  cavity,  covering 
the  cres a-sec t 1 onu 1  co-ordinates  b  <  x  <  a  and  -a  <  x  <  -b 
where  the  values  of  a  and  b  (<c)  are  determined  by  the  defect 
geometry  parameters,  i.e.  r  and  c,  and  by  the  load. 

This  contact  problem  can  be  solved  by  applying  the 
Muskheli shvi i i  theory  of  complex  variables  to  the  mixed  boundary 
value  problem  of  an  elastic  half  plane  with  multi-contact  zones 
(36). 

In  rolling  elements,  the  surface  defects  are  small  such  that 
the  characteristic  dimensions,  c  and  r,  of  the  defects  are  of  a 
smeller  order  then  the  linear  dimensions  of  the  rolling  eloaents. 
Therefore,  it  is  Justified  to  consider  a  limiting  case  wherein 
c/H^o  and  r/R^'O,  Using  this  assumption,  the  problem  reduces  to 
tho  compression  of  two  straight-edged  bodies,  one  of  which  con¬ 
tains  a  shallow  surface  defect,  as  shown  in  Figure  22.  Frictions! 
traction  on  the  surface  is  neglected.  The  contact  pressure 
at  the  interface  when  the  defect  is  absent  is  assumed  to  be  pa„ 

In  the  presence  of  0  defect,  the  contact  pressure  at  points 
removed  from  tho  defect  is  expected  to  approach  asymptotically 
tha  un distorted  value  p0  as  x  approaches  infinity, 
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k  olostd  fora  gelation  lag  boon  derived  in  IppggdlK  l¥  lor 
Ills  Hailing  ease,  XI©  expression  for  %  Is «  eeatgot  prs§sar@  lg 
as  feilegss 


p(s) 


(E’/2TTsr )  [  Ka)  -  1  C-x> 


<?,  n 


I(s) 


( x-c ) • log 


(  x~b )  ■  Tl~ 
( x-b)  *  TW 


aad  E ' 


(  1  =  V4 


)  /ff  g  j  4-  (  1  -=■ 


the  reduaed  Yuuag's  aodultis 


(For  §l©tli9  I!  =  51.8  ife/is8) 

1)  ^  tan  (0. 5  •  cos  «  ) 

c 

fb©  tain©  of  b  is  deie raised  fros  ike  followiag  forsaulet 

leg  (\/e 1 - 1  ■>  <sj  )  =  C0  (9.2) 

«.=  e/b 


iQd  £g  a  r?spe  r/cg'  «  a  diaeasionloii  psraaeter. 

figare  13  plots  the  r©laii  oasiiip  be6@@ea  m  sad  C@0 
gist©©  Ge  oii  b©  computed  free  tie  known  ’fgri tbl©s  Pa  ©ad  r/ae 
it  is  peiilbl©  t©  abliig  b  «  e/m  wies©  ©is  @bt©ia©d  for  e  gires 
^sla©  sf  6@  froa  Figurt  13. 

A  as@@rioi,l  esaapli  i§§  beta  oespuied  for  fist  distributive 
of  eoateet  pressure  ©a  the  surface,  sisasiag  s /h  ”  1.209  ©orr@s= 
pegdfitg  %@  g^  =  0,069.  wslues  egg  b©  ofetaistd  by  setting 

(1)  r/e  =  l.g  eiiSgS  3x10^  pai  or  (2)  r/e  3  o,9  sad  p@  “  4  s  10®p§ie 
Figart  14  plots  tfe©  diaossi ©siiggs  pressure  p/p®  as  &  fsaetioa  os  x. 
Ik©  results  sLas  that  there  is  e  pressure  rise  in  tie  vioiaity  of 
the  defeet  ©dye,  resebiag  e  ¥alu@  of  3.0&,  ifie=  there  is  a  sigal~- 
fioaat  oeaeenirati ©a  of  pressure  at  the  defect  edge. 
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Figure  24,  Pressure  Distribution  Near  s  Surface  Defect 
(far  c/b  =  1,2) 
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Figure  25  shows  the  distance  x0  of  the  pressure  peak  arising 
a  s  «  result  of  the  defect  as  measured  front  the  defect  center  snfl 
the  magnitude  of  the  maximum  contact  pressure  in  terms  o*  p0 
plotted  as  a  function  of  c/b„  The  result  shows  the  high  pressure 
peaks  occur  within  the  area  of  radius  v.  Based  on  Figures  24  and 
25,  Figure  26  plots  the  variation  of  PmaH  at  a  function  of  r/c  ^ 
for  steel,  using  four  values  of  P0<i.e.,  1  x  2  x  10Je  3  x  10^ 

and  4  x  10®  psi  as  parofactcrs.  it  can  be  readily  seen  that  the 
maximum  contact  pressure  Increases  with  decreasing  values  of  r/c. 


constant 


it  can  be  seen  from  Figure  2t>  that  nil  the  curves  approach 
nt  values  of  pmax , /p0(:“2.  U )  when  r/c  is  greater  than  2.  It 


ix  -3  .  o  t.  « 'i  •  ^  i  ['max  .  ‘  v  o  '  -  •  *■'  '  " 11  ^  u  ~ 

can  be  concluded  that  for  all  values  of  r  relatively  small  com¬ 
pared  with  the  rolling  element  size  and  Po  <4  x  10J  psi,  wo 
have  pMOV  /p,>2.5. 


^-p0  B  i  x  IQ0  lb/in§ 

_-p  3  2  slOS  lb/ln8 
. o 


=  3  s  10°  lb/ Itp 


Pn»  4  8  10s  lb/ln* 


Defect  Corner  Radius  /  Semi  Defect  Width  r/V 


Figure  26.  Variation  of  Maximum  Contsct  Pressure  as  a  function 
of  Defect  Geometric  Parameter  r/c  and  Nominal 
Pressure  p„ 
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SUBS  US  FACE  STRESS  f) ISTBIBUTIGM  IN  THE  VICINITY  OF  A  SURFACE 

DEFECT 


i 


It  is  of  interest  to  know  the  stresses  existing  near  the 
high  pressure  posk  nt  x  £  Q,  The  sub-surface  stress  distribu- 
tion  eon  fee  computed  using  a  numerical  integration  technique 
bnged  oa  tho  solution  for  the  stress  field  on  a  half  plane  under 
a  concentrated  normal  load.  The  numerical  method  requires  that 
the  region  of  surface  loading  be  finite  in  width  along  the  X 
axis.  This  can  be  arranged  by  resolving  the  surface  pressure 
into  one  component  h  dch  occupies  a  finite  width  and  §  uniform 
pressure  noting  on  the  surface  of  the  entire  half  plane  as 
shown  in  Figure  27. 

Figure  20  plots  the  contours  of  equal  von  Mines  yield 
stress  6(3  for  a  typical  surface  defect  with  parameters  e/b  -  1,2 
corresponding  to  r/c  =  0,9  for  p0  ~  4  x  10^  psi  or  r/c  =  1,8 
for  p0  ~  2  x  105  psi  in  steel.  It  can  be  seen  that  the  maximum 
value  Toct0r  octahedral  shear  stress  occurs  under  the  location 
of  maximum  surface  pressure.  Thus  (  0q  )max.  1®  about  1.36  p0 
which  is  considerably  higher  than  0.7  pQ  occurring  at  the  axis 
of  iyaastry  (x  =  o)  shown  in  Figure  28.  Furthermore,  as  shown 
in  Figure  28e  the  depth  of  (a q  )mr}S  is  *.  0,13  e  which  is  con¬ 
siderably  smaller  than  0,9  c  which  occurs  along  x  -  0. 

Sn  the  above  defect  originated  stress  analysis,  the  deter- 
(aiding  parameters  are  found  to  be  the  ratio  of  defect  width  2c 
and  corner  radius  r  and  the  undisturbed  surface  pressure  p0„ 

Tht  degree  of  stress  concentration  at  the  shoulder  increases 
with  increasing  defect  width  and  decreasing  defect  shoulder 
radiug,  Although  subsurface  maximum  shear  stress  has  not  been 
computed  for  »  wide  range  of  $■  values,  on  example  has  been  corn- 
puled,  corresponding  to  a  typical  realistic  defect  size.  The 
result  shows  that  th '  contact  pressure  reaches  e  peak  value  of 
3.1  times  p  sheresg  the  maximum  value  of  the  quantity  giving  the 
vcn^Miseg  yield  criterion  is  1.3  pg  (compared  with  0,32  p0  in 
tho  Hertzian  contact.)  It  is  expected  that  significant  plastic 
deformation  will  occur  at  the  defect  shoulder. 

Surface  plastic  deformation  has  been  observed  at  asperities 
(on  ground  end  at  times  on  honed  surfaces)  st  low  %/e  vale©®,  and 
at  @eny  surface  defect  shoulders,  It  is  believed  that  this  ig 
due  to  the  high  shear  stregg  predicted  above  for  both  of  these 
surface  failure  origination  point®.  In  the  cast  of  severe 
asperity  interaction  the  plastic  deformation  occurs  at  tha  tip 
0  aasiay  a  dee  re as©  of  asperity  slope.  In  the  case  of  surface 
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defects,  sharp  corner  radii  (compared  with  defect  width)  have  been 
observed  in  new  surfaces.  After  many  cycles  of  rolling,  i  .  h  «  s 
•'Sen  observed  by  = efface  tracing  that,  the  corner  rat.  rs  increases 
considerably  as  n  result  of  "rounding  off*. 

It  is  plausible  to  assume  that  the  surface  profile  of  as¬ 
perities  and  defect  shoulders  stabilizes,  i .e,  the  plastic  defor¬ 
mation  ceases  to  grow  after  a  certain  number  of  str-ss  cycles. 

This  is  a  kind  of  "shakedown”  process  on  the  microscopic  scale. 
Elo’redge  and  Tabor  (30)  and  Johnson  (32)  in  studying  macroscopic 
shakedown  of  ball  tracks  hsv«'  concluded  that  after  e  certain 
shakedown  is  reached  the  material  will  behave  elastically.  In 
spite  of  the  large  difference  in  scale,  the  basic  mechanism 
involving  shakedown  in  bearing  rolling  tracks  d/or  asperity 
tips  (or  defect  shoulders)  can  be  similar.  Using  this  argument 
it  may  be  assumed  that  the  asperity  tip  and  defect  shoulder  will 
behave  elastically  after  a  stabilized  surface  is  achieved.  Based 
on  Figures  17  -  19  and  25  it  is  seen  that  a  run-in  asperity  (or 
defect  surface)  having  acquired  a  smaller  asperity  slope  (or  a 
larger  defeat  corner  radius)  suffers  a  lower  degree  of  stress 
concentration  than  when  it  is  new.  From  this  fact  it  can  be  said 
that  surface  plastic  deformation  flows  in  the  direction  required 
to  reduce  the  degree  of  stress  concentration.  In  hard  steel,  work 
hardening  is  high  and  the  amount  of  deformation  at  *s>  nakedown " 
will  be  limited.  It  is  possible  that  the  iun~in  shape  of  asperity 
tips  and  defect  shoulders  will  continue  to  have  some  stress  con¬ 
centration. 
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APPENDIX  I 


FORMULAS  FOR  STRESSES  IN  A  HERTZIAN  STRESS  MELD 


The  stresses  at  an  arbitrary  point  below  the  surface  of  a  bearing  ring  induced  oy 
contact  of  a  ball  or  roller  arc  computed  under  the  assumption  th  it  areas  of  the  ring  and 
rolling  body  are  large  enough  compared  to  the  size  of  the  contacting  area  that  the  con¬ 
tacting  bodies  may  be  considered  infinite  in  extent. 

In  the  vicinity  of  the  contact  the  surfaces  are  assumed  to  be  dec  crib  able  by  second 
degree  polynomials. 

All  of  the  stresses  concerned  arc  referred  to  a  rectangular  Cartesian  coordinate 
system  with  the  xy-plane  fixed  on  the  t  oundary  surface  of  the  eomi-infinite  body  with 
the  z  -  axis  directed  into  the  body.  The  xz  and  yz  planes  coinage  with  the  symmetry 
pianos  of  the  contact  ellipse,  the  equation  of  which  is, 


2  2 


where  a  is  the  major  semi-axis  lying  on  the  x-axis  and  b  the  minor  semi-axis  lying 
on  the  y  axis. 

The  stresses  are  given  by  (34). 


ZM  ZM 

+  (1-2  V)  Nx-2  (l-y)~  +  2V  ~~ 


ZM  ZM 

(l-2tf)  N  -2  (l-v)  +  2v-~- 

y  L  L 


(Al-2) 


(Al-3) 


(Al-4) 


(1-2  V  )  N 


(Al-5) 
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xz  ^  ZX 
— —  y  ~~7T~~~, 7 
rr  » ,  * 

o  A  *  L 


1  +  1j 


in  which, 

A=f  X-i  y4  z4 

b ,  b,  b  ,  b 

-3P 

ao  2  TT  ab 

and  P  is  the  total  toad 

L  is  the  largest  positive  root  of  the  following  equation 


r 


Jii 


arctan 


Y  V  A  -  i 


l.LK  j-J^F)  (1  +  L")_ 


= 


a2-i 


arctanh  p 


:A2 


2  2  Z  /  2~  2  2 

A  +L  +f-y(A  +L  )  (1  +  L  ) 


N  =  -~  (X@~Yf) 
A  -1 


(A  l -6) 


(Al-7) 


(Al-8a) 


(Al  -8b) 


.2  „2  ?2 

(Al-8c) 

Z  A 

l 

LV(A^L2)  (1  +  L2) 

(  \2  ( 

LY  N2  /Z\2 

\a2  +  l2/  V 

1  +  L  /  U) 

(Al-Gd) 

(Al-8e) 

(Al-8f) 

(Al-8g) 


nn 
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M 

x 


L  (F  -  E) 
2 

A  -1 


M  - 

y 


L  (A2  E  -  F) 


L2  A 


/a2+l2)  (1+:.2) 


(Al  -8J) 


(Al-8k) 


{Al -81) 


where  F  and  E  are  the  ordinary  elliptic  integrals  of  the  first,  and  second  kind,  re- 

A 

8  =  arctan  —  ;  V  Is 

Poisson's  ratio,  assumed  to  be  0.3  in  all  of  the  numerical  calculations. 

For  points  in  the  contact,  i.  e.  Z  =  0  ,  the  largest  positive  root  L  is  zero.  It 
follows  from  equation  (Al-8c)  that 


speotively,  with  modulus  k 


-J'-h 


and  argument 


iim  /  X2  y.2 

Z-o  L  \/*  "  .2 


(Al-9) 


Using  this  relation,  the  stress  formulas  can  easily  be  obtained  from  equations  (Al-2) 
through  (A!  =7)  as  follows: 


ax 

—  =  (1-2W)  N 
°o  X 


(1-2  p)  N 


(Al-10) 


APPENDIX  II 


FORMULAS  FOR  STRESSES  CORRESPONDING  TO  AN  INFINITELY 
NARROW  CONTACT  ELLIPSE  IN  A  HERTSAN  STRESS  FIELD 


In  the  limiting;  ease,  A-5  w  or  A  ~  =  9,  the  formulas  in  Appendix  I  must  be  modi¬ 
fied  before  they  can  be  used  in  numerical  calculations.  The  following  approjdraationa 
must  be  made  In  deriving  the  stress  formulas. 

It  is  proper  that,  before  working  on  the  at  rows  formulae,  attention  be  directed 
to  the  modifications  of  equations  (Al=8).  Designating 

„>  X 

x  —  A  ■ 

A  a 


Equation  (A1  -Be)  becomea 


1 


(A2-1) 


Tuus  L  is  the  largest  positive  root  satisfying  equation  (A2-1).  Sines 


lim 

A-® 


JFF? 


=  1 


it  is  evident  that 


Eliminating  Y  by  means  of  equation  (A2-1),  this  can  further  bo  reduced  to 


(A2-2a) 


As  A  -  ®,  both  $  and  $  approach  zero.  Thus 

N  0  =X' tf  =  0  (A2-2b) 
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N  0 

X  A 


(A2-  2c) 


N  ^  _Z,  1  .  X.  ^ 

y_  L  ‘/TT?  '  ~L 


(A2-2d> 


For  the  elliptic  integrals  ths  modulus  k  will  approach  unity  while  the  argument 
(9  will  be  close  to  r'/2.  Therefore,  it  ears  ba  Been  that  (40) 

lim  U|A^  =  Urn  ki2F(^/k)s  Jf1  k'2K 

A-4®  „2  k'-^O  '  ?  k'-’O 

A 

.  lim  .,2,  4 

«  .  ,  n  k'  la-r- 
k’ -^O  k' 


.  lim  .  ,2  2  n 

*  ~k.-0  k  lnk  =0 


and  lim  E  =  1,  where  k’  =Jl  -  k  . 
It  then  follows 


(A2-2©) 


(A2-2f) 


a  +  L  -  L 


<A2~2g) 


Thus  ih&  stress  formulas  can  be  obtained  directly  from  equations  (Al-2)  through 
(Al-7).  Thseo  are 


jk  2  yZ 


ij  z 


(A2~3) 
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PLANE  CONTACT  OF  ASPEMTIEf 


The  governing  equation  relating  the  contact  pressure  p(x)  in  the  region  -&<  x<  a 
to  a  symmetrical  surface  profile  designated  by  £(x),  is  given  by  (36). 


-a  z-t 


dt  -  -  £'  (x) 


(A3-1) 


where  E' 


[(l-J'^rrE1 


4-(l-y2ynE2 


f‘(x)  =  df{x)/ds 


Equation  (A3-1)  shows  that  the  plane  contact  problem  of  two  bodies  having  moduli!  E^ 
and  E„,  respectively,  is  equivalent  to  that  in  which  on©  of  the  Bolide  is  rigid 

(E  -«)  and  the  other  has  a  reduced  Young's  modulus,  E1  (l-v  )  /»  ,  provided  that  f(s) 
remains  the  same. 

The  above  singular  integral  equation  can  be  integrated  for  p(x).  For  the  special 
case  that  p  (+  a)  =  0,  i,  e.  the  contact  pressure  is  aero  at  the  contact  edges  (this  is  valid 
if  the  profile  does  not  have  a  sharp  earner  or  ridge  at  the  edges),  the  contact  pressure 
is  given  by  (36) 


E'  /2  2 

PCs)  -rva  -x 


a  J&  -t  (t-x) 


For  the  particular  profile  depicted  in  Figure  12,  one  has 

t'  (t)  =  t/B  fori  1 1  c  b  ; 

=  b/R  for  a  >|tl 


(A3 -2) 


(A3  “3) 


Substituting  equation  (A3-3)  ir.to  equation  (A3=2)  and  carrying  out  the  integration, 
the  pressure  distribution  is  obtained  as  follows:*  \  /  A 


is*  r*  /  2  2 

pfe)  =  —  T-  2  /a  -x  0Q  -  (x+b)  log 


2rr  R 


JO 


-  (x-b)  log 


mv, 


-G£ 

1  O 


(A3 -4) 


wraniuimib 
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APPENDIX  IV 


COMPRESSION  OF1  A  HALF  PLANE  CO  NT  AIM  NG  AN  IDEALIZED  SURFACE  DEFECT 


As  depicted  in  Figure  22,  the  lielf  plane  is  straight  edged  ©seep*  at  the  dapr^isloa. 
The  line  of  symmetry  of  the  profile  is  chosen  as  the  y-axis  whereas  the  edge  of  the 
half  plane  is  chosen  as  the  s-asie.  The  comer  radii  of  the  depression,  assumed  to  be 
identical,  are  denoted  by  r.  Tho  approach  of  the  two  bodies  is  assumed  to  be  normal 
to  the  surface  of  tho  half  plane.  Thus  at  infinity  there  is  an  uniform  compression  pQ 
throughout  the  bodies.  The  contact  region  occupies  the  boundary  of  the  half  plane 
except  at  the  center  of  the  depression  or  !  xt  <b. 


Using  the  same  notations  m  in  Appendix  in,  the  governing  integral  equation, 
similar  to  equation  (A3=l)  in  form  exempt  for  the  limit  of  the  integral,  is  given  below 


JL 

E' 


t~  f  “b 


b 


P<*) 
J  x  -t 


dt  =  -  S'  (x) 


(A4~l) 


The  derivative  of  tbs  profile  expression  £  (x)  with  respect  to  x  for  the  present  problem 
is  given  by: 

f'(x)  =0  e>|xl>e 

£'  fe)  =  (e-s)/r  b  <2  <  c 

=  (-C"X)/r  -c  <  x  <  -b  (A4-2) 

According  to  (30)  the  complex  potential  in  Ms  problem,  for  the  case  p  (  +  b)  o, 
can  be  given  as  follows: 


-  ,  .  E'  L  2.2 
#(g)  =  — -S  yjQ  -b 

4tt 


~b  f  ® 

•  a 

T  I 

-»  b 


B  t  Q-K+iz 

t  =b  f i-C)  5 


(A4-3) 


The  contact  pressure  in  the  contact  repon  is  given  by  (36). 

e>  ~b  />  e  . 
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2rr 


2 

8  -b 


£M  <te_ 


t  -b  (t“X) 


(A 4-4: 
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mwU'mung  equation  (A4-2)  into  equation  (A4=4)  and  rear 


ranging  the  teens,  yields: 


P  (k)  = 


2  2  x 

s  "  b  •  '■== 

r 


t’=b2  (t2  -*2) 


2tr2? 


Tl  (x)  -  I  hcj] 


■wbsia  I  Cs)  =  (s-c)  lo| 


te-b )T) 

C?-b)T|  ,,//  -b2 


•  tea  [o.  5  ooa  -1  fb/e)  J 


(A4-5) 


nation  (A4-5)  ceB  be  ^2,  ^  . 

fU=^r  4  £  logfv^Tj  ^  _  /?  7  b  “1 


Po  2  r 
TT 


1  Vi 


(A4-6) 


uUli)3  p0,  e,  r  and  £  are  known  constant®,  it  is  Doasihia  tn  A ,  ,  . 

required  in  equation  (A4-5)  The  d&i&miina+in  /»,  ~  to  so-*v@  i0r  b  which  is 

na — on  of  e/b  can  be  made  by  using  Figure  2s 

P-uisiag  e/b  8f atnst  the  dixnensional  quantity  c  =  tv  2p  r/c  E 1 


CTDuumi  lUEDuiiiiciwnnuidg  iiia’iiMinjViia.'Mu 
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APPEN.uiX  V 

AVERAGE  SHEAR  RANGE  IN  THE  STRESSED  AREA,  S,  ENCLOSED  BY  A  CONTOUR 

OF  EQUAL  SHEAR  RANGE,  T_  ,  IN  A  HERTZIAN  ELLIPTICAL  CONTACT 

H 


The  relationship  between  the  area,  B,  enclosed  by  a  contour  of  equal  shear  range 
T  anti  Tie  given  in  equation  (7-5)  of  Section  VII  os  follows: 

rt'T'  3/4 


r  2T  - T  1 
o 

o  L  T 


S-  u  a  2 


(A5-1) 


in  which  a  -  the  major  axle  of  the  contact  ellipse,  z  =  the  depth  corresponding  to 
mariuium  reversing  shear  stress  T  and  c=  constant. 

°  =1 

From  (A4-1)  T  =  2  ?o  (cazo)  4/3  S4/<3  -k  (cazQ)  4/3  J 

’5  '  ""  °  (A5-2) 

For  a  specific  stressed  area  6  bounded  by  a  contour  of  equal  shear  range  the 
averagg  shear  stress  becomes 


(1E^av.  ”  S 


T  An  s 


2To  fS 


,(cazo) 

■37577“  4/3 
b  +  (cazQ) 


2T  eaz 
0  o 


3/cas 


w4/^3  + 1 


(A5-3) 


The  upper  limit  of  ths  integral  can  be  determined  from  equation  (Ah-1)  by  moans 
of  Figure  10,  and  therefore  the  integral  can  be  evaluated  using  a  numerical 
integral! oa  technique. 

Using  equation  (A5-3),  values  of — ’  have  been  calculated  and  plotted  in 

(Figure  SB  )  against  B/§5q  „  jt  can  be  sltowu  thee  for  ‘r^mn  smaller  than  4. 7,  the 
curve  can  be  approximated  by  a  sin14  gut  line  expressed  as  follows: 


a  <TR>av 

—  =  -  6.  48  +  12.  9 

as  T 

o  o 

Solving  for  (T  )  in  equation  (A5-4),  one  has 
!TS>,v./T0  >0-1“(“.»'8/«o) 


(A5-4) 


(A5-5) 
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Figure  32.  Variation  of  Average  5ft 
Range  in  a  Contour  of  E 
To  with  Enclosed  Area  S 


APPENDIX  VI 


plausbJjE  defect  sevekty  merniDUTiOMS 


In  what  follows  we  will  explore  the  behavior  of  the  function  T  necessary  for  the 
satisfaction  of  Equation  (6=9)  for  two  plausible  ehoie-  §  of  F(d). 

It  may  be  reasoned  that  the  distribution  of  defect  severity  is  such  that  most  of  the 
defects  have-very  small  severity  and  a  diminishing  pror-or  Uou  have  greater  severities. 
That  is,  the  proportion  of  defects  with  severity  in  the  interval  +  A  d  is  larger  than 
the  proportion  in  the  interval  as  long  as 

The  density  of  the  two  parameter  exponential  distribution  with  unit  location  para¬ 
meter  to  be  given  below  has  mode  at  unity  and  decreases  with  increasing  variate 
values  and  thus  satisfies  the  condition  poatulated  above. 

With  this  law  the  probability  that  a  randomly  selected  defect  has  severity  less 
than  a  value  d  is  ^ven  by 


d  <  d  = 


=  F(d)  •=  1  -  exp 


-(t)1 


=  0  ;  d  <  ] 


(A6-1) 


The  quantity  d„  is  a  constant  parameter  of  the  distribution  related  to  the  average 
severity  d  by  d0  =  d  -1.  d  may  vary  with  different  materials. 

Another  distrihutioii  that  satisfies  the  above  postulated  condition  is 


F{d)  =  1  -  d"C  ;  d  &  1,  c  >  0 

=  0  ;  d  <  1  (A6-2) 

where  c  is  a  constant  parameter. 

Equation  (A6-2)  is  a  form  of  Pareto's  distribution  (47),  and  has  boon  found  useful 
in  economic  studies.  The  parameter  c  is  expressible  in  terms  of  are  average  defect 
severity  5  as 

d 

c  = 

d  -1  (A6-3) 

The  distributions  of  Equations  (AG-1)  and  (AG  =2)  are  shown  ou  Figure  30,  plotted 
for  d  =  2.0.  The  corresponding  density  functions  are  shown  on  Figure  31. 


PROBABILITY  DENSITY  RJMCTIOM  OF  DEFECT  SEVERITY 


Figure  31,  Probability  Density  Function  of  Defect 
Severity  with  d  -  2,0 
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1 .  B1HAVIOH  OF  F  AS  N 


Using  the  espoaaniigl  distribution  of  Equation  (Ag-i)  as  the  defect  life  distribution 
in  EqugQoa  (6-1)  gives 


G  {Kj )  3  esp  '■  g 


Invoking  the  requirement  of  Equation  (S-9)  gives 
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TlaJUf)  r°-  amt  feahuve 


1  /  B  \  a 


as  N  -*  H0«  0 


Taking  th©  iaverse 
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(A8=8) 


Equation  |A6-6)  states  that  as  d  b©eoiaes  large*  F  (d)  behaves  as  an  gspoasntial 
fcaeiioa.  Th©  result  is  required  osly  Bsymptoilavily  gsd  sgasot  te  true  over  fee 
whole  range  of  d  sicse  according  to  Equation  (A6-8)  T  does  mst  approach  sero  as 
d  -5 1  bs  reqtdted  by  %isilea  (§-12)  $3  -  1  corresponds  to  bo  severity). 


<A6-U) 
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